Surface EMG signal processing during isometric contractions

This paper provides an overview of techniques suitable for the estimation, interpretation and understanding of time variations that affect the surface electromyographic (EMG) signal during sustained voluntary or electrically elicited contractions. These variations concern amplitude variables, spectral variables and muscle fiber conduction velocity, are interdependent and are referred to as the ‘fatigue plot'. The fatigue plot provides information suitable for the classification of muscle behavior. In addition, the information obtainable by means of linear electrode arrays is discussed, and applications of mathematical models for the interpretation of array signals are presented. The model approach provides additional ways for the classification of muscle behavior.     

Non-uniform electromyographic activity during fatigue and recovery of the vastus medialis and lateralis muscles

The aim of the study was to investigate EMG signal features during fatigue and recovery at three locations of the vastus medialis and lateralis muscles. Surface EMG signals were detected from 10 healthy male subjects with six 8-electrode arrays located at 10%, 20%, and 30% of the distance from the medial (for vastus medialis) and lateral (vastus lateralis) border of the patella to the anterior superior spine of the pelvic. Subjects performed contractions at 40% and 80% of the maximal force (MVC) until failure to maintain the target force, followed by 20 2-s contractions at the same force levels every minute for 20 min (recovery). Average rectified value, mean power spectral frequency, and muscle fiber conduction velocity were estimated from the EMG signals in 10 epochs from the beginning of the contraction to task failure (time to task failure, mean ± SD, 70.7 ± 25.8 s for 40% MVC; 27.4 ± 16.8 s for 80% MVC) and from the 20 2 s time intervals during recovery. During the fatiguing contraction, the trend over time of EMG average rectified value depended on location for both muscles (P < 0.05). After 20-min recovery, mean frequency and conduction velocity of both muscles were larger than in the beginning of the fatigue task (P < 0.05) (supernormal values). Moreover, the trend over time of mean frequency during recovery was affected by location and conduction velocity values depended on location for both muscles (P < 0.05). The results indicate spatial dependency of EMG variables during fatigue and recovery and thus the necessity of EMG spatial sampling for global muscle assessment.     

Reproducibility of surface EMG variables in isometric sub-maximal contractions of jaw elevator muscles.

The aims of this study were: (1) to develop and assess reproducibility of a new method for measuring masticatory force in the intercuspal position; (2) to test the reproducibility of surface EMG signal amplitude and spectral variables in constant force contractions of jaw elevator muscles and its dependency on inter-electrode distance. The study was performed on the masseter and temporalis anterior muscles of both sides of nine healthy volunteers. An intraoral compressive-force sensor was used to measure maximal voluntary contraction forces in the intercuspal position and to provide a visual feedback on sub-maximal forces to the subject. Three experimental sessions were performed in three days. In each session, three isometric contractions at 80% of the maximal force were sustained by the subjects for 30s. The intra-class correlation coefficient (ICC) of the maximal force measure was 71.9%. ICC of average rectified value and mean power spectral frequency of the EMG signal increased with inter-electrode distance, with values larger than 70% with 30 mm inter-electrode distance. It was concluded that surface EMG variables measured in isometric contractions of the jaw elevator muscles with the proposed force recording system show good reproducibility for clinical applications when a 30 mm inter-electrode distance is considered.     

Nonlinear surface EMG analysis to detect changes of motor unit conduction velocity and synchronization.

Amplitude and frequency content of the surface electromyographic (EMG) signal reflect central and peripheral modifications of the neuromuscular system. Classic surface EMG spectral variables applied to assess muscle functions are the centroid and median power spectral frequencies. More recently, nonlinear tools have been introduced to analyze the surface EMG; among them, the recurrence quantification analysis (RQA) was shown to be particularly promising for the detection of muscle status changes. The purpose of this work was to analyze the effect of motor unit short-term synchronization and conduction velocity (CV) on EMG spectral variables and two variables extracted by RQA, the percentage of recurrence (%Rec) and determinism (%Det). The study was performed on the basis of a simulation model, which allowed changing the degree of synchronization and mean CV of a number of motor units, and of an experimental investigation of the surface EMG signal properties detected during high-force-level isometric fatiguing contractions of the biceps brachii muscle. Simulations and experimental results were largely in agreement and show that 1) spectral variables, %Rec, and %Det are influenced by CV and degree of synchronization; 2) spectral variables are highly correlated with %Det (R = -0.95 in the simulations and -0.78 and -0.75 for the initial values and normalized slopes, respectively, in the experimental signals), and thus the information they provide on muscle properties is basically the same; and 3) variations of %Det and %Rec in response to changes in muscle properties are significantly larger than the variations of spectral variables. This study validates RQA as a means for fatigue assessment with potential advantages (such as the higher sensitivity to changes of muscle status) with respect to the classic spectral analysis.     

