Volume conductor models in surface electromyography: computational techniques

Nonlinear spatio-temporal filter to reduce crosstalk in bipolar electromyogram

Objective.The wide detection volume of surface electromyogram (EMG) makes it prone to crosstalk, i.e. the signal from other muscles than the target one. Removing this perturbation from bipolar recordings is an important open problem for many applications.Approach.An innovative nonlinear spatio-temporal filter is developed to estimate the EMG generated by the target muscle by processing noisy signals from two bipolar channels, placed over the target and the crosstalk muscle, respectively. The filter is trained on some calibration data and then can be applied on new signals. Tests are provided in simulations (considering different thicknesses of the subcutaneous tissue, inter-electrode distances, locations of the EMG channels, force levels) and experiments (from pronator teres and flexor carpi radialis of 8 healthy subjects).Main results.The proposed filter allows to reduce the effect of crosstalk in all investigated conditions, with a statistically significant reduction of its root mean squared of about 20%, both in simulated and experimental data. Its performances are also superior to those of a blind source separation method applied to the same data.Significance.The proposed filter is simple to be applied and feasible in applications in which single bipolar channels are placed over the muscles of interest. It can be useful in many fields, such as in gait analysis, tests of myoelectric fatigue, rehabilitation with EMG biofeedback, clinical studies, prosthesis control.

Investigating the spatial resolution of EMG and MMG based on a systemic multi-scale model

While electromyography (EMG) and magnetomyography (MMG) are both methods to measure the electrical activity of skeletal muscles, no systematic comparison between both signals exists. Within this work, we propose a novel in silico model for EMG and MMG and test the hypothesis that MMG surpasses EMG in terms of spatial selectivity, i.e. the ability to distinguish spatially shifted sources. The results show that MMG provides a slightly better spatial selectivity than EMG when recorded directly on the muscle surface. However, there is a remarkable difference in spatial selectivity for non-invasive surface measurements. The spatial selectivity of the MMG components aligned with the muscle fibres and normal to the body surface outperforms the spatial selectivity of surface EMG. Particularly, for the MMG’s normal-to-the-surface component the influence of subcutaneous fat is minimal. Further, for the first time, we analyse the contribution of different structural components, i.e. muscle fibres from different motor units and the extracellular space, to the measurable biomagnetic field. Notably, the simulations show that for the normal-to-the-surface MMG component, the contribution from volume currents in the extracellular space and in surrounding inactive tissues, is negligible. Further, our model predicts a surprisingly high contribution of the passive muscle fibres to the observable magnetic field.

Motor Nerve Conduction Block Estimation in Demyelinating Neuropathies by Deconvolution

A deconvolution method is proposed for conduction block (CB) estimation based on two compound muscle action potentials (CMAPs) elicited by stimulating a nerve proximal and distal to the region in which the block is suspected. It estimates the time delay distributions by CMAPs deconvolution, from which CB is computed. The slow afterwave (SAW) is included to describe the motor unit potential, as it gives an important contribution in case of the large temporal dispersion (TD) often found in patients. The method is tested on experimental signals obtained from both healthy subjects and pathological patients, with either Chronic Inflammatory Demyelinating Polyneuropathy (CIDP) or Multifocal Motor Neuropathy (MMN). The new technique outperforms the clinical methods (based on amplitude and area of CMAPs) and a previous state-of-the-art deconvolution approach. It compensates phase cancellations, allowing to discriminate among CB and TD: estimated by the methods of amplitude, area and deconvolution, CB showed a correlation with TD equal to 39.3%, 29.5% and 8.2%, respectively. Moreover, a significant decrease of percentage reconstruction errors of the CMAPs with respect to the previous deconvolution approach is obtained (from a mean/median of 19.1%/16.7% to 11.7%/11.2%). Therefore, the new method is able to discriminate between CB and TD (overcoming the important limitation of clinical approaches) and can approximate patients’ CMAPs better than the previous deconvolution algorithm. Then, it appears to be promising for the diagnosis of demyelinating polyneuropathies, to be further tested in the future in a prospective clinical trial.

Nonhomogeneous volume conduction effects affecting needle electromyography: an analytical and simulation study

