Journal of Applied Physiology Ad Instruments
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 98: 120-131, 2005. First published September 17, 2004; doi:10.1152/japplphysiol.00894.2004
8750-7587/05 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
98/1/120    most recent
00894.2004v1
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (85)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Keenan, K. G.
Right arrow Articles by Enoka, R. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Keenan, K. G.
Right arrow Articles by Enoka, R. M.

Influence of amplitude cancellation on the simulated surface electromyogram

Kevin G. Keenan,1 Dario Farina,2,3 Katrina S. Maluf,1 Roberto Merletti,2 and Roger M. Enoka1

1Department of Integrative Physiology, University of Colorado at Boulder, Boulder, Colorado; 2Laboratorio di Ingegneria del Sistema Neuromuscolare, Dipartimento di Elettronica, Politecnico di Torino, Torino, Italy; and 3Center for Sensory-Motor Interaction (SMI), Department of Health Science and Technology, Aalborg University, Aalborg, Denmark

Submitted 18 August 2004 ; accepted in final form 14 September 2004


    ABSTRACT
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
The purpose of the study was to quantify the influence of selected motor unit properties and patterns of activity on amplitude cancellation in the simulated surface electromyogram (EMG). The study involved computer simulations of a motor unit population with physiologically defined recruitment and rate coding characteristics that activated muscle fibers whose potentials were recorded on the skin over the muscle. Amplitude cancellation was quantified as the percent difference in signal amplitude when motor unit potentials were summed before and after rectification. The simulations involved varying the level of activation for the motor unit population, the recording configuration, the upper limit of motor unit recruitment, peak discharge rates, the amount of motor unit synchronization, muscle fiber length, the thickness of the subcutaneous tissue, and the motor unit properties that change with advancing age. The results confirmed a previous experimental report (Day SJ and Hulliger M, J Neurophysiol 86: 2144–2158, 2001) that amplitude cancellation in the surface EMG can reach 62% at maximal activation. A decrease in the range of amplitudes of the motor unit potentials, as can occur during fatiguing contractions, increased amplitude cancellation up to ~85%. Differences in the amount of amplitude cancellation were observed across all simulated conditions, and resulted in substantial changes in the absolute magnitude of the EMG signal. The most profound factors influencing amplitude cancellation were the number of active motor units and the duration of the action potentials. The effects of amplitude cancellation were minimal (<5%) when the EMG amplitude was normalized to maximal values, with the exception of variations in peak discharge rate and recruitment range, which resulted in differences up to 17% in the normalized EMG signal across conditions. These results indicate the amount of amplitude cancellation that can occur in various experimental conditions and its influence on absolute and relative measures of EMG amplitude.

computer simulations; peak discharge rate; fatigue; motor unit; normalization


THE SURFACE ELECTROMYOGRAM (EMG) is the algebraic sum of trains of action potentials generated by active motor units as detected by electrodes placed on the skin over a muscle. Because the surface EMG is influenced by the number of motor units that are active and the rates at which the motor neurons discharge action potentials, the EMG signal has been widely used to quantify motor unit activity and to estimate output from the spinal cord. However, the surface EMG underestimates the amount of motor unit activity due to the loss of information that occurs when overlapping positive and negative phases of motor unit potentials cancel one another and reduce the amplitude of the signal (5).

Many techniques are limited by amplitude cancellation in the EMG, including the decomposition of an interference EMG signal into its constituent motor unit potentials (24), the estimation of changes in central conduction from the amplitude of motor evoked potentials (52), and spike-triggered averaging of the interference EMG (19). Furthermore, amplitude cancellation likely confounds interpretation of changes in the surface EMG, such as assessing correlated activity in surface EMG signals from pairs of muscles (27, 35) and explaining the failure of EMG amplitude to reach maximal levels at the endurance limit during a submaximal fatiguing contraction (23).

Although problems and limitations due to amplitude cancellation have been recognized for several decades (1, 44, 49), only one previous experimental study quantified amplitude cancellation in the interference EMG (5). Groups of motor units were stimulated in a cat hindlimb muscle, and amplitude cancellation was quantified by computing the difference between EMG amplitude derived from rectified and unrectified trains of motor unit potentials. At maximal levels of imposed excitation, amplitude cancellation reduced EMG amplitude by 50%. Due to the simultaneous activation of many motor units, however, Day and Hulliger (5) suggested that the actual amount of amplitude cancellation might be greater during a voluntary contraction. Furthermore, to generalize the results of the experimental study it is necessary to determine the sensitivity of amplitude cancellation to variation in the relevant physiological parameters.

Computational models have proven useful in characterizing the sensitivity of the surface EMG to the parameters of the systems involved in its generation and detection (4, 15, 16), and the application of these models has been suggested as a strategy for understanding amplitude cancellation (55). The purpose of the study was to quantify the influence of selected motor unit properties and patterns of activity on amplitude cancellation in the simulated surface EMG. The expectation was that increases in the amount of overlap between action potentials, such as occurs with conditions that increase the numbers and durations of motor unit potentials, would enhance amplitude cancellation and confound the use of the surface EMG as a population index of motor unit activity. Some of these data were presented in abstract form (34).


    METHODS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
The study involved computer simulations that were based on models previously described in detail by Fuglevand and colleagues (21, 22), with the addition of a model of surface EMG that incorporated muscle fibers with a finite length and a volume conductor with layers of skin, subcutaneous, and muscle (15). These models represent the activation of a motor neuron pool from minimal to maximal levels of excitation, and they produce motor unit potentials and a surface-detected EMG consistent with the known physiological properties of first dorsal interosseus muscle. The simulations involved three main steps: 1) determination of the recruitment and discharge times of a population of 120 motor neurons in response to different levels of excitation; 2) simulation of motor unit potentials from estimates of the number and location of muscle fibers for each motor unit and the conduction velocities of the muscle fiber action potentials; and 3) simulation of the surface EMG by summing the trains of motor unit potentials. Amplitude cancellation was quantified as the percent change in signal amplitude when motor unit potentials were summed before and after rectification (5). The amplitude of the surface EMG and percent amplitude cancellation were compared across conditions that simulated different properties and patterns of activity of the motor neuron pool, as well as different recording configurations (Table 1).


View this table:
[in this window]
[in a new window]
 
Table 1. Description of simulation conditions

 
Surface EMG Simulation

The model was implemented in Matlab version 6.5 (The Mathworks, Natick, MA). The basic parameters in the model were similar to those published previously (56, 57). The distributions of properties across the motor unit pool were based on the Size Principle (29), and these included recruitment thresholds, innervation number, motor unit territory, and the conduction velocity of muscle fiber action potentials. Activation of the motor unit pool was modeled as a ramp-and-hold function, with a 1-s ramp increase in excitation to a mean level that was constant for 4 s. Input was uniformly distributed across the motor neuron pool and all neurons received the same level of excitation, thereby allowing simulated motor neuron input-output functions to be specified by well-established relations between discharge rate and injected current (28). Maximal excitation was denoted as the level of input necessary to bring the last recruited motor neuron to its assigned peak discharge rate and values of excitation were expressed as a percent of the calculated maximum. The distribution of recruitment thresholds for the motor neurons was determined from an exponential function with many low-threshold neurons and progressively fewer high-threshold neurons (12). Each motor unit began discharging at 8 pulses per second (pps) once excitation exceeded the assigned recruitment threshold of the unit, and discharge rates increased linearly with increased excitation (3 pps per 10% increase in excitation). As described by Fuglevand et al. (21), the first recruited unit (MU 1) had a maximal discharge rate of 35 pps, whereas peak discharge rate decreased linearly with increasing recruitment threshold, with the last recruited unit (MU 120) assigned a peak discharge rate of 25 pps. The discharge rate was modeled as a random process with a Gaussian distribution, varying the coefficient of variation for discharge rate from 10 to 40%.

