The effects of motor unit synchronization on the power spectrum of the electromyogram

A realistic model for two synchronized motor unit action potential trains (MUAPT) is presented in which the variability of the time difference between corresponding action potentials (hereafter denoted by delay) is taken into account. Specifically, this delay is modeled as a continuous random variable that may assume both positive and negative values.Expressions are derived for the auto- and cross-power spectra of two such trains using their relations with the auto- and cross-correlation functions, respectively, with which they form Fourier transform pairs.The results show that the auto- and the cross-power spectra of two such synchronized MUAPTs differ from the auto- and the cross-spectra of two independent MUAPTs. The contribution of the statistics of the interpulse intervals to one of the autopower spectra is smaller and the cross-power spectra no longer reduce to a Dirac -function at the origin but are now determined by the other auto-power spectrum and by the Fourier transform of the density function associated with the time difference between corresponding action potentials. As a consequence of this change in the cross-power spectra synchronization leads to an absolute increase of power at low frequencies and to a relative decrease of power at high frequencies.The results are then generalized to electromyograms (EMG) composed of more than just two MUAPTs and illustrated with simulated power spectra with which the theory shows excellent agreement.     

Changes in the surface electromyogram during increasing isometric shoulder forward flexions

When using electromyographic techniques in the evaluation of muscular load it is necessary to determine the mathematical relationship between the torque and the amplitude of the electromyographic signal. Isometric gradually increasing contractions up to 100% MVC can then be used. Often more than linear increases for the amplitude (RMS)--force regression have been reported. The present study was designed to test whether changes in power spectral density function take place during a gradually increasing isometric contraction (duration 10 s). Twenty-two clinically healthy females performed an increasing isometric shoulder forward flexion for 10 s using an isokinetic dynamometer. Electromyographic activity was measured in trapezius, deltoid, infraspinatus and biceps brachii using surface electrodes. Mean torque values were determined together with mean power frequency (MPF) and root mean square values (RMS) from the EMG signals for each 256 ms period. The RMS-torque regressions showed higher regression coefficients during the 6th to 9th sec than during the first 5 s. No significant correlation existed between MPF for the four muscles and the torque. A gradual decrease in MPF was generally found from the 6th s. It is concluded that this decrease in power spectral density function might have contributed to the significantly higher regression coefficient for the RMS torque regression at the high output part of the gradually increasing isometric contraction.     

Effects of immobilization on electromyogram power spectrum changes during fatigue.