Reproducibility of surface EMG variables and peak torque during three sets of ten dynamic contractions.

The interpretation of the electromyogram (EMG) of dynamic contractions might be difficult because the movement per se introduces additional factors that could affect its characteristics. There is a lack of studies concerning the reproducibility of surface EMG registrations during dynamic contractions. The aim was to investigate the during-the-day reproducibility (using intra-class correlation; ICC) of the peak torque (PT) and the EMG variables (without removing the electrodes) of dynamic contractions. Ten healthy subjects performed three sets of 10 dynamic maximum right-knee extensions with a one-hour interval in between, using an isokinetic dynamometer and the PT was determined. EMG signals were recorded from the right vastus lateralis, rectus femoris and vastus medialis muscles using surface electrodes and the mean frequency of the power spectrum (MNF [Hz]) and the signal amplitude (RMS [microV]), were computed. The ability to relax in-between the maximum extensions was calculated as a ratio of the RMS during the passive flexion phase and the RMS during the active extension phase of each contraction cycle: the signal amplitude ratio (SAR). Both PT (ICC = 0.99) and RMS (ICC = 0.83-0.98) had good reproducibility. The reproducibility of MNF was good for all muscles when the mean of contraction nos.: 1-10 was used. Vastus lateralis had the highest ICC among the three muscles. The reproducibility of SAR was generally poor (ICC < 0.60). The present study showed good reproducibility for common EMG variables (MNF and RMS) obtained during maximum isokinetic contractions.     

Myoelectric manifestations of fatigue in voluntary and electrically elicited contractions

The time course of muscle fiber conduction velocity and surface myoelectric signal spectral (mean and median frequency of the power spectrum) and amplitude (average rectified and root-mean-square value) parameters was studied in 20 experiments on the tibialis anterior muscle of 10 healthy human subjects during sustained isometric voluntary or electrically elicited contractions. Voluntary contractions at 20% maximal voluntary contraction (MVC) and at 80% MVC with duration of 20 s were performed at the beginning of each experiment. Tetanic electrical stimulation was then applied to the main muscle motor point for 20 s with surface electrodes at five stimulation frequencies (20, 25, 30, 35, and 40 Hz). All subjects showed myoelectric manifestations of muscle fatigue consisting of negative trends of spectral variables and conduction velocity and positive trends of amplitude variables. The main findings of this work are 1) myoelectric signal variables obtained from electrically elicited contractions show fluctuations smaller than those observed in voluntary contractions, 2) spectral variables are more sensitive to fatigue than conduction velocity and the average rectified value is more sensitive to fatigue than the root-mean-square value, 3) conduction velocity is not the only physiological factor affecting spectral variables, and 4) contractions elicited at supramaximal stimulation and frequencies greater than 30 Hz demonstrate myoelectric manifestations of muscle fatigue greater than those observed at 80% MVC sustained for the same time.     

Rectus femoris surface myoelectric signal cross-talk during static contractions