Motor unit territories.   The simulated muscle had a circular cross-section with a radius of 8.67 mm, derived from physiological cross-sectional areas calculated by Keen et al. (33). The number of muscle fibers was 66,000, based on an average fiber diameter of 56 µm (8), a muscle radius of 8.67 mm, and an assumption that the noncontractile tissue accounted for 20% of the cross-sectional area. These values are similar to those used originally by Fuglevand et al. (21) (71,747 muscle fibers and a muscle radius of 7.5 mm). The number of fibers innervated by a single motor neuron and the cross-sectional area of the motor unit increased exponentially from MU 1 (26 muscle fibers, 1.3 mm2) to MU 120 (2,510 muscle fibers, 125.5 mm2). An exponential increase in motor unit territories reflected the skewed distribution of motor unit forces, with more motor units that exert small forces (42) and the high correlation between innervation number and tetanic force of a motor unit (32). The fibers of a motor unit were scattered over a broad region of the muscle cross-section and intermingled with fibers belonging to many other units. Motor unit territories were randomly distributed within the muscle and were circular, except when constrained by the muscle boundary (e.g., MU 120 in Fig. 1). The density of fibers within the territory of the motor unit was assumed to be ~20 fibers/mm2 but was increased when a portion of the motor unit territory was constrained by the muscle boundary. The result was little change in fiber density for small motor units, but a greater fiber density for the largest motor units, consistent with experimental findings (32).



View larger version (35K):
[in this window]
[in a new window]
 
Fig. 1. Schematic of the model used to generate surface-recorded electromyogram (EMG) signals. A: cross-section of the volume conductor model. Three tissue layers (skin, subcutaneous, and muscle) were simulated with the recording system positioned directly over the center of the muscle. Four of the 120 motor unit territories are shown, with each filled circle indicating an individual muscle fiber of a given motor unit. Innervation number increased exponentially from MU 1 to MU 120. B: lateral view of a motor unit with 4 fibers. Innervation zone was set at 40% of the length of the muscle, and the recording system was positioned halfway between the innervation zone and the distal tendon. Insertion of each fiber into the tendons varied randomly over a range of 5 mm.

 
Generation of motor unit potentials.   The surface EMG model simulated action potentials generated by muscle fibers with a finite length and detected by surface electrodes with physical dimensions. The volume conductor was an anisotropic and nonhomogeneous medium representing the muscle, subcutaneous (1.5 mm thick), and skin (1 mm thick) tissues. Each tissue layer was homogeneous and the conductivities of the tissue layers were the same as in Farina and Merletti (15): muscle was anisotropic and more conductive in the longitudinal fiber direction, whereas the subcutaneous and skin tissues were isotropic. The analytic description of the intracellular action potential in the spatial domain was based on Rosenfalck (51) and did not include the negative afterpotential. There was no electrical interaction between the activities of adjacent muscle fibers. The surface-recorded motor unit potential was determined as the sum of the contributions generated by the muscle fibers assigned to the motor unit. The conduction velocities of muscle fiber action potentials belonging to a motor unit were related to motor unit size, with values increasing exponentially from 2.5 m/s for MU 1 to 5.5 m/s for MU 120.

Simulation Conditions

Selected conditions were simulated to determine the influence of specific physiological and signal-detection parameters on EMG amplitude and amplitude cancellation. The model parameter values are shown in Table 1; the condition marked with asterisks is referred to as the default condition. Although the default parameters characterized the first dorsal interosseus muscle, the parameters of fiber length, subcutaneous tissue thickness, peak discharge rate, and recruitment range were varied to investigate the effect of changes within this muscle and across other muscles. For each condition, 11 levels of excitatory drive to the motor unit population were simulated: 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, and 100% maximal excitation. The EMG was sampled at a rate of 4,096 samples/s.

Recording configuration.   The simulated signals in the default condition were detected by using a bipolar configuration with an interelectrode distance of 10 mm and circular electrodes (4 mm diameter). The center of the recording configuration was halfway between the innervation zone and the distal tendon. The simulated bipolar recordings were compared with a belly-tendon configuration in which one electrode was placed over the belly of the muscle and the other near the metacarpophalangeal joint (53, 63). The belly-tendon configuration had one recording electrode halfway along the muscle and the other electrode was 32.4 mm away, past the distal tendon region.

Activation pattern.   Most motor units in first dorsal interosseus are recruited at forces <50% maximal voluntary contraction (MVC) force (42). In other muscles, such as biceps brachii, recruitment has been observed at forces up to 85% MVC (7, 38). It has been suggested that variation in recruitment range may contribute to the different EMG-force relations observed experimentally (3, 21). The possibility that amplitude cancellation might vary with recruitment range was examined by simulating narrow and broad recruitment ranges. The upper limit of recruitment was 41% of maximal excitation for the narrow range and 68% of maximal excitation for the broad range.

As excitation to the motor neuron pool increases, the discharge rates of motor neurons increase up to a level where little further change in discharge rate is observed. Peak discharge rates vary across muscles (2, 59) and are sometimes greater for low-threshold motor units (6, 30), whereas others reported that high-threshold motor units reach greater discharge rates (25, 43). The influence of an increase in peak discharge rate on amplitude cancellation was examined by increasing peak discharge rates from the default condition (35–25 pps for MU 1 and MU 120, respectively) to 55–45 pps. In addition, a condition where high-threshold motor units achieved greater peak rates was simulated (25–35 pps), as well as a condition where the peak discharge rate was the same (30 pps) for the entire motor unit pool.

The increase in EMG amplitude that occurs when motor units discharge action potentials at about the same time, a condition known as motor unit synchronization, might be caused by a reduction in amplitude cancellation (37, 61). To assess this possibility, motor unit synchronization was examined by adjusting the timing of selected action potentials to impose a temporal association with action potentials discharged by different motor neurons (61). A function was applied that selected between 20 and 30% of the action potentials discharged by each motor unit to be synchronized with approximately eight other motor units with similar recruitment thresholds. Motor unit synchronization was quantified with the common input strength (CIS) index as the frequency of extra synchronous discharges (48). This function fixed CIS values near a value of two across the motor unit pool at all levels of excitation. To determine the level of motor unit synchronization across the motor unit pool, CIS values were calculated for every motor unit with the 15 preceding units and the 15 succeeding motor units.

Fiber length and subcutaneous tissue thickness.   Fiber length and average location of the innervation zone in the first dorsal interosseus were determined experimentally using a linear electrode array consisting of 16-pin electrodes with 2.5 mm distance between electrodes. The bipolar recordings were amplified, band-pass filtered (10–500 Hz), digitized at 2,048 samples/s, and the resultant 15 channels were displayed on a computer monitor. The innervation zone was marked by visual inspection as the location where motor unit potentials began propagating in two directions. The tendon ending was denoted as the location where the propagation of the motor unit potential ceased (40). The Local Ethics Committee of the Health Department of Region Piemonte in Italy approved the measurements. On the basis of measurements from five men, average fiber length in the default condition was set at 40 mm and the center of the innervation zone was located 40% along the length of the fibers distal to the proximal attachment. Insertion of each fiber into the tendons varied randomly (uniform distribution) over a range of 5 mm, resulting in fiber lengths of 35–45 mm.

The properties of the tissues separating the muscle fibers and the recording electrodes are known to influence surface recordings (14, 50). An increase in subcutaneous tissue thickness attenuates the amplitude of the EMG signal and changes the frequency content of the signal (18, 50). To determine the default value for subcutaneous tissue thickness in the model, ultrasound recordings (FFsonic UF-4000L, Fukuda Denshi) were taken from the first dorsal interosseus of 12 individuals (11 men, age range: 22–59 yr), under approval of the Local Ethics Committee, with the probe (7.5 MHz transducer) positioned between the innervation zone and distal tendon. Thickness values ranged from 1.3 to 2.4 mm. On the basis of these recordings, the thickness of the subcutaneous tissue was set at 1.5 mm in the default condition.

Surface recordings of muscle-fiber action potentials are influenced by the propagation of the intracellular action potential along the fiber and its extinction at the end of the fiber. The end-of-fiber components influence the shape and the power spectrum of surface detected potentials (16). The relative influence of end-of-fiber components with respect to the propagating signal components depends on fiber length, thickness of the subcutaneous layers, and depth of the fiber (17). To evaluate the effect of surface-potential shape on amplitude cancellation, fiber lengths of 20 mm and 100 mm and subcutaneous tissue thickness of 1.0 and 3.5 mm were simulated and compared with the default condition (Table 1).