The maximal force and median frequency (MF) of the electromyogram (EMG) power density spectrum (PDS) have been compared in disused (6 weeks' immobilization) and control (contralateral) human adductor pollicis muscles during fatigue induced by voluntary or electrically-triggered (30 Hz) contractions. The results indicated that after 6 weeks' immobilization, MF was not significantly different in disused and control muscles although the force and integrated EMG were drastically reduced during a maximal voluntary contraction (MVC; by 55% and 45%, respectively, n = 8). During sustained 60 s MVC, the force decreased at the same rate in immobilized and control muscles, but the shift of MF towards lower frequency values was smaller (P less than 0.05) in disused muscle as compared to control by (14% vs 28%, respectively). In electrically-induced fatigue, the force decrease and the MF shift were larger after inactivity (41% and 43% in one subject, and 50% and 54% in the other subject, respectively) as compared to control (29% and 34% in one subject, and 37% and 38% in the other subject, respectively). These results emphasize the caution that should be exercised when EMG signals are quantified by computing the power density spectrum. The different effects of fatigue during voluntary and electrically-imposed contractions in disused and control muscles indicated that immobilization induced changes in the neural command for the contraction which compensated, at least in part, for its decreased contractile efficiency and resistance to fatigue.     

Motor unit recruitment strategy of knee antagonist muscles in a step-wise, increasing isometric contraction

The purpose of this study was to determine if differences exist between the control strategies of two antagonist thigh muscles during knee flexion and extension muscular coactivation. Surface myoelectric signal (MES) of the quadriceps (rectus femoris) and the hamstrings (semitendinosus) were obtained from both muscles while performing step-wise increasing contractions during flexion and extension with the knee at 1.57 rad of flexion (90 degrees). The median frequency of the power density spectrum, which is related to the average muscle fiber action potential conduction velocity and therefore to motor unit recruitment, was calculated from each MES. The results suggest that, in all the subjects tested, when the muscle acts as antagonist most motor units are recruited up to 50% of the maximal voluntary force, whereas when the muscle acts as antagonist motor units are recruited up to 40% of the maximal voluntary force. The force range past 40–50% of the maximal force is also characterized by differences between the agonist/antagonist.     

Effects of oxygen fraction in inspired air on force production and electromyogram activity during ergometer rowing

Six male rowers rowed maximally for 2500 m in ergometer tests during normoxia (fractional concentration of oxygen in inspired air, F _IO₂ 0.209), in hyperoxia (F _IO₂ 0.622) and in hypoxia (F _IO₂ 0.158) in a randomized single-blind fashion. Oxygen consumption (V˙O₂), force production of strokes as well as integrated electromyographs (iEMG) and mean power frequency (MPF) from seven muscles were measured in 500-m intervals. The iEMG signals from individual muscles were summed to represent overall electrical activity of these muscles (sum-iEMG). Maximal force of a stroke (F _max) decreased from the 100% pre-exercise maximal value to 67 (SD 12)%, 63 (SD 15)% and 76 (SD 13)% (P<0.05 to normoxia, ANOVA) and impulse to 78 (SD 4)%, 75 (SD 14)% and 84 (SD 7)% (P<0.05) in normoxia, hypoxia and hyperoxia, respectively. A strong correlation between F _max and V˙O₂ was found in normoxia but not in hypoxia and hyperoxia. The mean sum-iEMG tended to be lower (P<0.05) in hypoxia than in normoxia but hyperoxia had no significant effect on it. In general, F _IO₂ did not affect MPF of individual muscles. In conclusion, it was found that force output during ergometer rowing was impaired during hypoxia and improved during hyperoxia when compared with normoxia. Moreover, the changes in force output were only partly accompanied by changes in muscle electrical activity as sum-iEMG was affected by hypoxic but not by hyperoxic gas. The lack of a significant correlation between F _max and V˙O₂ during hypoxia and hyperoxia may suggest a partial uncoupling of these processes and the existence of other limiting factors in addition to V˙O₂. Accepted: 2 June 1997     

Intramuscular and surface electromyogram changes during muscle fatigue

Twelve male subjects were tested to determine the effects of motor unit (MU) recruitment and firing frequency on the surface electromyogram (EMG) frequency power spectra during sustained maximal voluntary contraction (MVC) and 50% MVC of the biceps brachii muscle. Both the intramuscular MU spikes and surface EMG were recorded simultaneously and analyzed by means of a computer-aided intramuscular spike amplitude-frequency histogram and frequency power spectral analysis, respectively. Results indicated that both mean power frequency (MPF) and amplitude (rmsEMG) of the surface EMG fell significantly (P less than 0.001) together with a progressive reduction in MU spike amplitude and firing frequency during sustained MVC. During 50% MVC there was a significant decline in MPF (P less than 0.001), but this decline was accompanied by a significant increase in rmsEMG (P less than 0.001) and a progressive MU recruitment as evidenced by an increased number of MUs with relatively large spike amplitude. Our data suggest that the surface EMG amplitude could better represent the underlying MU activity during muscle fatigue and the frequency powers spectral shift may or may not reflect changes in MU recruitment and rate-coding patterns.     

Reproducibility of the mean power frequency of the surface electromyogram

There is evidence in the literature that the decrease of mean power frequency (MPF) during exercise is greater as a muscle become more fatigued. After strenuous exercise this phenomenon can last several days. It is usually assumed, however, that the MPF has a good reproducibility. In this study the reproducibility of the MPF of the surface electromyogram of the biceps brachii muscle was investigated for five subjects on 5 successive days. Force level, muscle length and skin temperature were kept constant. The results show that interindividual differences in MPF were large (SD 11.5 Hz). However, during these 5 days, the range in MPF for individual subjects was small. The SD of the trials within subjects and days was 2.0 Hz, while the SD trials excluded). It is hypothesized that this SD may be due to variations in the electrode replacement. It is concluded that the variability in MPF for a subject is small compared to the decrease of the MPF associated with muscle fatigue and which can therefore be determined reliably during longitudinal studies.     