The clinical application of EMG requires that the recorded signal is representative of the muscle of interest and is not contaminated with signals from adjacent muscles. Some authors report that surface EMG is not suitable for obtaining information on a single muscle but rather reflects muscle group function [J. Perry, C.S. Easterday, D.J. Antonelli, Surface versus intramuscular electrodes for electromyography of superficial and deep muscles. Physical Therapy 61 (1981) 7–15]. Other authors report however, that surface EMG is adequate to determine individual muscle function, once guidelines pertaining to data acquisition are followed [D.A. Winter, A.J. Fuglevand, S.E. Archer. Cross-talk in surface electromyography: theoretical and practical estimates. Journal of Electromyography and Kinesiology 4 (1994) 15–26]. The aim of this study was to determine whether surface EMG was suitable for monitoring rectus femoris (RF) activity during static contractions. Five healthy subjects, having given written informed consent, participated in this trial. Surface and fine wire EMG from the rectus femoris and the vastus lateralis (VL) muscles were recorded simultaneously during a protocol of static contractions consisting of knee extensions and hip flexions. Ratios were used to quantify the relationship between the surface EMG amplitude value and the fine wire EMG amplitude value for the same contraction. The results showed that hip flexion contractions elicited RF activation only and that knee extension contractions elicited fine wire activity in VL only. When the relationship between RF surface and RF fine wire electrodes was compared for hip flexion and knee extension contractions, it was observed that for all subjects, there was a tendency for increased RF surface activity in the absence of RF fine wire activity during knee extensions. It was concluded that the activity recorded by the RF surface electrode arrangement during knee extension consisted of EMG from the vastii, i.e., cross-talk and that vastus intermedius was the most likely origin of the erroneous signal. Therefore it is concluded that for accurate EMG information from RF, fine wire electrodes are necessary during a range of static contractions.     

Motor unit recruitment strategies investigated by surface EMG variables.

During isometric contractions of increasing strength, motor units (MUs) are recruited by the central nervous system in an orderly manner starting with the smallest, with muscle fibers that usually show the lowest conduction velocity (CV). Theory predicts that the higher the velocity of propagation of the action potential, the higher the power at high frequencies of the detected surface signal. These considerations suggest that the power spectral density of the surface detected electromyogram (EMG) signal may give indications about the MU recruitment process. The purpose of this paper is to investigate the potential and limitations of spectral analysis of the surface EMG signal as a technique for the investigation of muscle force control. The study is based on a simulation approach and on an experimental investigation of the properties of surface EMG signals detected from the biceps brachii during isometric linearly increasing torque contractions. Both simulation and experimental data indicate that volume conductor properties play an important role as confounding factors that may mask any relation between EMG spectral variables and estimated CV as a size principle parameter during ramp contractions. The correlation between spectral variables and CV is thus significantly lower when the MU pool is not stable than during constant-torque isometric contractions. Our results do not support the establishment of a general relationship between spectral EMG variables and torque or recruitment strategy.     

The effects of interelectrode distance over the innervation zone and normalization on the electromyographic amplitude and mean power frequency versus concentric,eccentric, and isometric torque relationships for the vastus lateralis muscle

The purpose of this study was to examine the influence of interelectrode distance (IED) over the estimated innervation zone (IZ) for the vastus lateralis muscle and normalization on the torque-related patterns of responses for electromyographic (EMG) amplitude and mean power frequency (MPF) during concentric isokinetic, eccentric isokinetic, and isometric muscle actions of the leg extensors. Eight men performed submaximal to maximal concentric isokinetic, eccentric isokinetic, and isometric muscle actions of the dominant leg extensors. Surface EMG signals were recorded simultaneously with two bipolar electrode arrangements in single differential configuration (20 and 40 mm IEDs) placed over the estimated IZ for the vastus lateralis muscle and a third electrode arrangement in single differential configuration (20 mm IED) placed distal to the estimated IZ. The results indicated that there were only a few (six of 90 statistical comparisons) significant (p < 0.05) mean differences among the three electrode arrangements for absolute EMG amplitude. There were no mean differences among the three electrode arrangements for absolute or normalized EMG MPF values or normalized EMG amplitude for the three types of muscle actions. Thus, it may be possible to reduce the potential influence of the IZ on amplitude and spectral parameters of the EMG signal through normalization.     

Oxygen availability and motor unit activity in humans.

Six men were studied to determine the interrelationships among blood supply, motor unit (MU) activity and lactate concentrations during intermittent isometric contractions of the hand grip muscles. The subjects performed repeated contractions at 20% of maximal voluntary contraction (MVC) for 2 s followed by 2-s rest for 4 min with either unhindered blood circulation or arterial occlusion given between the 1st and 2nd min. The simultaneously recorded intramuscular MU spikes and surface electromyogram (EMG) data indicated that mean MU spike amplitude, firing frequency and the parameters of surface EMG power spectra (mean power frequency and root mean square amplitude) remained constant during the experiment with unhindered circulation, providing no electrophysiological signs of muscle fatigue. Significant increases in mean MU spike amplitude and frequency were, however, evident during the contractions with arterial occlusion. Similar patterns of significant changes in the surface EMG spectra parameters and venous lactate concentration were also observed, while the integrated force-time curves remained constant. These data would suggest that the metabolic state of the active muscles may have played an important role in the regulation of MU recruitment and rate coding patterns during exercise.     