Objective.Needle electromyography (EMG) is used to study the electrical behavior of myofiber properties in patients with neuromuscular disorders. However, due to the complexity of electrical potential spatial propagation in nonhomogeneous diseased muscle, a comprehensive understanding of volume conduction effects remains elusive. Here, we develop a framework to study the conduction effect of extracellular abnormalities and electrode positioning on extracellular local field potential (LFP) recordings.Methods.The framework describes the macroscopic conduction of electrical potential in an isotropic, nonhomogeneous (i.e. two tissue) model. Numerical and finite element model simulations are provided to study the conduction effect in prototypical monopolar EMG measurements.Results.LFPs recorded are influenced in amplitude, phase and duration by the electrode position in regards to the vicinity of tissue with different electrical properties.Conclusion.The framework reveals the influence of multiple mechanisms affecting LFPs including changes in the distance between the source-electrode and tissue electrical properties.Clinical significance.Our modeled predictions may lead to new ways for interpreting volume conduction effects on recorded EMG activity, for example in neuromuscular diseases that cause structural and compositional changes in muscle tissue. These change will manifest itself by changing the electric properties of the conductor media and will impact recorded potentials in the area of affected tissue.

Inverse modelling to reduce crosstalk in high density surface electromyogram

Surface electromyogram (EMG) has a relatively large detection volume, so that it could include contributions both from the target muscle of interest and from nearby regions (i.e., crosstalk). This interference can prevent a correct interpretation of the activity of the target muscle, limiting the use of surface EMG in many fields. To counteract the problem, selective spatial filters have been proposed, but they reduce the representativeness of the data from the target muscle. A better solution would be to discard only crosstalk from the signal recorded in monopolar configuration (thus, keeping most information on the target muscle). An inverse modelling approach is here proposed to estimate the contributions of different muscles, in order to focus on the one of interest. The method is tested with simulated monopolar EMGs from superficial nearby muscles contracted at different force levels (either including or not model perturbations and noise), showing statistically significant improvements in information extraction from the data. The median over the entire dataset of the mean squared error in representing the EMG of the muscle under the detection electrode was reduced from 11.2% to 4.4% of the signal energy (5.3% if noisy data were processed); the median bias in conduction velocity estimation (from 3 monopolar channels aligned to the muscle fibres) was decreased from 2.12 to 0.72 m/s (1.1 m/s if noisy data were processed); the median absolute error in the estimation of median frequency was reduced from 1.02 to 0.67 Hz in noise free conditions and from 1.52 to 1.45 Hz considering noisy data.

Implementation and experimental validation of surface electromyogram and force model of Tibialis Anterior muscle for examining muscular factors

This study reports a surface electromyogram and force of contraction model. The objective was to investigate the effect of changes in the size, type and number of motor units in the Tibialis Anterior muscle to surface electromyogram and force of dorsiflexion. A computational model to simulate surface electromyogram and associated force of contraction by the Tibialis Anterior muscle was developed. This model was simulated for isometric dorsiflexion, and comparative experiments were conducted for validation. Repeated simulations were performed to investigate the different parameters and evaluate inter-experimental variability. An equivalence statistical test and the Bland–Altman method were used to observe the significance between the simulated and experimental data. Simulated and experimentally recorded data had high similarity for the three measures: maximal power of power spectral density ( p < 0.0001), root mean square of surface electromyogram ( p < 0.0001) and force recorded at the footplate ( p < 0.03). Inter-subject variability in the experimental results was in-line with the variability in the repeated simulation results. This experimentally validated computational model for the surface electromyogram and force of the Tibialis Anterior muscle is significant as it allows the examination of three important muscular factors associated with ageing and disease: size, fibre type and number of motor units.

Modelling the electrical activity of skeletal muscle tissue using a multi-domain approach

Electromyography (EMG) can be used to study the behaviour of the motor neurons and thus provides insights into the physiology of the central nervous system. However, due to the high complexity of neuromuscular control, EMG signals are challenging to interpret. While the exact knowledge of the excitation patterns of a specific muscle within an in vivo experimental setting remains elusive, simulations allow to systematically investigate EMG signals in a controlled environment. Within this context, simulations can provide virtual EMG data, which, for example, can be used to validate and optimise signal analysis methods that aim to estimate the relationship between EMG signals and the output of motor neuron pools. However, since existing methods, which are employed to compute EMG signals, exhibit deficiencies with respect to the physical model itself as well as with respect to numerical aspects, we propose a novel homogenised continuum model that closely resolves the electro-physiological behaviour of skeletal muscle tissue. The proposed model is based on an extension of the well-established bidomain model and includes a biophysically detailed description of the electrical activity within the tissue, which is due to the depolarisation of the muscle fibre membranes. In contrast to all other published EMG models, which assume that the electrical potential field for each muscle fibre can be calculated independently, the proposed model assumes that the electrical potential in the muscle fibres is coupled to the electrical potential in the extracellular space. We show that the newly proposed model is able to simulate realistic EMG signals and demonstrate the potential to employ the predicted virtual EMG signal in order to evaluate the goodness of automated decomposition algorithms.

