Artificial forms of mechanical limb stimulation are used within multiple fields of study to determine the level of cortical excitability and to map the trajectory of neuronal recovery from cortical damage or disease. Square-wave mechanical or electrical stimuli are often used in these studies, but a characterization of sensory-evoked response properties to square-waves with distinct fundamental frequencies but overlapping harmonics has not been performed. To distinguish between somatic stimuli, the primary somatosensory cortex must be able to represent distinct stimuli with unique patterns of activity, even if they have overlapping features. Thus, mechanical square-wave stimulation was used in conjunction with regional and cellular imaging to examine regional and cellular response properties evoked by different frequencies of stimulation. Flavoprotein autofluorescence imaging was used to map the somatosensory cortex of anaesthetized C57BL/6 mice, and in vivo two-photon Ca2+ imaging was used to define patterns of neuronal activation during mechanical square-wave stimulation of the contralateral forelimb or hindlimb at various frequencies (3, 10, 100, 200, and 300 Hz). The data revealed that neurons within the limb associated somatosensory cortex responding to various frequencies of square-wave stimuli exhibit stimulus-specific patterns of activity. Subsets of neurons were found to have sensory-evoked activity that is either primarily responsive to single stimulus frequencies or broadly responsive to multiple frequencies of limb stimulation. High frequency stimuli were shown to elicit more population activity, with a greater percentage of the population responding and greater percentage of cells with high amplitude responses. Stimulus-evoked cell-cell correlations within these neuronal networks varied as a function of frequency of stimulation, such that each stimulus elicited a distinct pattern that was more consistent across multiple trials of the same stimulus compared to trials at different frequencies of stimulation. The variation in cortical response to different square-wave stimuli can thus be represented by the population pattern of supra-threshold Ca2+ transients, the magnitude and temporal properties of the evoked activity, and the structure of the stimulus-evoked correlation between neurons.