Motor unit properties.   To examine the combined influence of motor unit properties on amplitude cancellation, simulations included factors that are known to change with advancing age: a reduction in the number of motor units and of the innervated muscle fibers, increased density of motor unit fibers and increased motor unit territory, increased area of the innervation zone, increased variability of the diameter of muscle fibers within a motor unit, and decreased diameter of large muscle fibers (11, 31, 36, 39, 46, 58). Table 1 summarizes the values that were used to simulate advancing age.

Data Analysis

The dependent variables of this study were the amplitude of the simulated EMG and the percent amplitude cancellation as a function of excitation level. Average values for the rectified EMG (average EMG) were determined over the interval from 1 to 5 s at each level of excitation. Amplitude cancellation involves the loss of signal when unrectified action potentials are summed together and there is an overlap of positive and negative phases of the potentials. Cancellation does not occur when individual action potentials are rectified before summation (Fig. 2). Percent loss in amplitude due to cancellation was quantified at each excitation level by computing the difference between EMG amplitude derived from rectified (No Cancellation) and unrectified (Cancellation) trains of motor unit potentials, expressed relative to the EMG amplitude in the No-Cancellation condition (5).



View larger version (25K):
[in this window]
[in a new window]
 
Fig. 2. Demonstration of amplitude cancellation. A, top: 2 motor unit potentials were summed together and then rectified to determine an average EMG value. Bottom: rectified EMG signal at 50% excitation formed by summing motor unit potentials and rectifying the interference EMG. B, top: 2 motor unit potentials were rectified and then summed (not applicable to experimental EMG), ensuring that amplitude cancellation did not occur. Bottom: rectified EMG signal at 50% excitation was formed using the same motor unit potentials and discharge patterns as in A, but the potentials were summed after being rectified. au, Arbitrary units.

 
For each condition, 20 libraries of 120 motor unit potentials were simulated. Motor units were randomly assigned a position within the muscle in each library. EMG amplitude and percent amplitude cancellation were calculated as the mean of 20 libraries. The purpose for including 20 random motor unit locations was to assess the influence of location on EMG amplitude and percent amplitude cancellation. The simulations were insensitive to this parameter, however, and the standard deviations in percent amplitude cancellation were generally 2% across excitation levels. Therefore, only mean values are reported in the text.

The activation pattern was identical for all simulations with the exception of simulations for motor unit synchronization, recruitment range, peak discharge rate, and reduced motor unit number. Therefore, only the characteristics of the motor unit potentials changed between conditions. Four parameters describing the motor unit potential were calculated for each condition: peak-to-peak amplitude, area, duration, and normalized duration [area divided by amplitude (47)]. Correlation coefficients were calculated to determine the relation between these parameters of the motor unit potentials and the percent amplitude cancellation at 10, 50, and 100% excitation. Each parameter was calculated as the mean value for all 120 motor units across the 20 random motor unit locations.


    RESULTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Selected parameters of the default model were independently adjusted to examine their influence on amplitude cancellation. Additionally, the combined effect of age-associated changes in motor unit properties on amplitude cancellation was examined by adjusting several parameters concurrently.

Quantification of Amplitude Cancellation in the Surface EMG

The difference between average EMGs computed from the sum of the rectified potentials (No Cancellation) and unrectified potentials (Cancellation) at each excitation level indicated that the amount of amplitude cancellation was substantial (Fig. 3A). In the default model, percent amplitude cancellation increased up to 61.7% at 100% maximal excitation. The values across all excitation levels for the default model are reported in Table 2.



View larger version (13K):
[in this window]
[in a new window]
 
Fig. 3. Simulated relation between excitation and average EMG. A: average EMG was calculated by summing motor unit potentials before (Cancellation) and after (No Cancellation) rectification at 11 excitation levels. Difference between the Cancellation and No-Cancellation conditions at each excitation level indicates the amount of amplitude cancellation. B: Cancellation and No-Cancellation conditions for the default condition of the model were normalized to their respective maximal values. C: difference between the 2 normalized conditions is depicted at each excitation level, with the greatest difference occurring near the excitation level where motor unit recruitment was complete (40% maximal excitation).

 

View this table:
[in this window]
[in a new window]
 
Table 2. Percent amplitude cancellation across conditions

 
Due to the many factors that influence the amplitude of the EMG signal, it is common to normalize EMG to the amplitude of the EMG at maximal levels of excitation. The relation between excitation and EMG amplitude for the No-Cancellation condition was linear (r2 = 0.99, Fig. 3, A and B). However, the increase in EMG amplitude with increasing excitation for the Cancellation condition (Fig. 3, A and B) had a stronger relation for a sigmoidal function (r2 = 0.99) compared with a linear function (r2 = 0.95). The difference between the normalized EMG for the Cancellation and No-Cancellation conditions at each level of excitation (Fig. 3C) resulted in the normalized EMG values overestimating motor unit activity at intermediate levels of excitation.

The convexity in normalized EMG values at intermediate levels of excitation may be caused by the recruitment of progressively larger motor units and therefore motor unit potentials with greater amplitudes. The possibility that the range of potential amplitudes influenced the amount of cancellation was examined by simulating a condition where each action potential had the same amplitude. Specifically, the potential for MU 1, MU 60, or MU 120 in the default condition was assigned to all 120 motor units. Identical motor unit potentials resulted in the same nonlinear relation depicted in Fig. 3. However, the percent amplitude cancellation increased up to 88% for MU 1, 85% for MU 60, and 80% for MU 120 at maximal levels of excitation (Table 2). Thus decreased variability in the range of amplitudes of the motor unit potentials increased amplitude cancellation, but it did not contribute to the nonlinear relation between excitation and EMG amplitude.

Amplitude Estimation and Recording Configuration

The EMG signal is often recorded using either the bipolar or belly-tendon recording configuration, and amplitude is commonly estimated by calculating either the average EMG or the root mean square (RMS) value. The RMS values were larger (Fig. 4A) and amplitude cancellation was less (Fig. 4B) for the RMS estimates compared with average EMG, regardless of the recording configuration that was simulated. However, when the absolute EMG values were normalized to the value at maximal excitation, the largest difference between the normalized values was <5% at 40% excitation (Fig. 4C).



View larger version (18K):
[in this window]
[in a new window]
 
Fig. 4. Variation in EMG amplitude and percent amplitude cancellation due to the methods used to estimate the amplitude and record the signals. A: average EMG and root mean square (RMS) values were calculated for simulated bipolar and belly-tendon recordings. RMS amplitude and belly-tendon recordings both had greater amplitudes than the default condition (calculated as the average EMG from the bipolar recordings). B: percent amplitude cancellation was less when computed using RMS values and was greater for the simulated belly-tendon recordings compared with the bipolar recordings. C: when absolute EMG values were normalized to the value at maximal excitation, largest difference between the normalized values was <5% at intermediate levels of excitation.

 
Activation of the Motor Neuron Pool

The influence of discharge rate variability on amplitude cancellation was examined by varying the coefficient of variation for discharge rate from 10 to 40%. The percent amplitude cancellation, averaged across excitation levels, was 49.8, 49.7, 49.7, and 49.5% when the coefficient of variation for discharge rate was 10, 20, 30, and 40%, respectively. Due to the minor change in percent cancellation with discharge rate variability, only one activation pattern (coefficient of variation = 20%) was used across those conditions that did not involve a direct examination of activation pattern.

The influence of recruitment range was investigated by simulating narrow and broad upper limits of motor unit recruitment. A broad recruitment range decreased normalized EMG values (Fig. 5A) and increased amplitude cancellation (Fig. 5C) at intermediate levels of excitation, with the largest difference (16%) in normalized values occurring near the point where motor unit recruitment was complete for the narrow recruitment range (41% maximal excitation). The influence of peak discharge rate was examined by varying the peak rate based on the recruitment threshold of the motor unit (Fig. 5, B and D). In contrast to the decrease in normalized values at intermediate levels of excitation with a broad recruitment range, normalized EMG values were 17% greater than the default condition when high-threshold motor units achieved higher discharge rates, with the largest difference occurring at 30% maximal excitation (Fig. 5B).



View larger version (22K):
[in this window]
[in a new window]
 
