Theoretical analysis of surface EMG in voluntary isometric contraction.

A new stochastic model of the surface EMG is suggested and the spectral density of the surface EMG is studied theoretically and experimentally to confirm the validity of this model. Theoretical results show that while the contraction level is not so high, the shape of the spectral density (distribution) does not change and its amplitude is directly proportional to the motor unit firing frequency and recruitment. To illustrate the theoretical results, experiments were carried out for rectus femoris and biceps brachii. The surface EMG was lead off by bipolar surface electrodes. And the spectral density of the surface EMG was calculated using FFT algorithm. From these experimental results, it was confirmed that our theoretical results were almost valid.     

Intramuscular pressure and EMG relate during static contractions but dissociate with movement and fatigue.

Intramuscular pressure (IMP) and electromyography (EMG) mirror muscle force in the nonfatigued muscle during static contractions. The present study explores whether the constant IMP-EMG relationship with increased force may be extended to dynamic contractions and to fatigued muscle. IMP and EMG were recorded from shoulder muscles in three sessions: 1). brief static arm abductions at angles from 0 to 90 degrees, with and without 1 kg in the hands; 2). dynamic arm abductions at angular velocities from 9 to 90 degrees /s, with and without 1 kg in the hands; and 3). prolonged static arm abduction at 30 degrees for 30 min followed by recovery. IMP and EMG increased in parallel with increasing shoulder torque during brief static tasks. During dynamic contractions, peak IMP and EMG increased to values higher than those during static contractions, and EMG, but not IMP, increased significantly with speed of abduction. In the nonfatigued supraspinatus muscle, a linear relationship was found between IMP and EMG; in contrast, during fatigue and recovery, significant timewise changes of the IMP-to-EMG ratio occurred. The results indicate that IMP should be included along with EMG when mechanical load sharing between muscles is evaluated during dynamic and fatiguing contractions.     

Visual feedback on the maximal momentum of the voluntary tibio-tarsal dorsal flexion in hemiplegic patients

In both sides (healthy and plegic) of an hemiplegic group the evaluation of tibio-tarsic isometric torque (with and without the aid of a visual feed-back) was carried out. No significative differences in the "total biological work" during a series of ten contractions, in the two cases were found. A marked increase of the pin-point single contraction in the hemiplegic side was noted. A smoothing in the decrement of the total biological work during serial contractions was also found.     

Inter-electrode spacing of surface EMG sensors: reduction of crosstalk contamination during voluntary contractions

We investigated the influence of inter-electrode spacing on the degree of crosstalk contamination in surface electromyographic (sEMG) signals in the tibialis anterior (target muscle), generated by the triceps surae (crosstalk muscle), using bar and disk electrode arrays. The degree of crosstalk contamination was assessed for voluntary constant-force isometric contractions and for dynamic contractions during walking. Single-differential signals were acquired with inter-electrode spacing ranging from 5 mm to 40 mm. Additionally, double differential signals were acquired at 10 mm spacing using the bar electrode array. Crosstalk contamination at the target muscle was expressed as the ratio of the detected crosstalk signal to that of the target muscle signal. The crosstalk contamination ratio approached a mean of 50% for the 40 mm spacing for triceps surae muscle contractions at 80% MVC and tibialis anterior muscle contractions at 10% MVC. For single differential recordings, the minimum crosstalk contamination was obtained from the 10 mm spacing. The results showed no significant differences between the bar and disk electrode arrays. During walking, the crosstalk contamination on the tibialis anterior muscle reached levels of 23% for a commonly used 22 mm spacing single-differential disk sensor, 17% for a 10 mm spacing single-differential bar sensor, and 8% for a 10 mm double-differential bar sensor. For both studies the effect of electrode spacing on crosstalk contamination was statistically significant. Crosstalk contamination and inter-electrode spacing should therefore be a serious concern in gait studies when the sEMG signal is collected with single differential sensors. The contamination can distort the target muscle signal and mislead the interpretation of its activation timing and force magnitude.     

Repeatability of surface EMG variables during voluntary isometric contractions of the biceps brachii muscle.

