This paper presents a novel method to quantify spatial changes in muscle activation pattern by multi-channel surface electromyography (MCSEMG) in order to investigate motor unit recruitment variation. The method is based on non-uniform distributions of motor units that cause spatial inhomogeneous muscle activation. To evaluate the method, 15 subjects performed three isometric elbow flexion contractions consisting of slow sinusoidal changes in force ranging from 0% to 80% of the maximal voluntary contraction. MCSEMG electrodes were placed in a 10 x 13 grid over the biceps brachii muscle. From all channels, root mean square (RMS) values of bipolar leadings were computed over 0.5 s epochs over the whole recording. Thereafter, correlation coefficients were calculated between the RMS values at one epoch, with the RMS values at another epoch. Results showed consistent spatial changes in the distribution of RMS at different contraction levels up to 80% of maximal voluntary contraction and when comparing increasing and decreasing contractions at the same force level. These findings are reproducible within and between subjects, and in agreement with physiological phenomena and therefore indicate that the spatial inhomogeneities of motor unit properties in the biceps brachii muscle can be used to study changes in motor unit recruitment.