Fig. 5. Variation in recruitment range and peak discharge rates changed the normalized EMG values at intermediate levels of excitation. A and C: normalized EMG values and percent amplitude cancellation varied with recruitment range. Default condition involved motor unit recruitment up to 41% of maximal excitation (narrow recruitment), whereas the upper limit of recruitment was 68% of maximal excitation for the broad range. Although cancellation for both conditions was greatest at maximal excitation, the largest difference in cancellation between the 2 conditions occurred near the point where motor unit recruitment was complete for the narrow recruitment range (41% maximal excitation). B and D: increased peak discharge rate for the high-threshold motor units resulted in an increase in normalized EMG values at intermediate excitation levels. Default condition involved peak discharge rates of 35–25 pulses per second (pps) for MU 1 and MU 120, respectively. Conditions were simulated that increased discharge rates (55–45 pps), increased preferentially the discharge rates of the high-threshold motor units (25–35 pps), and assigned the same discharge rate to all motor units (30 pps). Normalized EMG values were 17% greater than the default condition when high-threshold motor units achieved higher discharge rates than the low-threshold motor units, with the largest difference occurring at 30% maximal excitation.

 
The influence of motor unit synchronization on amplitude cancellation was examined by imposing a high level of synchrony across the motor unit pool (CIS = 2). The result was a slight reduction in percent amplitude cancellation due to a small increase in average EMG in the Cancellation condition (9% averaged across excitation levels), with no change in average EMG in the No-Cancellation condition (Table 2). Because these findings contrasted with those reported previously (61), the influence of variability in the timing of action potentials was examined. A representative motor unit potential was generated with 100 muscle fibers located 3 mm deep in the muscle, beneath a layer of skin (1 mm) and subcutaneous tissue (1.5 mm). Each fiber was assigned a conduction velocity of 4 m/s, and there was no variability in the location of the innervation zone or the end of the fibers. Three identical motor unit potentials were summed together, with no variability in the activation times (Fig. 6A). The amplitude of the sum was three times greater than the contributing potentials (perfect synchrony). However, when activation times of these same units were varied (MU 2 activated 1.7 ms after MU 1 and MU 3 activated 1.7 ms before MU 1), the average rectified amplitude decreased by 39.1% (Fig. 6B). When the three units were activated concurrently, but with variability in muscle fiber conduction velocities (conduction velocity for MU 1 = 4 m/s, MU 2 = 3.75 m/s, and MU 3 = 4.25 m/s), the arrival times of the potentials at the recording electrodes varied, and the average rectified amplitude was reduced by 49.1% (Fig. 6C). When the three units were activated concurrently, but with variability in the location of the innervation zone (innervation zone location for MU 1 = 40% along the length of the muscle, MU 2 = 5 mm distal to MU 1, and MU 3 = 5 mm proximal to MU 1), the average rectified amplitude was reduced by 23.5% (Fig. 6D). Because synchronized motor units have variable activation times, conduction velocities, and locations of the innervation zone, motor unit synchronization did not result in perfect summation of motor unit potentials and large increases in EMG amplitude.



View larger version (25K):
[in this window]
[in a new window]
 
Fig. 6. Variability in timing reduces the summation of action potentials. A, perfect synchrony: a representative motor unit was modeled by summing the action potentials from 100 muscle fibers positioned 3 mm below the surface of the muscle, with no variability in innervation zone, and all fibers assigned a conduction velocity of 4 m/s. In column 1, 3 identical motor unit potentials are summed together and the result is shown at bottom. B, imperfect synchrony: same motor units as in column 1, but with MU 2 activated 1.7 ms after MU 1 and MU 3 activated 1.7 ms before MU 1, which caused the times of arrival of the potentials at the recording electrodes to differ. Average rectified value of the resulting sum of the 3 motor units decreased by 39.1%. C, conduction velocity: each motor unit was activated at the same time, but MU 2 was assigned a conduction velocity of 3.75 m/s, and MU 3 was assigned a conduction velocity of 4.25 m/s. Average rectified value decreased by 49.1%. D, innervation zone: each motor unit was activated at the same time, but the innervation zone for MU 2 and MU 3 were located 5 mm distal and proximal to MU 1, respectively. Average rectified value decreased by 23.5%.

 
Influence of Selected Parameters on Amplitude Cancellation

Compared with the default condition, longer muscle fibers and thicker subcutaneous tissue decreased average EMG at maximal excitation by 10 and 42%, respectively (Fig. 7A). Amplitude cancellation was 10–15% greater for the longer fibers compared with shorter fibers, and thicker subcutaneous tissue caused only modest increases (≤5%) in percent amplitude cancellation (Fig. 7B). The substantial increase in amplitude cancellation with increased fiber length appears to be attributable to corresponding increases in motor unit potential duration (results presented below). Normalized values showed the same relation depicted in Fig. 4C, with the greatest difference being 3% at 30% excitation.



View larger version (20K):
[in this window]
[in a new window]
 
Fig. 7. Influence of muscle fiber length and subcutaneous tissue thickness on EMG amplitude and percent amplitude cancellation. A: increased fiber length (100 vs. 20 mm) and subcutaneous tissue thickness (3.5 vs. 1 mm) decreased average EMG. B: longer fibers and thicker subcutaneous tissue increased amplitude cancellation.

 
Those parameters associated with advancing age that had the largest influence on percent amplitude cancellation included the number of motor units and mean conduction velocity. Simulated old motor units (Table 2) and simulated decreases in the number of motor units from 120 to 60 decreased average EMG at maximal excitation by 34 and 13%, respectively. In contrast, a reduction in mean conduction velocity, which would accompany a decrease in fiber diameter, resulted in a 19% increase in EMG amplitude at maximal excitation (Fig. 8A). A reduction in the number of motor units decreased the amount of amplitude cancellation at maximal excitation from 61.7% in the default condition to 51.7% (Fig. 8B, Table 2). Because the simulated decrease in motor unit number involved reducing the number of muscle fibers from 66,000 to 46,200 fibers, simulations were performed in which motor unit number and muscle fiber number were varied independently. The results indicated that the reduction in cancellation was largely due to the decrease in motor unit number. In contrast, a reduction in mean conduction velocity resulted in greater amplitude cancellation (Fig. 8B, Table 2). Also, increased variability in muscle fiber conduction velocities within a motor unit resulted in greater amplitude cancellation (Table 2). The combined influence of age-associated changes in motor unit properties resulted in a net decrease in amplitude cancellation at low excitation levels (Fig. 8B, Table 2). Normalized values showed the same relation depicted in Fig. 4C, with the greatest difference being 4% at 50% excitation.



View larger version (18K):
[in this window]
[in a new window]
 
Fig. 8. Influence of age-associated changes in motor unit properties on EMG amplitude and percent amplitude cancellation (see Table 2 for data). A: loss of motor units (120 to 60) decreased average EMG, whereas a reduction in mean conduction velocity (4 to 3 m/s) increased average EMG. Simulated old motor units (Table 1) decreased average EMG. B: percent amplitude cancellation declined with a loss of motor units and increased with a reduction in mean conduction velocity. Simulated old motor units decreased percent amplitude cancellation at low excitation levels.

 
Figure 9 depicts the relations between different parameters of the motor unit potentials and the amount of cancellation that occurred at 50% maximal excitation in the 13 conditions when the activation pattern remained constant. Each data point denotes the mean value (amplitude, duration, or area) for all active motor units. The greatest correlations occurred for the duration measurements: potential duration explained 70% and normalized duration accounted for 79% of the variability in percent amplitude cancellation. Increased duration of the motor unit potentials resulted in greater overlap between the action potentials of different motor units and, therefore, greater cancellation. The influence of potential amplitude was weaker (r2 = 0.52) and caused small changes in average EMG in the Cancellation condition relative to the No-Cancellation condition. Because the calculation of percent amplitude cancellation was normalized to the values of average EMG in the No-Cancellation condition, increases in average EMG for the No-Cancellation condition, but not the Cancellation condition, decreased the percent amplitude cancellation.



View larger version (18K):
[in this window]
[in a new window]
 