The repeatability of initial value and rate of change of mean spectral frequency (MNF), average rectified values (ARV) and muscle fiber conduction velocity (CV) was investigated in the dominant biceps brachii of ten normal subjects during sustained isometric voluntary contractions. Four levels of contraction were studied: 10%, 30%, 50% and 70% of the maximal voluntary contraction level (MVC). Each contraction was repeated three times in each of three different days for a total of nine contractions/level/subject and 90 contractions per level across the ten subjects. Repeatability was investigated using the Intraclass Correlation Coefficient (ICC) and the standard error of the mean (SEM) of the estimates for each subject. Contrary to observations in other muscles, CV estimates appeared to be very repeatable both within and between subjects. CV showed a small but significant increase when contraction force increased from 10% to 50% MVC but no change for further increase of force. As force increased, MNF showed a slight decrease possibly related to a wider spreading of the CV values. The rate of time decrement of MNF and CV increased with the level of contraction. The normalized decrement (% of initial value per second) was in general higher for MNF than for CV and was more repeatable between subjects at 10% MVC than at 70% MVC. A final observation is that a resting time of 5 minutes may not be sufficient after a contraction at 50% or 70% MVC.     

Repeatability of maximal voluntary force and of surface EMG variables during voluntary isometric contraction of quadriceps muscles in healthy subjects

The repeatability of initial values and rate of change of EMG signal mean spectral frequency (MNF), average rectified values (ARV), muscle fiber conduction velocity (CV) and maximal voluntary contraction (MVC) was investigated in the vastus medialis obliquus (VMO) and vastus lateralis (VL) muscles of both legs of nine healthy male subjects during voluntary, isometric contractions sustained for 50 s at 50% MVC. The values of MVC were recorded for both legs three times on each day and for three subsequent days, while the EMG signals have been recorded twice a day for three subsequent days. The degree of repeatability was investigated using the Fisher test based upon the ANalysis Of VAriance (ANOVA), the Standard Error of the Mean (SEM) and the Intraclass Correlation Coefficient (ICC).

Data collected showed a high level of repeatability of MVC measurement (normalized SEM from 1.1% to 6.4% of the mean). MNF and ARV initial values also showed a high level of repeatability (ICC>70% for all muscles and legs except right VMO). At 50% MVC level no relevant pattern of fatigue was observed for the VMO and VL muscles, suggesting that other portions of the quadriceps might have contributed to the generated effort. These observations seem to suggest that in the investigation of muscles belonging to a multi-muscular group at submaximal level, the more selective electrically elicited contractions should be preferred to voluntary contractions.     

Toward an objective interpretation of surface EMG patterns: a voluntary response index (VRI).

Individuals with incomplete spinal cord injuries (SCI) retain varying degrees of voluntary motor control. The complexity of the motor control system and the nature of the recording biophysics have inhibited efforts to develop objective measures of voluntary motor control. This paper proposes the definition and use of a voluntary response index (VRI) calculated from quantitative analysis of surface electromyographic (sEMG) data recorded during defined voluntary movement as a sensitive measure of voluntary motor control in such individuals. The VRI is comprised of two numeric values, one derived from the total muscle activity recorded for the voluntary motor task (magnitude), and the other from the sEMG distribution across the recorded muscles (similarity index (SI)). Calculated as a vector, the distribution of sEMG from the test subject is compared to the average vector calculated from sEMG recordings of the same motor task from 10 neurologically intact subjects in a protocol called brain motor control assessment (BMCA). To evaluate the stability of the VRI, a group of five healthy subjects were individually compared to the prototype, average healthy-subject vectors for all of the maneuvers. To evaluate the sensitivity of this method, the VRI was obtained from two SCI subjects participating in other research studies. One was undergoing supported treadmill ambulation training, and the other a controlled withdrawal of anti-spasticity medications. The supported treadmill training patient's VRI, calculated from pre- and post-training BMCA recordings, reflected the qualitative changes in sEMG patterns and functional improvement of motor control. The VRI of the patient followed by serial BMCA during medication withdrawal also reflected changes in the motor control as a result of changes in anti-spasticity medication. To validate this index for clinical use, serial studies using larger numbers of subjects with compromised motor control should be performed.     

Differential contribution of central command to the cardiovascular responses during static exercise of ankle dorsal and plantar flexion in humans