Electromyography and fatigue during prolonged, low-level static contractions

Findings from five separate studies of EMG changes and muscle fatigue during prolonged low-level static contractions are summarized, and the possible mechanisms behind the changes are briefly discussed. Sustained static contractions (10%, 7% and 5% MVC) of up to 1 h duration were performed by finger flexors, elbow flexors and extensors, and knee extensors. In one experiment, intermittent static arm pulling (triceps) (10 s contraction and 5 s rest, average work load 14% and 10% MVC) was performed for 7 h. The endurance time for the sustained contractions was around one hour for 10% MVC, and it was shown--all in all--that the concept of "indefinite" endurance times at contractions below 15-20% MVC cannot be maintained. After 5% MVC sustained contractions for one hour a 12% reduction in MVC was seen, and significant increases in EMG amplitude and decreases in the mean spectral frequency of the EMG-power spectrum were found. Marked differences were also seen in the EMG changes in the elbow flexors and extensors, and transcutaneous electrical stimulation of the knee extensors showed that low frequency fatigue was present after the contraction. With intermittent contractions similar changes in the EMG parameters were seen after 2-3 h of contractions at 14% MVC. On average, during contractions of 10% MVC no EMG changes were detected. Increased extracellular potassium concentration in the contracting muscles is suggested as a possible explanation of these findings.     

Criterion validation of surface EMG variables as fatigue indicators using peak torque: a study of repetitive maximum isokinetic knee extensions.

A number of studies have been published that have used variables of the electromyogram (EMG) power spectrum during dynamic exercise. Despite these studies there is a shortage of studies of the validity of surface EMG registrations during repetitive dynamic contractions with respect to fatigue. The aim of this study was to investigate if the surface EMG variables mean frequency (MNF [Hz]) and the signal amplitude (RMS [microV]) are valid indicators of muscular fatigue (defined as "any exercise-induced reduction in the capacity to generate force or power output") during maximum repeated isokinetic knee extensions (i.e. criterion validity using peak torque). Twenty-one healthy volunteers performed 100 isokinetic knee extensions at 90 degrees s(-1). EMG signals were recorded from the vastus lateralis, the rectus femoris and the vastus medialis of the right thigh by surface electrodes. MNF and RMS of the EMG together with peak torque (PT [Nm]) were determined for each contraction. MNF showed consequently higher correlation coefficients with PT than RMS did. Positive correlations generally existed between MNF and PT. The majority of the subjects had positive correlations between RMS and PT (i.e. decreases both in PT and in RMS).In conclusion, at the individual level MNF generally - in contrast to RMS - showed good criterion validity with respect to biomechanical fatigue during dynamic maximum contractions.     

Low-level activity of the trunk extensor muscles causes electromyographic manifestations of fatigue in absence of decreased oxygenation

This study was designed to determine whether trunk extensor fatigue occurs during low-level activity and whether this is associated with a drop in muscle tissue oxygenation. Electromyography (EMG) feedback was used to impose constant activity in a part of the trunk extensor muscles. We hypothesized that electromyographic manifestations of fatigue and decreased oxygenation would be observed at the feedback site and that EMG activity at other sites would be more variable without fatigue manifestations. Twelve volunteers performed 30-min contractions at 2% and 5% of the maximum EMG amplitude (EMGmax) at the feedback site. EMG was recorded from six sites over the lumbar extensor muscles and near-infrared spectroscopy was used to measure changes in oxygenation at the feedback site (left L3 level, 3 cm paravertebral). In both conditions, mean EMG activity was not significantly different between electrode sites, whereas the coefficient of variation was lower at the feedback site compared to other recording sites. The EMG mean power frequency (MPF) decreased consistently at the feedback site only. At 5% EMGmax, the decrease in MPF was significant at the group level at all sites ipsilateral to the feedback site. These results suggest that the limited variability of muscle activity at the EMG feedback site and at ipsilateral locations enhances fatigue development. No decreases in tissue oxygenation were detected. In conclusion, even at mean activity levels as low as 2% EMGmax, fatigue manifestations were found in the trunk extensors. These occured in absence of changes in oxygenation of the muscle tissue.     