Multiscale modeling of the neuromuscular system: Coupling neurophysiology and skeletal muscle mechanics

Mathematical models and computer simulations have the great potential to substantially increase our understanding of the biophysical behavior of the neuromuscular system. This, however, requires detailed multiscale, and multiphysics models. Once validated, such models allow systematic in silico investigations that are not necessarily feasible within experiments and, therefore, have the ability to provide valuable insights into the complex interrelations within the healthy system and for pathological conditions. Most of the existing models focus on individual parts of the neuromuscular system and do not consider the neuromuscular system as an integrated physiological system. Hence, the aim of this advanced review is to facilitate the prospective development of detailed biophysical models of the entire neuromuscular system. For this purpose, this review is subdivided into three parts. The first part introduces the key anatomical and physiological aspects of the healthy neuromuscular system necessary for modeling the neuromuscular system. The second part provides an overview on state-of-the-art modeling approaches representing all major components of the neuromuscular system on different time and length scales. Within the last part, a specific multiscale neuromuscular system model is introduced. The integrated system model combines existing models of the motor neuron pool, of the sensory system and of a multiscale model describing the mechanical behavior of skeletal muscles. Since many sub-models are based on strictly biophysical modeling approaches, it closely represents the underlying physiological system and thus could be employed as starting point for further improvements and future developments. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Analytical and Computational Methods > Computational Methods Models of Systems Properties and Processes > Organ, Tissue, and Physiological Models.

Enabling Detailed, Biophysics-Based Skeletal Muscle Models on HPC Systems

Realistic simulations of detailed, biophysics-based, multi-scale models often require very high resolution and, thus, large-scale compute facilities. Existing simulation environments, especially for biomedical applications, are typically designed to allow for high flexibility and generality in model development. Flexibility and model development, however, are often a limiting factor for large-scale simulations. Therefore, new models are typically tested and run on small-scale compute facilities. By using a detailed biophysics-based, chemo-electromechanical skeletal muscle model and the international open-source software library OpenCMISS as an example, we present an approach to upgrade an existing muscle simulation framework from a moderately parallel version toward a massively parallel one that scales both in terms of problem size and in terms of the number of parallel processes. For this purpose, we investigate different modeling, algorithmic and implementational aspects. We present improvements addressing both numerical and parallel scalability. In addition, our approach includes a novel visualization environment which is based on the MegaMol framework and is capable of handling large amounts of simulated data. We present the results of a number of scaling studies at the Tier-1 supercomputer HazelHen at the High Performance Computing Center Stuttgart (HLRS). We improve the overall runtime by a factor of up to 2.6 and achieve good scalability on up to 768 cores.

Motor Unit Action Potential Clustering-Theoretical Consideration for Muscle Activation during a Motor Task

During dynamic or sustained isometric contractions, bursts of muscle activity appear in the electromyography (EMG) signal. Theoretically, these bursts of activity likely occur because motor units are constrained to fire temporally close to one another and thus the impulses are "clustered" with short delays to elicit bursts of muscle activity. The purpose of this study was to investigate whether a sequence comprised of "clustered" motor unit action potentials (MUAP) can explain spectral and amplitude changes of the EMG during a simulated motor task. This question would be difficult to answer experimentally and thus, required a model to study this type of muscle activation pattern. To this end, we modeled two EMG signals, whereby a single MUAP was either convolved with a randomly distributed impulse train (EMG-rand) or a "clustered" sequence of impulses (EMG-clust). The clustering occurred in windows lasting 5-100 ms. A final mixed signal of EMG-clust and EMG-rand, with ratios (1:1-1:10), was also modeled. A ratio of 1:1 would indicate that 50% of MUAP were randomly distributed, while 50% of "clustered" MUAP occurred in a given time window (5-100 ms). The results of the model showed that clustering MUAP caused a downshift in the mean power frequency (i.e., ~30 Hz) with the largest shift occurring with a cluster window of 10 ms. The mean frequency shift was largest when the ratio of EMG-clust to EMG-rand was high. Further, the clustering of MUAP also caused a substantial increase in the amplitude of the EMG signal. This model potentially explains an activation pattern that changes the EMG spectra during a motor task and thus, a potential activation pattern of muscles observed experimentally. Changes in EMG measurements during fatiguing conditions are typically attributed to slowing of conduction velocity but could, per this model, also result from changes of the clustering of MUAP. From a clinical standpoint, this type of muscle activation pattern might help describe the pathological movement issues in people with Parkinson's disease or essential tremor. Based on our model, researchers moving forward should consider how MUAP clustering influences EMG spectral and amplitude measurements and how these changes influence movements.