To examine whether central command contributes differently to the cardiovascular responses during voluntary static exercise engaged by different muscle groups, we encouraged healthy subjects to perform voluntary and electrically evoked involuntary static exercise of ankle dorsal and plantar flexion. Each exercise was conducted with 25% of the maximum voluntary force of the right ankle dorsal and plantar flexion, respectively, for 2 min. Heart rate (HR) and mean arterial blood pressure (MAP) were recorded, and stroke volume, cardiac output (CO), and total peripheral resistance were calculated. With voluntary exercise, HR, MAP, and CO significantly increased during dorsal flexion (the maximum increase, HR: 12 ± 2.3 beats/min; MAP: 14 ± 2.0 mmHg; CO: 1 ± 0.2 l/min), whereas only MAP increased during plantar flexion (the maximum increase, 6 ± 2.0 mmHg). Stroke volume and total peripheral resistance were unchanged throughout the two kinds of voluntary static exercise. With involuntary exercise, there were no significant changes in all cardiovascular variables, irrespective of dorsal or plantar flexion. Furthermore, before the force onset of voluntary static exercise, HR and MAP started to increase without muscle contraction, whereas they had no significant changes with involuntary exercise at the moment. The present findings indicate that differential contribution of central command is responsible for the different cardiovascular responses to static exercise, depending on the strength of central control of the contracting muscle.     

Normalization procedures using maximum voluntary isometric contractions for the serratus anterior and trapezius muscles during surface EMG analysis.

The serratus anterior and trapezius muscles are considered to be the only upward rotators of the scapula and are very important for normal shoulder function. A variety of methods have been used to produce a maximum voluntary isometric contraction (MVIC) of these muscles for normalization of EMG data. The purpose of this study was to quantify the surface EMG activity of the serratus anterior muscle and the upper, middle, and lower parts of the trapezius during 9 manual muscle tests performed with maximum effort in 30 subjects. It was found that no one muscle test produced a MVIC for all individuals. Therefore, to perform normalization within each subject, it is suggested that the 2 or 3 tests identified in this study that produce high levels of EMG activity for each muscle be performed. The scapular protraction muscle test that is often used to normalize data for the serratus anterior muscle produced relatively low levels of EMG activity and was not found to be an optimal test. Muscle tests in which an attempt was made to de-rotate the scapula from an upwardly rotated position produced much higher levels of EMG activity in the serratus anterior muscle.     

EMG and kinematics of normal subjects performing trunk flexion/extensions freely in space.

Ten normal subjects performed continuous trunk flexion/extensions (F/E) without any restraining apparatus at free, 3, 2.25 and 1.5 s periods and a fatiguing task of F/E at 1.5 s period during 45 s. Kinematics of the trunk was obtained with bilateral electromyographic (EMG) activity of the erector spinae (three levels), the abdominal oblique muscles and the rectus abdominis muscles. The free period chosen by the subjects was found to vary between 3.05 and 1.47 s. Lateral flexion of the spine was similar in each task but rotation about its longitudinal axis increased as the F/E period shortened. When left and right side EMG signals were grouped by level of recording, a significant difference in activity was found. Subjects who produced the slowest free F/E displayed larger fatigue indexes derived from the EMG signals for some of their back muscles than for other subjects. The flexion/relaxation phenomenon was considered present in a muscle if a level <10% of the maximum signal recorded during extension was detected. The phenomenon was found in >50% of the observations and occurred at a similar angle in each task. Kinematics and several characteristics of the EMG signals of the trunk were statistically independent of the speed of motion.     

Automated estimation of ankle muscle EMG envelopes and resulting plantar-dorsi flexion torque from 64 garment-embedded electrodes uniformly distributed around the human leg

The design of personalized movement training and rehabilitation pipelines relies on the ability of assessing the activation of individual muscles concurrently with the resulting joint torques exerted during functional movements. Despite advances in motion capturing, force sensing and bio-electrical recording technologies, the estimation of muscle activation and resulting force still relies on lengthy experimental and computational procedures that are not clinically viable. This work proposes a wearable technology for the rapid, yet quantitative, assessment of musculoskeletal function. It comprises of (1) a soft leg garment sensorized with 64 uniformly distributed electromyography (EMG) electrodes, (2) an algorithm that automatically groups electrodes into seven muscle-specific clusters, and (3) a EMG-driven musculoskeletal model that estimates the resulting force and torque produced about the ankle joint sagittal plane. Our results show the ability of the proposed technology to automatically select a sub-set of muscle-specific electrodes that enabled accurate estimation of muscle excitations and resulting joint torques across a large range of biomechanically diverse movements, underlying different excitation patterns, in a group of eight healthy individuals. This may substantially decrease time needed for localization of muscle sites and electrode placement procedures, thereby facilitating applicability of EMG-driven modelling pipelines in standard clinical protocols.     