Modeling nonlinear errors in surface electromyography due to baseline noise: A new methodology

The surface electromyographic (EMG) signal is often contaminated by some degree of baseline noise. It is customary for scientists to subtract baseline noise from the measured EMG signal prior to further analyses based on the assumption that baseline noise adds linearly to the observed EMG signal. The stochastic nature of both the baseline and EMG signal, however, may invalidate this assumption. Alternately, “true” EMG signals may be either minimally or nonlinearly affected by baseline noise. This information is particularly relevant at low contraction intensities when signal-to-noise ratios (SNR) may be lowest. Thus, the purpose of this simulation study was to investigate the influence of varying levels of baseline noise (approximately 2–40% maximum EMG amplitude) on mean EMG burst amplitude and to assess the best means to account for signal noise. The simulations indicated baseline noise had minimal effects on mean EMG activity for maximum contractions, but increased nonlinearly with increasing noise levels and decreasing signal amplitudes. Thus, the simple baseline noise subtraction resulted in substantial error when estimating mean activity during low intensity EMG bursts. Conversely, correcting EMG signal as a nonlinear function of both baseline and measured signal amplitude provided highly accurate estimates of EMG amplitude. This novel nonlinear error modeling approach has potential implications for EMG signal processing, particularly when assessing co-activation of antagonist muscles or small amplitude contractions where the SNR can be low.     

Correlations between short-time Fourier- and continuous wavelet transforms in the analysis of localized back and hip muscle fatigue during isometric contractions.

The aims of the current study were to examine the stationarities of surface electromyographic (EMG) signals obtained from eight bilateral back and hip muscles during a modified Biering-Sørensen test, and to investigate whether short-time Fourier (STFT) and continuous wavelet transforms (CWT) provided similar information with regard to EMG spectral parameters in the analysis of localized muscle fatigue. Twenty healthy subjects participated in the study after giving their informed consent. Reverse arrangement tests showed that 91.6% of the EMG signal epochs demonstrated no significant trends (all p > 0.05), meaning 91.6% of the EMG signal epochs could be considered as stationary signals. Pearson correlation coefficients showed that STFT and CWT in general provide similar information with respect to the EMG spectral variables during isometric back extensions, and as a consequence STFT can still be used.     

Reliability of surface EMG during sustained contractions of the quadriceps.

The purpose of this study was to determine test-retest reliability for median frequency (MDF) and amplitude of surface EMG during sustained fatiguing contractions of the quadriceps. Twenty-two healthy subjects (11 males and 11 females) were tested on two days held one week apart. Surface EMG was recorded from rectus femoris (RF), vastus lateralis (VL) and vastus medialis (VM) during sustained isometric contractions at 80% and 20% of maximal voluntary contraction (MVC) held to exhaustion. Quadriceps fatigue was described using four measures for both MDF and amplitude of EMG: initial, final, normalized final and slope. For both MDF and amplitude, the initial, final and normalized EMG showed moderate to high reliability for all three muscle groups at both contraction levels (ICC=0.59-0.88 for MDF; ICC=0.58-0.99 for amplitude). Slope of MDF and amplitude was associated with a large degree of variability and low ICCs for the 80% but not the 20% MVC. MDF and amplitude of EMG during sustained contractions of the quadriceps are reproducible; normalized final values of MDF and amplitude show better reliability than slope.     

Are the myoelectric manifestations of fatigue distributed regionally in the human medial gastrocnemius muscle?