Fig. 9. Association between motor unit action potential parameters and percent amplitude cancellation. Data correspond to the mean values for all active motor units in each of the 13 conditions in which the activation pattern remained constant. A: normalized motor unit action potential duration (area/amplitude). B: motor unit action potential duration. C: motor unit action potential area. D: motor unit action potential amplitude.

 

    DISCUSSION
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
This study identified the parameters of the surface EMG model that influenced amplitude cancellation and quantified the contribution of these parameters to the loss of signal information. Similar to a previous report (5), the results indicated that the surface EMG underestimated motor unit activity due to the loss of up to 62% of the amplitude of the signal when overlapping positive and negative phases of motor unit potentials cancelled one another. Also, changes in specific parameters of the model resulted in variation in the amount of amplitude cancellation, further confounding the interpretation of the surface EMG. The main characteristic of the motor unit potentials influencing amplitude cancellation was the duration of the potentials (Fig. 9). In addition, the relation between the average EMG and excitation was nonlinear such that the amplitude of the EMG increased at a gradually diminishing rate with an increase in excitation. This effect caused the normalized EMG values to overestimate motor unit activity by up to 13% at intermediate levels of excitation. However, when normalized values were compared across conditions, the difference was generally <5%, only increasing up to 17% at intermediate levels of excitation with variation in recruitment range and peak discharge rate. These findings have significant implications for the interpretation of experimental EMG data.

Identifying Neural Strategies from the Surface EMG

As a population index of motor unit activity, the surface EMG is used to examine neural strategies in a variety of conditions. For example, the surface EMG is used to compare the capabilities of young and old individuals and to evaluate adaptations in the same muscle before and after various interventions. In addition to the many factors known to influence EMG amplitude (16), however, the results of this study indicate that amplitude cancellation imposes a significant impediment to the identification of neural strategies from absolute surface EMG values.

The absolute amplitude of the surface EMG is commonly used in studies that compare young and old adults to determine whether differences exist in the ability of the nervous system to activate muscle (33, 45). For example, the finding of similar increases in maximal EMG amplitude in young and old adults after a training program has been interpreted as an increase in the neural drive to the trained muscles in both young and old adults (26). Such an interpretation is problematic because cancellation causes EMG amplitude to change due to factors that are independent of variation in neural drive, and EMG amplitude may be relatively insensitive to increases in motor unit activity at high excitation levels. For example, one would expect that an increase in the number of active motor units and muscle fibers would be accompanied by a large increase in the amplitude of the surface EMG. However, when the number of motor units was increased from 60 to 120 and the corresponding number of muscle fibers increased from 46,200 to 66,000, there was only a 15% increase in the amplitude of the EMG at maximal levels of excitation (Fig. 8A). The small increase in EMG amplitude with a doubling in motor unit number and a 43% increase in fiber number was due to an increase in amplitude cancellation (Fig. 8B). The results of this study demonstrate that amplitude cancellation confounds the interpretation of changes in absolute EMG amplitudes and limits the usefulness of measures that rely on absolute EMG values such as neuromuscular efficiency (9, 45).

Many factors influence EMG amplitude during muscle fatigue and confound its use as a measure of neural drive to muscle (10). It has been reported previously that although EMG amplitude increases during submaximal fatiguing contractions, the amplitude of the EMG is significantly less than maximum at the endurance limit for most subjects (23). It was proposed that the reduced EMG amplitude was due to a deficit in the ability to activate the muscle maximally, but two additional factors that enhance amplitude cancellation during a fatiguing contraction could also contribute to this observation. First, there is a decrease in muscle fiber conduction velocity that results in an increase in the duration of motor unit potentials during sustained contractions (54), leading to greater overlap between potentials (20) and an increase in amplitude cancellation (Figs. 8B and 9). However, the increase in potential duration can also augment EMG amplitude as the area of each motor unit potential increases (Fig. 8A). Second, sustained activation preferentially decreases the amplitude and area of the largest motor unit potentials (13), resulting in a more narrow range of amplitudes for the motor unit potentials. The current study demonstrated that a reduction in the range of potential amplitudes causes large increases in amplitude cancellation (Table 2). Hence, changes in the duration and range of amplitudes for the motor unit potentials can contribute to the inability of EMG amplitude to reach maximal levels after a fatiguing contraction, independent of any deficit in the capacity to activate the muscle.

Normalized EMG Values

It is common practice to normalize EMG values to the amplitude of the EMG at maximal levels of excitation so that comparisons can be made across muscles, between subjects, and between days. Normalized EMG values are less variable than absolute values, and therefore considered a more reliable index of muscle activation (60). Although the nonlinear increase in EMG values with increasing excitation has been demonstrated previously (5, 44, 49), the impact of this nonlinear relation on the normalized surface EMG values has not been addressed.

The greater probability of cancellation with increased motor unit activity has been referred to as a "saturation" effect (49) and as a downward nonlinearity (5) due to the failure of the EMG to increase at the same rate as motor unit activity. The current study demonstrates that amplitude cancellation did not increase linearly across excitation levels (Fig. 3), which resulted in an overestimation of motor neuron activity at intermediate levels of excitation when the EMG was normalized to maximal levels. Furthermore, variation in the range of amplitudes for the motor unit potentials did not alter the normalized EMG values and neither the method used to estimate the amplitude (RMS or average EMG) nor the recording technique (bipolar or belly-tendon) influenced the relation between normalized EMG and excitation (Fig. 4B). Therefore, the nonlinear relation between EMG and excitation was due to an increase in the probability of cancellation with increasing levels of motor unit activity.

The importance of this nonlinearity between normalized EMG and excitation likely depends on the type of comparison being made. Examination of those figures depicting percent amplitude cancellation across excitation levels (Figs. 4B; 5, B and D; 7B; and 8B) demonstrates that whereas the absolute amount of cancellation may vary across conditions, variation in specific conditions resulted in little change in normalized EMG values. This finding likely explains the reduced variability and increased reliability in normalized EMG values compared with absolute EMG values (41, 60). The two conditions that demonstrated the largest change in normalized EMG values (variation in recruitment range and peak discharge rates) resulted in a difference with the default condition of up to 17%, and this difference was limited to intermediate levels of excitation. The results of the study suggest caution in interpreting different levels of normalized surface EMG in muscles where recruitment range or peak discharge rates may vary, particularly at intermediate levels of excitation.

Motor Unit Synchronization

The current study found that EMG amplitude increased with motor unit synchronization, but the increase was much less than that reported previously (37, 61). There are at least two factors that may explain this difference. First, variability in the timing of synchronization between surface-recorded motor unit potentials is influenced by variability in motor unit conduction velocities, innervation zone locations, and tendon endings (Fig. 6). Each source of variability contributes to the failure of EMG amplitude to increase with motor unit synchronization. Second, constraining the imposed synchrony to motor units with similar recruitment thresholds reduced the chance of all motor units in the pool being synchronized with one another. Because high levels of synchrony involve many discharges that are nearly coincidental, the increase in overlap between potentials generates both in-phase and out-of-phase alignment between potentials and causes both increased and decreased amounts of amplitude cancellation with increases in motor unit synchronization. Even if higher levels of motor unit synchronization were imposed, a minimal increase in EMG amplitude would be expected due to the generation of both in-phase and out-of-phase alignments. In contrast, Zhou and Rymer (62) observed that modeling only variation in the timing of the synchronized motor units resulted in less out-of-phase alignments as the duration of the motor unit potentials increased and, therefore, there was less cancellation. Variation in conduction velocities and the random location of the innervation zones simulated in the current study influenced the increase in amplitude with synchronization. Thus systematic evaluation of the variability in the timing of the motor unit potentials suggests that motor unit synchronization has only a modest effect on EMG amplitude.

Limitations

There is a lack of information regarding the range of motor unit properties in a muscle (12), with limited information on the distribution of motor unit numbers, conduction velocities, peak discharge rates, and motor unit synchronization across the entire population. The information that does exist suggests great variability in many of the relevant physiological parameters that influence the EMG signal. The current study was limited in scope and identified selected motor unit properties that are likely to vary in different muscles and across different populations. Accordingly, measurements were made to provide realistic values for subcutaneous tissue thickness and fiber length. The surface EMG model of Fuglevand et al. (22) was updated with a model that has been validated to produce realistic motor unit potentials (15). Also the modeling approach allowed the assessment of the contribution of individual parameters to the surface EMG, something rarely possible during experimental studies. Despite the assumptions used in the current modeling study, however, the results are consistent with the one experimental study that systematically quantified amplitude cancellation using an experimental protocol (5).