Influence of anatomical,physical, and detection-system parameters on surface EMG

 Many previous studies were focused on the influence of anatomical, physical, and detection-system parameters on recorded surface EMG signals. Most of them were conducted by simulations. Previous EMG models have been limited by simplifications which did not allow simulation of several aspects of the EMG generation and detection systems. We recently proposed a model for fast and accurate simulation of the surface EMG. It characterizes the volume conductor as a non-homogeneous and anisotropic medium, and allows simulation of EMG signals generated by finite-length fibers without approximation of the current-density source. The influence of thickness of the subcutaneous tissue layers, fiber inclination, fiber depth, electrode size and shape, spatial filter transfer function, interelectrode distance, length of the fibers on surface, single-fiber action-potential amplitude, frequency content, and estimated conduction velocity are investigated in this paper. Implications of the results on electrode positioning procedures, spatial filter design, and EMG signal interpretation are discussed. Received: 23 August 2000 / Accepted in revised form: 18 December 2001     

Quantitative EMG analysis to investigate synergistic coactivation of ankle and knee muscles during isokinetic ankle movement. Part 2: time frequency analysis.

Fundamental to intralimb coordination in the lower extremity, ankle-knee synergy induced by motor irradiation has long been employed to secure facilitation of paralyzed muscles. This study, a companion research subsequent to the time amplitude analysis of surface electromyography in part 1, was to investigate the recruitment strategy of irradiated muscles and prime movers during ankle isokinetic contraction at different contraction speeds (30, 60, 120 and 240 degrees/s) with time frequency analysis. The results indicated the recruitment strategies of the major irradiated muscles (ipsilateral rectus femoris/ipsilateral biceps femoris) and prime movers (anterior tibialis/gastrocnemius) were time-dependent and significantly different in terms of the instantaneous median frequency. In general, the prime movers for ankle isokinetic concentric contraction demonstrated a similar recruitment strategy, irrespective of different contraction speeds. This finding is consistent with the idea of generalized motor programs that speed is one of the constraint parameters supplied to motor programs. Nevertheless, the recruitment strategies of the irradiated muscles were highly inconsistent, varying across trials at different contraction speeds, and were not relevant to those of the prime movers. In addition, the recruitment in the irradiated muscles seemly limited to motor units of low threshold, in spite of maximal voluntary contraction of the prime movers.     

Quantitative EMG analysis to investigate synergistic coactivation of ankle and knee muscles during isokinetic ankle movement. Part 1: time amplitude analysis.

Synergy generally refers to the coordinated action of several motor elements to produce a specific motor task, either intentionally or automatically. One example is motor irradiation, a sudden spread of synergistic muscular coactivation resulting from a forceful single joint movement. To investigate this type of synergy pattern, a quantitative EMG approach was employed to characterize explicit neuromuscular synergy in the ankle-knee complex during maximal ankle isokinetic contraction. In the present study, isokinetic ankle contractions, both dorsiflexion and plantarflexion, at four different speeds (30, 60, 120, and 240 degrees/s) were studied in a normal adult population (N=11) to assess synergistic coactivation of the prime movers (tibialis anterior and gastrocnemius) and irradiated muscles (ipsilateral and contralateral rectus femoris and biceps femoris) of the ankle-knee complex. Electromyographic signals were collected with surface EMG electrodes and processed with traditional time-amplitude analysis to examine specific neural control strategies. The data generally supported several empirical assumptions common to neurological facilitation techniques. (1) Motor irradiation to the knee muscles due to ankle muscle isokinetic contraction was strongly directionally dependent. (2) Motor irradiation to the ipsilateral knee muscles due to ankle isokinetic contraction was speed dependent. (3) The prime movers demonstrated a similar control strategy, irrespective of different contraction speeds.     

Simulation analysis of interference EMG during fatiguing voluntary contractions. Part II--changes in amplitude and spectral characteristics.