Inter-day reliability of surface electromyography recordings of the lumbar part of erector spinae longissimus and trapezius descendens during box lifting

Background

Low back pain and neck-shoulder pain are the most reported types of work-related musculoskeletal disorders, and performing heavy lifting at work and working with trunk rotation increase the risk of developing work-related musculoskeletal disorders. Surface electromyography (sEMG) provides information about the electrical activity of muscles. Thus it has the potential to retrieve indirect information about the physical exposure of specific muscles of workers during their actual work. This study aimed to investigate the inter-day reliability of absolute and normalized amplitude of sEMG measurements obtained during repeated standardized reference lifts.

Methods

The inter-day reliability of sEMG of the erector spinae longissimus and trapezius descendens muscles was tested during standardized box lifts. The lifts were performed with loads of 3, 15 and 30 kg from floor to table and from table to table in three conditions, i.e., forearm length (short reaching distance), ¾ arm length (long reaching distance) and forearm length with trunk rotation. Absolute and normalized root mean square (absRMS and normRMS) values were extracted. In line with the guidelines for reporting reliability and agreement studies, we reported relative and absolute reliability estimated by intra class correlation (ICC_3,K), standard error of measurement (SEM) and minimal detectable change in percent (MDC).

Results

The ICC_3,K was higher for absRMS compared with normRMS while SEM and maximal voluntary contraction (MVC) were similar. A total of 50 out of 56, i.e., 89%, and 41 out of 56, i.e., 73%, of the lifting situations were in the range from moderate to almost perfect for absRMS and normRMS, respectively. The SEM and MDC shoved more variation in the lifting situations performed from floor to table and in the trapezius descendens muscle than in the erector spinae longissimus muscle.

Conclusion

This reliability study showed that maximum absRMS and normRMS were found to have a fair to substantial relative inter-day reliability for most lifts but were more reliable when lifting from table to table than from floor to table for both trapezius descendens and erector spinae muscles. The relative inter-day reliability was higher for absolute compared with normalized sEMG amplitudes while the absolute reliability was similar.

Predicting electromyographic signals under realistic conditions using a multiscale chemo-electro-mechanical finite element model

This paper presents a novel multiscale finite element-based framework for modelling electromyographic (EMG) signals. The framework combines (i) a biophysical description of the excitation-contraction coupling at the half-sarcomere level, (ii) a model of the action potential (AP) propagation along muscle fibres, (iii) a continuum-mechanical formulation of force generation and deformation of the muscle, and (iv) a model for predicting the intramuscular and surface EMG. Owing to the biophysical description of the half-sarcomere, the model inherently accounts for physiological properties of skeletal muscle. To demonstrate this, the influence of membrane fatigue on the EMG signal during sustained contractions is investigated. During a stimulation period of 500 ms at 100 Hz, the predicted EMG amplitude decreases by 40% and the AP propagation velocity decreases by 15%. Further, the model can take into account contraction-induced deformations of the muscle. This is demonstrated by simulating fixed-length contractions of an idealized geometry and a model of the human tibialis anterior muscle (TA). The model of the TA furthermore demonstrates that the proposed finite element model is capable of simulating realistic geometries, complex fibre architectures, and can include different types of heterogeneities. In addition, the TA model accounts for a distributed innervation zone, different fibre types and appeals to motor unit discharge times that are based on a biophysical description of the α motor neurons.

Influence of different geometric representations of the volume conductor on nerve activation during electrical stimulation

Volume conductor models with different geometric representations, such as the parallel layer model (PM), the cylindrical layer model (CM), or the anatomically based model (AM), have been employed during the implementation of bioelectrical models for electrical stimulation (FES). Evaluating their strengths and limitations to predict nerve activation is fundamental to achieve a good trade-off between accuracy and computation time. However, there are no studies aimed at clarifying the following questions. (1) Does the nerve activation differ between CM and PM? (2) How well do CM and PM approximate an AM? (3) What is the effect of the presence of blood vessels and nerve trunk on nerve activation prediction? Therefore, in this study, we addressed these questions by comparing nerve activation between CM, PM, and AM models by FES. The activation threshold was used to evaluate the models under different configurations of superficial electrodes (size and distance), nerve depths, and stimulation sites. Additionally, the influences of the sciatic nerve, femoral artery, and femoral vein were inspected for a human thigh. The results showed that the CM and PM had a high error rate, but the variation of the activation threshold followed the same tendency for electrode size and interelectrode distance variation as AM.