Myoelectric fatigue typically manifests as variations in the amplitude and spectrum of surface electromyograms (EMGs). Interestingly, these variations seem to be represented locally in different muscles. In this study, we ask whether such a regional distribution of myoelectric fatigue extends to the medial gastrocnemius (MG) muscle. If the MG muscle is activated locally during fatiguing contractions, or if the most fatigable MG fibers are located at distinct muscle regions, then, the myoelectric manifestations of MG fatigue are expected to appear locally in a grid of surface electrodes. With a matrix of surface electrodes (7 × 15 single-differential EMGs) we show that myoelectric fatigue, indeed, manifests regionally in the MG muscle of 12 subjects, who exerted intermittent, fatiguing plantar flections at 50% of their maximal effort. Contrary to the root mean square amplitude, the median frequency of surface EMGs varied consistently across subjects throughout the plantar flections (P = 0.002). On average, changes in EMG spectrum were represented at 78–93 (interquartile interval) out of the 105 channels in the matrix, though with different degrees across channels. For all participants, about 29% of the channels detected significantly greater reductions in median frequency when compared to all channels in the matrix (P < 0.003). Strikingly, these channels were not sparsely distributed; they rather occupied localized skin regions across subjects. Physiologically, our results suggest that, during sub-maximal fatiguing tasks, myoelectric manifestations of MG fatigue are represented in spatially localized muscle regions. Technically, the possibility of studying myoelectric fatigue in the MG muscle appears to depend on the electrode location.     

The relationship between contraction and relaxation during fatiguing isokinetic shoulder flexions. An electromyographic study

Knowledge of the strength, endurance and coordination of the shoulder muscles during dynamic contractions in healthy women would contribute to the understanding of symptoms in that part of the body in patients with myalgia. Twenty clinically healthy women performed single maximal forward shoulder flexions at four different angular velocities (0.57-3.14 rad.s-1). The same subjects also took part in two endurance tests (at angular velocities of 0.57 and 2.09 rad.s-1, respectively) consisting of 150 repeated maximal shoulder flexions. Electromyographic activity (EMG) was registered from four shoulder flexors using surface electrodes. Work was used as the mechanical variable. During the endurance tests subjects rated their perception of fatigue in the shoulder muscles. Work and the amplitude of the EMG signals decreaesd with angular velocity. The mean power frequency of the EMG was constant in the span of angular velocities investigated. During the endurance tests, work and the mean power frequency decreased during the initial 40-60 contractions followed by stable levels. The relative work level was higher at 2.09 than at 0.57 rad.s-1. Greater relative increases of the signal amplitudes of EMG occurred at 2.09 than at 0.57 rad.s-1. The EMG activity between the flexions (during the supposed passive extension) was higher at 2.09 than at 0.57 rad.s-1. Such a high activity was associated with a low mechanical performance at 2.09 rad.s-1. It is suggested that the initial sharp decreases in work and in mean power frequency reflect the fatiguing of the fast twitch motor units. Dynamic work consisting of continuous activity could predispose to muscle complaints.     

Adaptive filtering for ECG rejection from surface EMG recordings.

Surface electromyograms (EMG) of back muscles are often corrupted by electrocardiogram (ECG) signals. This noise in the EMG signals does not allow to appreciate correctly the spectral content of the EMG signals and to follow its evolution during, for example, a fatigue process. Several methods have been proposed to reject the ECG noise from EMG recordings, but seldom taking into account the eventual changes in ECG characteristics during the experiment. In this paper we propose an adaptive filtering algorithm specifically developed for the rejection of the electrocardiogram corrupting surface electromyograms (SEMG). The first step of the study was to choose the ECG electrode position in order to record the ECG with a shape similar to that found in the noised SEMGs. Then, the efficiency of different algorithms were tested on 28 erector spinae SEMG recordings. The best algorithm belongs to the fast recursive least square family (FRLS). More precisely, the best results were obtained with the simplified formulation of a FRLS algorithm. As an application of the adaptive filtering, the paper compares the evolutions of spectral parameters of noised or denoised (after adaptive filtering) surface EMGs recorded on erector spinae muscles during a trunk extension. The fatigue test was analyzed on 16 EMG recordings. After adaptive filtering, mean initial values of energy and of mean power frequency (MPF) were significantly lower and higher respectively. The differences corresponded to the removal of the ECG components. Furthermore, classical fatigue criteria (increase in energy and decrease in MPF values over time during the fatigue test) were better observed on the denoised EMGs. The mean values of the slopes of the energy-time and MPF-time linear relationships differed significantly when established before and after adaptive filtering. These results account for the efficacy of the adaptive filtering method proposed here to denoise electrophysiological signals.