Comparison of muscle force,muscle endurance,and electromyogram activity during an expedition at high altitude

Handgrip force (HF), maximal pinch force (MF), muscle endurance (ME), and the median power frequency (MdPF) of the activity shown in the electromyogram (EMG) were studied at various altitudes in eight normal healthy subjects. MF and ME were measured between the index finger and thumb, and all measurements were obtained at altitudes ranging from 610 to 4860 m during an expedition in the Qinghai Plateau in China. With the change in altitude HF, ME, and MF showed no significant change. Compared to the MdPF at 2260 m on ascent, the MdPF at other altitudes showed a significant decrease (P<0.01). Thus, we conclude that muscle performance (HF, MF, and ME) was not affected by the environment at high altitude. However, MdPF was affected and the mean MdPF at 610 m after the expedition did not recover to initial values of MdPF. We suggest these results may have been affected by fatigue and chronic exposure to the hypobaric hypoxic environment, since the members of the expedition party expressed feelings of sluggishness and fatigue after the expedition.     

Motor unit activity during long-lasting intermittent muscle contractions in humans

Changes accompanying long-lasting intermittent muscle contractions (30%–50% of the maximal) were investigated by tracing the activity of 38 motor units (MU) of the human biceps brachii muscle recorded from fine-wire branched electrodes. The motor task was a continuous repetition of ramp-and-hold cycles of isometric flexion contractions. During ramp-up phases a significant decline in recruitment thresholds was found with no changes in the discharge pattern. During ramp-down phases the unchanged mean value of derecruitment thresholds during the task was accompanied by increased duration of the last two interspike intervals (ISI). These findings would suggest that during fatigue development the main compensatory mechanism during ramp-up contractions is space coding while for ramp-down contractions it is rate coding. During the steady-state phases the mean value of ISI, as well as the firing variability, had increased by the end of the task in most of the MU investigated . In addition 17 recruited MU were also investigated. These units revealed a lower initial discharge rate and a faster decrease in the mean discharge rate with the development of fatigue. The gradual reduction of the recruitment threshold of already active MU and the recruitment of new units demonstrated an increased excitability of the motorneuron pool during fatigue. A typical recruitment pattern (a first short ISI followed by a long one) was observed during ramp-up contractions in units active from the very beginning of the task, as well as during sustained contractions at the onset of the stable discharge of the additionally recruited MU. Accepted: 23 September 1997     

Characteristics of power spectrum density function of EMG during muscle contraction below 30%MVC

The aim of the study was to quantify changes in PSDF frequency bands of the EMG signal and EMG parameters such as MF, MPF and zero crossing, with an increase in the level of muscle contractions in the range from 0.5% to 30% RMS_max and to determine the frequency bands with the lowest dependency on RMS level so that this could be used in investigating muscle fatigue.Sixteen men, aged from 23 to 33 years old (mean 26.1), who participated in the study performed two force exertion tests. Fragments of EMG which corresponded to the levels of muscle contraction of 0.5%, 1%, 2.5%, 5%, 10%, 15%, 20%, 25%, 30% RMS_max registered from left and right trapezius pars descendents (TP) and left and right extensor digitorum superficialis (ED) muscles were selected for analysis. The analysis included changes in standard parameters of the EMG signal and changes in PSDF frequency bands, which occurred across muscle contraction levels. To analyze changes in PSDF across the level of muscle contraction, the spectrum was divided into six frequency bandwidths. The analysis of parameters focused on the differences in those parameters between the analyzed muscles, at different levels of muscle contraction.The study revealed that, at muscle contraction levels below 5% RMSmax, contraction level influences standard parameters of the EMG signal and that at such levels of muscle contraction every change in muscle contraction level (recruitment of additional MUs) is reflected in PSDF. The frequency band with the lowest dependency on contraction level was 76–140 Hz for which in both muscles no contraction level effect was detected for contraction levels above 5% RMS_max. The reproducibility of the results was very high, since the observations in of the left and right muscles were almost equal. The other factor, which strongly influences PSDF of the EMG signal, is probably the examined muscle structure (muscle morphology, size, function, subcutaneous layer, cross talk). It seems that low frequency bands up to 25 Hz are especially feasible for type of muscle.     