In summary, up to 62% of the surface EMG signal amplitude was lost due to cancellation. The degree of signal loss depended on selected physiological parameters, especially the number of active motor units and the duration of the action potentials. Nonetheless, normalization of the surface EMG amplitude to the values obtained with maximal activation increases the reliability of the measurement. Although normalized EMG values overestimate the amount of motor unit activity at intermediate levels of activation by up to 13%, the importance of this nonlinearity is likely minimized as normalized values are largely invariant to changes in different parameters. These results have implications for identifying neural strategies from the surface EMG, particularly when comparing absolute EMG values or when assessing fatigue-related changes in EMG.


    GRANTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
This work was supported by National Institute of Neurological Disorders and Stroke Award NS-42734 and National Institute on Aging (NIA) Award AG-09000 to R. M. Enoka and European Space Agency Project C15097 [GenBank] /01/NL/SN to D. Farina and R. Merletti. K. G. Keenan was supported by NIA Award T32-AG-00279.


    FOOTNOTES
 

Address for reprint requests and other correspondence: R. M. Enoka, Dept. of Integrative Physiology, Univ. of Colorado, Boulder, CO 80309-0354 (E-mail: enoka{at}colorado.edu)

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.


    REFERENCES
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 

  1. Adrian ED. Interpretation of the electromyogram. Lancet, \?\June 13: \?\1229–1233; June 20: 1282–1286, 1925.
  2. Bellemare F, Woods JJ, Johansson R, and Bigland-Ritchie B. Motor-unit discharge rates in maximal voluntary contractions of three human muscles. J Neurophysiol 50: 1380–1392, 1983.[Abstract/Free Full Text]
  3. Bigland-Ritchie B. EMG/force relations and fatigue of human voluntary contractions. Exerc Sport Sci Rev 9: 75–117, 1981.[Medline]
  4. Blok JH, Stegeman DF, and van Oosterom A. Three-layer volume conductor model and software package for applications in surface electromyography. Ann Biomed Engin 30: 566–577, 2002.[CrossRef][Web of Science][Medline]
  5. Day SJ and Hulliger M. Experimental simulation of cat electromyogram: evidence for algebraic summation of motor-unit action-potential trains. J Neurophysiol 86: 2144–2158, 2001.[Abstract/Free Full Text]
  6. De Luca CJ, LeFever RS, McCue MP, and Xenakis AP. Behaviour of human motor units in different muscles during linearly varying contractions. J Physiol 329: 113–128, 1982.[Abstract/Free Full Text]
  7. De Luca CJ, LeFever RS, McCue MP, and Xenakis AP. Control scheme governing concurrently active human motor units during voluntary contractions. J Physiol 329: 129–142, 1982.[Abstract/Free Full Text]
  8. Dennett X and Fry HJ. Overuse syndrome: a muscle biopsy study. Lancet 331: 905–908, 1988.[CrossRef]
  9. Deschenes MR, Giles JA, McCoy RW, Volek JS, Gomez AL, and Kraemer WJ. Neural factors account for strength decrements observed after short-term muscle unloading. Am J Physiol Regul Integr Comp Physiol 282: R578–R583, 2002.[Abstract/Free Full Text]
  10. Dimitrova NA and Dimitrov GV. Interpretation of EMG changes with fatigue: facts, pitfalls, and fallacies. J Electromyogr Kinesiol 13: 13–36, 2003.[CrossRef][Web of Science][Medline]
  11. Doherty TJ, Vandervoort AA, Taylor AW, and Brown WF. Effects of motor unit losses on strength in older men and women. J Appl Physiol 74: 868–874, 1993.[Abstract/Free Full Text]
  12. Enoka RM and Fuglevand AJ. Motor unit physiology: some unresolved issues. Muscle Nerve 24: 4–17, 2001.[CrossRef][Web of Science][Medline]
  13. Enoka RM, Trayanova N, Laouris Y, Bevan L, Reinking RM, and Stuart DG. Fatigue-related changes in motor unit action potentials of adult cats. Muscle Nerve 15: 138–150, 1992.[CrossRef][Web of Science][Medline]
  14. Farina D, Cescon C, and Merletti R. Influence of anatomical, physical, and detection-system parameters on surface EMG. Biol Cybern 86: 445–456, 2002.[CrossRef][Web of Science][Medline]
  15. Farina D and Merletti R. A novel approach for precise simulation of the EMG signal detected by surface electrodes. IEEE Trans Biomed Eng 48: 637–646, 2001.[CrossRef][Web of Science][Medline]
  16. Farina D, Merletti R, and Enoka RM. The extraction of neural strategies from the surface EMG. J Appl Physiol 96: 1486–1495, 2004.[Abstract/Free Full Text]
  17. Farina D, Merletti R, Indino B, Nazzaro M, and Pozzo M. Surface EMG crosstalk between knee extensor muscles: experimental and model results. Muscle Nerve 26: 681–695, 2002.[CrossRef][Web of Science][Medline]
  18. Farina D and Rainoldi A. Compensation of the effect of sub-cutaneous tissue layers on surface EMG: a simulation study. Med Eng Phys 21: 487–497, 1999.[CrossRef][Web of Science][Medline]
  19. Fortier PA. Use of spike triggered averaging of muscle activity to quantify inputs to motoneuron pools. J Neurophysiol 72: 248–265, 1994.[Abstract/Free Full Text]
  20. Fuglevand AJ. The role of the sarcolemma action potential in fatigue. In: Fatigue: Neural and Muscular Mechanisms, edited by Gandevia SC, Enoka RM, McComas AJ, Stuart DG, and Thomas CK. New York: Plenum, 1995, p. 101–108.
  21. Fuglevand AJ, Winter DA, and Patla AE. Models of recruitment and rate coding organization in motor-unit pools. J Neurophysiol 70: 2470–2488, 1993.[Abstract/Free Full Text]
  22. Fuglevand AJ, Winter DA, Patla AE, and Stashuk D. Detection of motor unit action potentials with surface electrodes: influence of electrode size and spacing. Biol Cybern 67: 143–153, 1992.[CrossRef][Web of Science][Medline]
  23. Fuglevand AJ, Zackowski KM, Huey KA, and Enoka RM. Impairment of neuromuscular propagation during human fatiguing contractions at submaximal forces. J Physiol 460: 549–572, 1993.[Abstract/Free Full Text]
  24. Gazzoni M, Farina D, and Merletti R. A new method for the extraction and classification of single motor unit action potentials from surface EMG signals. J Neurosci Methods 136: 165–177, 2004.[CrossRef][Web of Science][Medline]
  25. Gydikov A and Kosarov D. Some features of different motor units in human biceps brachii. Pflügers Arch 347: 75–88, 1974.[CrossRef][Web of Science][Medline]
  26. Hakkinen K, Newton RU, Gordon SE, McCormick M, Volek JS, Nindl BC, Gotshalk LA, Campbell WW, Evans WJ, Hakkinen A, Humphries BJ, and Kraemer WJ. Changes in muscle morphology, electromyographic activity, and force production characteristics during progressive strength training in young and older men. J Gerontol B Psychol Sci Soc Sci 53: B415–423, 1998.
  27. Hansen NL, Hansen S, Christensen LO, Petersen NT, and Nielsen JB. Synchronization of lower limb motor unit activity during walking in human subjects. J Neurophysiol 86: 1266–1276, 2001.[Abstract/Free Full Text]
  28. Heckman CJ and Binder MD. Computer simulation of the steady-state input-output function of the cat medial gastrocnemius motoneuron pool. J Neurophysiol 65: 952–967, 1991.[Abstract/Free Full Text]
  29. Henneman E. Relation between size of neurons and their susceptibility to discharge. Science 126: 1345–1347, 1957.[Free Full Text]
  30. Kamen G, Sison SV, Du CC, and Patten C. Motor unit discharge behavior in older adults during maximal-effort contractions. J Appl Physiol 79: 1908–1913, 1995.[Abstract/Free Full Text]
  31. Kanda K and Hashizume K. Changes in properties of the medial gastrocnemius motor units in aging rats. J Neurophysiol 61: 737–746, 1989.[Abstract/Free Full Text]
  32. Kanda K and Hashizume K. Factors causing difference in force output among motor units in the rat medial gastrocnemius muscle. J Physiol 448: 677–695, 1992.[Abstract/Free Full Text]
  33. Keen DA, Yue GH, and Enoka RM. Training-related enhancement in the control of motor output in elderly humans. J Appl Physiol 77: 2648–2658, 1994.[Abstract/Free Full Text]
  34. Keenan KG, Farina D, Maluf KS, Merletti R, and Enoka RM. Age-associated changes in motor unit properties reduce signal cancellation in the simulated electromyogram. Soc Neurosci Abstr 914.11, 2003.
  35. Kilner JM, Baker SN, Salenius S, Jousmaki V, Hari R, and Lemon RN. Task-dependent modulation of 15–30 Hz coherence between rectified EMGs from human hand and forearm muscles. J Physiol 516: 559–570, 1999.[Abstract/Free Full Text]
  36. Klein CS, Marsh GD, Petrella RJ, and Rice CL. Muscle fiber number in the biceps brachii muscle of young and old men. Muscle Nerve 28: 62–68, 2003.[CrossRef][Web of Science][Medline]
  37. Kleine BU, Stegeman DF, Mund D, and Anders C. Influence of motoneuron firing synchronization on SEMG characteristics in dependence of electrode position. J Appl Physiol 91: 1588–1599, 2001.[Abstract/Free Full Text]
  38. Kukulka CG and Clamann HP. Comparison of the recruitment and discharge properties of motor units in human brachial biceps and adductor pollicis during isometric contractions. Brain Res 219: 45–55, 1981.[CrossRef][Web of Science][Medline]
  39. Lexell J, Henriksson-Larsen K, Winblad B, and Sjostrom M. Distribution of different fiber types in human skeletal muscles: effects of aging studied in whole muscle cross sections. Muscle Nerve 6: 588–595, 1983.[CrossRef][Web of Science][Medline]
  40. Masuda T and Sadoyama T. The propagation of single motor unit action potentials detected by a surface electrode array. Electroencephalogr Clin Neurophysiol 63: 590–598, 1986.[CrossRef][Web of Science][Medline]
  41. Mathiassen SE, Winkel J, and Hagg GM. Normalization of surface EMG amplitude from the upper trapezius in ergonomic studies—a review. J Electromyogr Kinesiol 5: 197–226, 1995.
  42. Milner-Brown HS, Stein RB, and Yemm R. The orderly recruitment of human motor units during voluntary isometric contractions. J Physiol 230: 359–370, 1973.[Abstract/Free Full Text]
  43. Monster AW and Chan H. Isometric force production by motor units of extensor digitorum communis muscle in man. J Neurophysiol 40: 1432–1443, 1977.[Free Full Text]
  44. Moore AD. Synthesized EMG waves and their implications. Am J Phys Med 46: 1302–1316, 1967.[Medline]
  45. Moritani T and deVries HA. Potential for gross muscle hypertrophy in older men. J Gerontol A Biol Sci Med Sci 35: 672–682, 1980.
  46. Nandedkar SD and Sanders DB. Simulation of myopathic motor unit action potentials. Muscle Nerve 12: 197–202, 1989.[CrossRef][Web of Science][Medline]
  47. Nandedkar SD, Sanders DB, and Stalberg EV. Automatic analysis of the electromyographic interference pattern. Part I: development of quantitative features. Muscle Nerve 9: 431–439, 1986.[CrossRef][Web of Science][Medline]
  48. Nordstrom MA, Fuglevand AJ, and Enoka RM. Estimating the strength of common input to human motoneurons from the cross-correlogram. J Physiol 453: 547–574, 1992.[Abstract/Free Full Text]
  49. Person RS and Libkind MS. Simulation of electromyograms showing interference patterns. Electroencephalogr Clin Neurophysiol 28: 625–632, 1970.[CrossRef][Web of Science][Medline]
  50. Roeleveld K, Blok JH, Stegeman DF, and van Oosterom A. Volume conduction models for surface EMG; confrontation with measurements. J Electromyogr Kinesiol 7: 221–232, 1997.[CrossRef][Web of Science][Medline]
  51. Rosenfalck P. Intra- and extracellular potential fields of active nerve and muscle fibers. Acta Physiol Scand Suppl 47: 239–246, 1969.
  52. Rosler KM, Petrow E, Mathis J, Aranyi Z, Hess CW, and Magistris MR. Effect of discharge desynchronization on the size of motor evoked potentials: an analysis. Clin Neurophysiol 113: 1680–1687, 2002.[CrossRef][Web of Science][Medline]
  53. Shinohara M, Keenan KG, and Enoka RM. Contralateral activity in a homologous hand muscle during voluntary contractions is greater in old adults. J Appl Physiol 94: 966–974, 2003.[Abstract/Free Full Text]
  54. Stalberg E. Propagation velocity in human muscle fibers in situ. Acta Physiol Scand Suppl 287: 1–112, 1966.[Medline]
  55. Stegeman DF, Blok JH, Hermens HJ, and Roeleveld K. Surface EMG models: properties and applications. J Electromyogr Kinesiol 10: 313–326, 2000.[CrossRef][Web of Science][Medline]
  56. Taylor AM, Christou EA, and Enoka RM. Multiple features of motor-unit activity influence force fluctuations during isometric contractions. J Neurophysiol 90: 1350–1361, 2003.[Abstract/Free Full Text]
  57. Taylor AM, Steege JW, and Enoka RM. Motor-unit synchronization alters spike-triggered average force in simulated contractions. J Neurophysiol 88: 265–276, 2002.[Abstract/Free Full Text]
  58. Tomlinson BE and Irving D. The numbers of limb motor neurons in the human lumbosacral cord throughout life. J Neurol Sci 34: 213–219, 1977.[CrossRef][Web of Science][Medline]
  59. Van Cutsem M, Feiereisen P, Duchateau J, and Hainaut K. Mechanical properties and behaviour of motor units in the tibialis anterior during voluntary contractions. Can J Appl Physiol 22: 585–597, 1997.[Web of Science][Medline]
  60. Yang JF and Winter DA. Electromyography reliability in maximal and submaximal isometric contractions. Arch Phys Med Rehabil 64: 417–420, 1983.[Web of Science][Medline]
  61. Yao W, Fuglevand AJ, and Enoka RM. Motor-unit synchronization increases EMG amplitude and decreases force steadiness of simulated contractions. J Neurophysiol 83: 441–452, 2000.[Abstract/Free Full Text]
  62. Zhou P and Rymer WZ. Factors governing the form of the relation between muscle force and the electromyogram (EMG): a simulation study. J Neurophysiol 92: 2878–2886, 2004.[Abstract/Free Full Text]
  63. Zijdewind I, de Groot MC, and Kernell D. Task-related variations in motoneuronal drive to a human intrinsic hand muscle. Neurosci Lett 242: 139–142, 1998.[CrossRef][Web of Science][Medline]