Capabilities of amplitude and spectral methods for information extraction from interference EMG signals were assessed through simulation and preliminary experiment. Muscle was composed of 4 types of motor units (MUs). Different hypotheses on changes in firing frequency of individual MUs, intracellular action potential (IAP) and muscle fibre propagation velocity (MFPV) during fatigue were analyzed. It was found that changes in amplitude characteristics of interference signals (root mean square, RMS, or integrated rectified value, IEMG) detected by intramuscular and surface electrodes differed. RMS and IEMG of surface detected interference signals could increase even under MU firing rate reduction and without MU synchronisation. IAP profile lengthening can affect amplitude characteristics more significantly than MU firing frequency. Thus, an increase of interference EMG amplitude is unreliable to reflect changes in the neural drive. The ratio between EMG amplitude and contraction response can hardly characterise the so-called 'neuromuscular efficiency'. The recently proposed spectral fatigue indices can be used for quantification of interference EMG signals. The indices are practically insensitive to MU firing frequency. IAP profile lengthening and decrease in MFPV enhanced the index value, while recruitment of fast fatigable MUs reduced it. Sensitivity of the indices was higher than that of indices traditionally used.     

Active stiffness of the ankle in response to inertial and elastic loads.

Effective stiffness of the musculoskeletal system was examined as a function of the characteristics of an external load. Thirteen healthy subjects provided active contraction of the ankle plantarflexion musculature in a neutral ankle posture to support an external load. Musculoskeletal stiffness was computed from kinetic data recorded in response to dorsiflexion/plantarflexion perturbations. Ankle dynamics were recorded while supporting external loads of 19 and 38 kg with and without antagonistic co-contraction. External loads were applied using pure gravitational mass. In separate trials external loads were applied from stretch of steel springs in parallel with the plantarflexion musculature that also provided added parallel stiffness to the system. Adding external stiffness of 4.9 and 8.1 kN/m surprisingly failed to significantly change the stiffness of the ankle-plus-spring system. This suggests contributions from intrinsic muscle stiffness and reflex stiffness declined in response to added external stiffness. This could not be explained by load magnitudes, ankle postures, or co-activation as these were similar between the inertial and elastic loading conditions. However, non-linear parametric analyses suggest mean intrinsic stiffness of 35.5 kN/m and reflex gain of 11.6 kN/m with a constant reflex delay of 70 ms accurately described the empirical results. The phase response between the mechanical dynamics of the musculoskeletal system and delayed neuromotor feedback combine to provide robust control of system behavior.     

Using two-dimensional spatial information in decomposition of surface EMG signals.

Recently, high-density surface EMG electrode grids and multi-channel amplifiers became available for non-invasive recording of human motor units (MUs). We present a way to decompose surface EMG signals into MU firing patterns, whereby we concentrate on the importance of two-dimensional spatial differences between the MU action potentials (MUAPs). Our method is exemplified with high-density EMG data from the vastus lateralis muscle of a single subject. Bipolar and Laplacian spatial filtering was applied to the monopolar raw signals. From the single recording in this subject six different simultaneously active MUs could be distinguished using the spatial differences between MUAPs in the direction perpendicular to the muscle fiber direction. After spike-triggered averaging, 125-channel two-dimensional MUAP templates were obtained. Template-matching allowed tracking of all MU firings. The impact of spatial information was measured by using subsets of the MUAP templates, either in parallel or perpendicular to the muscle fiber direction. The use of one-dimensional spatial information perpendicular to the muscle fiber direction was superior to the use of a linear array electrode in the longitudinal direction. However, to detect the firing events of the MUs with a high accuracy, as needed for instance for estimation of firing synchrony, two-dimensional information from the complete grid electrode appears essential.     

A comparison of surface and fine wire EMG recordings of gluteus medius during selected maximum isometric voluntary contractions of the hip

Electromyographic (EMG) studies into gluteus medius (GMed) typically involve surface EMG electrodes. Previous comparisons of surface and fine wire electrode recordings in other muscles during high load isometric tasks suggest that recordings between electrodes are comparable when the muscle is contracting at a high intensity, however, surface electrodes record additional activity when the muscle is contracting at a low intensity. The purpose of this study was to compare surface and fine wire recordings of GMed at high and low intensities of muscle contractions, under high load conditions (maximum voluntary isometric contractions, MVICs). Mann–Whitney U tests compared median electrode recordings during three MVIC hip actions; abduction, internal rotation and external rotation, in nine healthy adults. There were no significant differences between electrode recordings in positions that evoked a high intensity contraction (internal rotation and abduction, fine wire activity >77% MVIC; effect size, ES < 0.42; p > 0.277). During external rotation, the intensity of muscle activity was low (4.2% MVIC), and surface electrodes recorded additional myoelectric activity (ES = 0.67, p = 0.002). At low levels of muscle activity during high load isometric tasks, the use of surface electrodes may result in additional myoelectric recordings of GMed, potentially reflective of cross-talk from surrounding muscles.