Muscle sound and electromyogram spectrum analysis during exhausting contractions in man

The changes in the soundmyogram (SMG) and electromyogram (EMG) frequency content during exhausting contractions at 20%, 40%, 60% and 80% of the maximal voluntary contraction (MVC) were investigated by the spectral analysis of the SMG and EMG detected from the biceps brachii muscles of 13 healthy men. The root mean squares (rms) of the two signals were also calculated. Throughout contraction the EMG rms always increased while this was true only at 20% MVC for the SMG. A marked decrease was detected at 60% and 80% MVC. With fatigue the EMG spectra presented a compression towards the lower frequencies at all exercise intensities. The SMG showed a more complex behaviour with a transient increase in its frequency content, followed by a continuous compression of the spectra, at 60% and 80% MVC, and a nearly stable frequency content at lower contraction intensities. This study suggested that different aspects of the changes in the motor unit's activation strategy at different levels of exhausting contractions can be monitored by SMG and EMG signals.     

Motor unit activity during stereotyped finger tasks and computer mouse work.

Motor unit (MU) activity pattern was examined in the right-hand extensor digitorum communis muscle (EDC) during standardised finger movements simulating actual computer mouse tasks. Intramuscular recordings were performed with a quadripolar needle electrode. Nine women performed four lifts of their right-hand index finger, middle finger or both as well as a number of double clicks. Additionally, the subjects performed contra lateral activity with their left-hand fingers and for three subjects recordings were also obtained during an interview with no physical activity. Besides the expected close coupling of MU activity with finger movement, activity was observed in three different situations with no physical requirements. Attention related activity was found before or after performance of the finger movement task, contra lateral activity in right EDC during left-hand finger tasks, and activity during mental activity without any finger movements involved. A relatively large number of doublet occurrences suggest they are a natural part of the activation pattern during performance of the rapid finger movement required to perform an efficient double click on the computer mouse.     

Motor unit recruitment and EMG power spectra during ramp contractions of a bifunctional muscle.

Surface electromyograms (EMGs) were analysed on the short and long head of the biceps brachii (BBSH and BBLH) during single (F and S) or dual (F+S) flexion and supination tasks. It was confirmed, by the analysis of EMG root-mean-square (RMS) values, that the highest activations of BBSH and BBLH were obtained during a maximal dual task. This study was essentially concerned with the analysis of power spectra data obtained during progressive or ramp contractions (RCs). The shape of the power spectra established during the first second of the RCs differs between F, S and F+S tasks. Differences in mean power frequency (MPF) calculated during RCs would be representative of a recruitment of motor units (MUs) that is, at least partly, task-dependent. In order to compare MPF values calculated from RCs performed under different mechanical conditions (F, S and F+S), MPF-RMS(PSD) relationships have been established (RMS(PSD) being defined as the power spectrum density RMS). Both BBSH and BBLH exhibited initial MPF values higher in supination RC than in flexion RC. Because of plateau values reached at the same level of muscle activation whatever the task performed, the slope of the MPF-RMS(PSD) relationship was lower in S than in F. These results are in favour of MU recruitment that is, at least partly, different in F and in S conditions. Dual submaximal tasks seem to mix the activation of the F and S subpopulations of MUs as revealed by the spectral parameters obtained during F+S ramp contractions. This study could find some implication in the field of muscle rehabilitation or reinforcement.     

Motor unit activation patterns during isometric, concentric and eccentric actions at different force levels.