This article has been cited by other articles:


Home page
J. Neurophysiol.Home page
K. Schmitt, C. DelloRusso, and R. F. Fregosi
Force-EMG Changes During Sustained Contractions of a Human Upper Airway Muscle
J Neurophysiol, February 1, 2009; 101(2): 558 - 568.
[Abstract] [Full Text] [PDF]


Home page
J. Physiol.Home page
M. Amann, L. T. Proctor, J. J. Sebranek, D. F. Pegelow, and J. A. Dempsey
Opioid-mediated muscle afferents inhibit central motor drive and limit peripheral muscle fatigue development in humans
J. Physiol., January 1, 2009; 587(1): 271 - 283.
[Abstract] [Full Text] [PDF]


Home page
J. Physiol.Home page
J. Duchateau and R. M. Enoka
Neural control of shortening and lengthening contractions: influence of task constraints
J. Physiol., December 15, 2008; 586(24): 5853 - 5864.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
M. Amann, L. T. Proctor, J. J. Sebranek, M. W. Eldridge, D. F. Pegelow, and J. A. Dempsey
Somatosensory feedback from the limbs exerts inhibitory influences on central neural drive during whole body endurance exercise
J Appl Physiol, December 1, 2008; 105(6): 1714 - 1724.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
Z. A. Riley, S. Baudry, and R. M. Enoka
Reflex Inhibition in Human Biceps Brachii Decreases With Practice of a Fatiguing Contraction
J Neurophysiol, November 1, 2008; 100(5): 2843 - 2851.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
J. M. Dundon, J. Cirillo, and J. G. Semmler
Low-frequency fatigue and neuromuscular performance after exercise-induced damage to elbow flexor muscles
J Appl Physiol, October 1, 2008; 105(4): 1146 - 1155.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
D. Farina, F. Negro, M. Gazzoni, and R. M. Enoka
Detecting the Unique Representation of Motor-Unit Action Potentials in the Surface Electromyogram
J Neurophysiol, September 1, 2008; 100(3): 1223 - 1233.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
S. K. Hunter, T. Yoon, J. Farinella, E. E. Griffith, and A. V. Ng
Time to task failure and muscle activation vary with load type for a submaximal fatiguing contraction with the lower leg
J Appl Physiol, August 1, 2008; 105(2): 463 - 472.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
M. Post, S. Bayrak, D. Kernell, and I. Zijdewind
Contralateral muscle activity and fatigue in the human first dorsal interosseous muscle
J Appl Physiol, July 1, 2008; 105(1): 70 - 82.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
D. Farina, C. Cescon, F. Negro, and R. M. Enoka
Amplitude Cancellation of Motor-Unit Action Potentials in the Surface Electromyogram Can Be Estimated With Spike-Triggered Averaging
J Neurophysiol, July 1, 2008; 100(1): 431 - 440.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
T. Rudroff, D. Staudenmann, and R. M. Enoka
Electromyographic measures of muscle activation and changes in muscle architecture of human elbow flexors during fatiguing contractions
J Appl Physiol, June 1, 2008; 104(6): 1720 - 1726.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
M. Klass, M. Levenez, R. M. Enoka, and J. Duchateau
Spinal Mechanisms Contribute to Differences in the Time to Failure of Submaximal Fatiguing Contractions Performed With Different Loads
J Neurophysiol, March 1, 2008; 99(3): 1096 - 1104.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
M. Klass, S. Baudry, and J. Duchateau
Age-related decline in rate of torque development is accompanied by lower maximal motor unit discharge frequency during fast contractions
J Appl Physiol, March 1, 2008; 104(3): 739 - 746.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
M. Levenez, S. J. Garland, M. Klass, and J. Duchateau
Cortical and Spinal Modulation of Antagonist Coactivation During a Submaximal Fatiguing Contraction in Humans
J Neurophysiol, February 1, 2008; 99(2): 554 - 563.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
T. J. Dartnall, M. A. Nordstrom, and J. G. Semmler
Motor Unit Synchronization Is Increased in Biceps Brachii After Exercise-Induced Damage to Elbow Flexor Muscles
J Neurophysiol, February 1, 2008; 99(2): 1008 - 1019.
[Abstract] [Full Text] [PDF]


Home page
J. Physiol.Home page
M. Amann and J. A. Dempsey
Locomotor muscle fatigue modifies central motor drive in healthy humans and imposes a limitation to exercise performance
J. Physiol., January 1, 2008; 586(1): 161 - 173.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
K. G. Keenan and F. J. Valero-Cuevas
Experimentally Valid Predictions of Muscle Force and EMG in Models of Motor-Unit Function Are Most Sensitive to Neural Properties
J Neurophysiol, September 1, 2007; 98(3): 1581 - 1590.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
C. Suetta, P. Aagaard, S. P. Magnusson, L. L. Andersen, S. Sipila, A. Rosted, A. K. Jakobsen, B. Duus, and M. Kjaer
Muscle size, neuromuscular activation, and rapid force characteristics in elderly men and women: effects of unilateral long-term disuse due to hip-osteoarthritis
J Appl Physiol, March 1, 2007; 102(3): 942 - 948.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Regul. Integr. Comp. Physiol.Home page
K. Katayama, M. Amann, D. F. Pegelow, A. J. Jacques, and J. A. Dempsey
Effect of arterial oxygenation on quadriceps fatigability during isolated muscle exercise
Am J Physiol Regulatory Integrative Comp Physiol, March 1, 2007; 292(3): R1279 - R1286.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
J. Duchateau, J. G. Semmler, and R. M. Enoka
Training adaptations in the behavior of human motor units
J Appl Physiol, December 1, 2006; 101(6): 1766 - 1775.
[Abstract] [Full Text] [PDF]


Home page
J. Physiol.Home page
B. Pasquet, A. Carpentier, and J. Duchateau
Specific modulation of motor unit discharge for a similar change in fascicle length during shortening and lengthening contractions in humans
J. Physiol., December 1, 2006; 577(2): 753 - 765.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
M. Amann, L. M. Romer, D. F. Pegelow, A. J. Jacques, C. J. Hess, and J. A. Dempsey
Effects of arterial oxygen content on peripheral locomotor muscle fatigue
J Appl Physiol, July 1, 2006; 101(1): 119 - 127.
[Abstract] [Full Text] [PDF]


Home page
Br. J. Sports. Med.Home page
J P Weir, T W Beck, J T Cramer, T J Housh, T D Noakes, A St Clair Gibson, and E V Lambert
Is fatigue all in your head? A critical review of the central governor model * Commentary.
Br. J. Sports Med., July 1, 2006; 40(7): 573 - 586.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
K. G. Keenan, D. Farina, R. Merletti, and R. M. Enoka
Amplitude cancellation reduces the size of motor unit potentials averaged from the surface EMG
J Appl Physiol, June 1, 2006; 100(6): 1928 - 1937.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
C. K. Thomas, R. S. Johansson, and B. Bigland-Ritchie
EMG Changes in Human Thenar Motor Units With Force Potentiation and Fatigue
J Neurophysiol, March 1, 2006; 95(3): 1518 - 1526.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
B. Pasquet, A. Carpentier, and J. Duchateau
Change in Muscle Fascicle Length Influences the Recruitment and Discharge Rate of Motor Units During Isometric Contractions
J Neurophysiol, November 1, 2005; 94(5): 3126 - 3133.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
M. Shinohara, C. T. Moritz, M. A. Pascoe, and R. M. Enoka
Prolonged muscle vibration increases stretch reflex amplitude, motor unit discharge rate, and force fluctuations in a hand muscle
J Appl Physiol, November 1, 2005; 99(5): 1835 - 1842.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
M. Levenez, C. Kotzamanidis, A. Carpentier, and J. Duchateau
Spinal reflexes and coactivation of ankle muscles during a submaximal fatiguing contraction
J Appl Physiol, September 1, 2005; 99(3): 1182 - 1188.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
K. S. Maluf and R. M. Enoka
Task failure during fatiguing contractions performed by humans
J Appl Physiol, August 1, 2005; 99(2): 389 - 396.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
M. Gazzoni, F. Camelia, and D. Farina
Conduction Velocity of Quiescent Muscle Fibers Decreases During Sustained Contraction
J Neurophysiol, July 1, 2005; 94(1): 387 - 394.
[Abstract] [Full Text] [PDF]


Home page
Journals of Gerontology Series A: Biological Sciences and Medical SciencesHome page
E. Simoneau, A. Martin, and J. Van Hoecke
Muscular Performances at the Ankle Joint in Young and Elderly Men
J. Gerontol. A Biol. Sci. Med. Sci., April 1, 2005; 60(4): 439 - 447.
[Abstract] [Full Text] [PDF]


Home page
J. Physiol.Home page
C Westad and R. H Westgaard
The influence of contraction amplitude and firing history on spike-triggered averaged trapezius motor unit potentials
J. Physiol., February 1, 2005; 562(3): 965 - 975.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
98/1/120    most recent
00894.2004v1
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (85)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Keenan, K. G.
Right arrow Articles by Enoka, R. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Keenan, K. G.
Right arrow Articles by Enoka, R. M.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2005 by the American Physiological Society.