Motor unit activation patterns were studied during four different force levels of concentric and eccentric actions. Eight male subjects performed concentric and eccentric forearm flexions with the movement range from 100 degrees to 60 degrees in concentric and from 100 degrees to 140 degrees elbow angle in eccentric actions. The movements were started either from zero preactivation or with isometric preactivation of the force levels of 20, 40, 60 and 80% MVC. The subjects were then instructed to maintain the corresponding relative force levels during the dynamic actions. Intramuscular and surface EMG was recorded from biceps brachii muscle. Altogether 28 motoneuron pools were analyzed using the intramuscular spike-amplitude frequency (ISAF) analysis technique of Moritani et al. The mean spike amplitude was lower and the mean spike frequency higher in the isometric preactivation phase than in the consequent concentric and eccentric actions. When the movements started with isometric preactivation the mean spike amplitude increased significantly (P<0.001) up to 80% in isometric and concentric actions but in eccentric actions the increase continued only up to 60% (P<0.01). The mean spike frequency in isometric preactivation and in concentric action with preactivation was lower only at the 20% force level (P<0.01) as compared to the other force levels while in eccentric action with preactivation the increase between the force levels was significant (P<0.01) up to 60%. When the movement was started without preactivation the mean spike amplitude at 20% and at 40% force level was higher (P<0.01) in eccentric action than in concentric actions. It was concluded that the recruitment threshold may be lower in dynamic as compared to isometric actions. The recruitment of fast motor units may continue to higher force levels in isometric and in concentric as in eccentric actions which, on the other hand, seems to achieve the higher forces by increasing the firing rate of the active units. At the lower force levels mean spike amplitude was higher in eccentric than in concentric actions which might indicate selective activation of fast motor units. This was, however, the case only when the movements were started without isometric preactivation.     

Motor unit tracking with high-density surface EMG

Following (tracking) individual motor units over time can provide important new insights, both into the relationships among various motor unit (MU) morphological and functional properties and into how these properties are influenced by neuromuscular disorders or interventions. The present study aimed to determine whether high-density surface EMG (HD-sEMG) recordings, which use an array of surface electrodes over a muscle, can increase the yield of MU tracking studies in terms of the number of MUs that can be tracked. For that purpose, four HD-sEMG recording sessions were performed on the thenar muscles of ten healthy subjects. Decomposition of the recorded composite responses yielded a study total of 2849 motor unit action potentials (MUAPs). MUAPs that were found in both of the first two sessions, performed on the same day, were defined as trackable MUAPs. Our results show that 22 (median value; range, 13–34) MUAPs per nerve were trackable, which represented approximately 5% of the total MU population. Of these trackable MUAPs, 16 (11–26) could also be found in one or both of the third and fourth sessions, which were performed between 1 and 13 weeks after the initial studies. Nine (4–18) MUAPs were found in all four sessions. Many of the characteristic MUAP shapes matched well between sessions, even when these sessions were several weeks apart. However, some MUAPs seem very sensitive to changes in arm position or in the muscle’s morphology (e.g., to changes in muscle fiber length due to variable degrees of thumb flexion or extension), particularly those from larger and/or superficial MUs. Standardization is, therefore, essential to detect even small MUAP changes, as may occur with pathology or interventions. If this is accomplished, MU tracking with HD-sEMG may prove to be a powerful tool for a promising type of neurophysiological investigation.     

Relation between electromyogram and force in fatigue

The relationship between the surface electromyogram (SEMG) and force was examined during maximal voluntary contraction (MVC). Isometric MVC of elbow flexors were studied in 18 subjects who performed 27 trials, each consisting of six MVCs lasting 45 s at intervals of 30 s. There was a decrease in the median frequency (Fm) of the SEMG and of the compound action potentials (CAP) during MVC. The CAPs demonstrated that the fall in Fm was associated with a proportional increase in signal power, whereas CAP amplitude did not decrease, indicating intact neuromuscular transmission. The SEMG root-mean-square amplitude remained fairly constant, progressively deviating from force with time of contraction (r = 0.40). When SEMG amplitude was corrected for the Fm change, it tracked force more closely (r = 0.68), indicating a fall in motoneuron drive during MVC. The corrected SEMG was used to calculate the change in the generalized firing rate of motoneurons. The firing rate decreased 60% in the first and sixth contractions, tracked force closely, and corresponded to the firing rate fall seen in late adaptation of motoneurons (r = 0.90, P less than 0.001).