- 1School of Aerospace Engineering, Beijing Institute of Technology, Beijing, China
- 2Department of Exercise Sciences, The University of Auckland, Auckland, New Zealand
- 3Sport Biomechanics Center, Institute of Artificial Intelligence in Sports, Capital University of Physical Education and Sports, Beijing, China
- 4Emerging Interdisciplinary Platform for Medicine and Engineering in Sports (EIPMES), Beijing Municipal Education Commission, Beijing, China
Introduction: Cross-country sit-skiers are often individuals with spinal cord injuries, cerebral palsy, or lower limb disabilities, relying heavily on upper limb strength to generate propulsion during skiing. However, frequent shoulder joint movements significantly increase the incidence of shoulder joint disorders. Therefore, quantifying muscle forces during movement is crucial for understanding upper limb force generation patterns. Currently, electromyography (EMG)-driven neuromusculoskeletal (NMSK) models are the predominant method for calculating muscle forces and joint moments. However, this approach heavily depends on the quality and quantity of EMG data. Surface electrodes are typically used to collect activity data from superficial muscles, but during dynamic movements, factors such as skin stretching, sweating, or friction may cause electrode detachment or poor contact, leading to EMG signal acquisition failures or data loss. In this study, we propose a synergy-assisted EMG-driven NMSK model to predict the activation patterns of missing muscles for cross-country sit-skiing double poling.
Methods: This method is based on individualized EMG-driven NMSK models constructed for each participant, incorporating data from 10 muscles. By utilizing the activation data of 9 known muscles, the model predicts the activation of one missing muscle through synergy analysis. For synergy method selection, we systematically compared four approaches: Non-negative Matrix Factorization (NMF), Principal Component Analysis (PCA), Independent Component Analysis (ICA), and Factor Analysis (FA).
Results: The results demonstrated NMF’s superior performance at 5 synergies, accurately predicting any missing muscle activation among 10 muscles (
Conclusion: This finding validates the effectiveness of the proposed method in predicting upper limb muscle activation during coupled shoulder and elbow joint movements.
1 Introduction
Cross-country sit-skiing, as a high-intensity aerobic endurance sport, was officially included in the competition system at the 4th winter Paralympic Games held in Innsbruck, Austria, in 1988 (Gastaldi et al., 2012). Based on the skier’s ability to control their trunk and pelvis, the International Ski and Snowboard Federation (FIS) classifies athletes into five levels (LW10 to LW12), ranging from low to high. Athletes of all levels compete in the same race, with final rankings determined by adjusting their actual race times using level-specific time factors (FIS, 2024). Cross-country sit-skiers, commonly individuals with spinal cord injuries, growth defects, or cerebral palsy, utilize a sled-based double poling (DP) technique involving trunk flexion and coordinated upper limb muscle activation to generate propulsive force (Ohlsson and Laaksonen, 2017; Rosso et al., 2017; Chen et al., 2023). Existing research has shown that LW10 athletes can generate significant propulsive force during the early phase of the DP cycle by rapidly pressing down with their arms and swinging their upper body (Gastaldi et al., 2016). Liu et al., through an analysis of upper limb isokinetic muscle strength during the DP cycle, further confirmed the critical role of enhancing upper limb strength and coordination in improving athletic performance (Liu et al., 2022). However, studies on upper limb muscle strength in this sport remain relatively scarce. To date, only one study has examined the peak forces of six major muscle groups under different poling camber angles during the DP cycle (Tian et al., 2023).
Methods for quantifying muscle force fall into two primary categories: direct and indirect approaches. Direct measurement methods require surgical implantation of force sensors within human tissue to collect data. Although this approach provides precise measurements, its invasive nature carries inherent risks such as infection and tissue damage, limiting its application primarily to clinical research rather than widespread testing environments (Beidokhti et al., 2017; Trepczynski et al., 2018). Indirect measurement methods utilize neuromusculoskeletal (NMSK) modeling, which integrates computational approaches for calculating muscle activity with musculoskeletal geometries and contact models, to estimate muscle forces and joint moments during movement effectively (Rabbi et al., 2024). A fundamental challenge in this field stems from the anatomical complexity of the human musculoskeletal system: the number of muscles exceeds the skeletal degrees of freedom (DOFs), resulting in muscle redundancy. To resolve this redundancy, researchers have implemented optimization algorithms, notably static optimization (SO) and dynamic optimization (DO), for estimating muscle activation levels (Anderson and Pandy, 2001). These algorithms assume unique muscle force distribution patterns during movement and optimal muscle function. Nevertheless, without experimental electromyography (EMG) data (Manal and Buchanan, 2003) or established muscle contribution ratios (Ackermann and van den Bogert, 2010), these optimization methods cannot generate definitive solutions. Additionally, the predicted muscle activation patterns often fail to accurately represent physiological muscle activity or account for muscle co-contraction. Lloyd and Besier addressed these limitations by proposing an EMG-driven method that has gained widespread acceptance (Lloyd and Besier, 2003). This approach utilizes experimentally measured EMG data and musculotendon unit kinematic data from the NMSK model as inputs. The optimization objective minimizes the discrepancy between joint moments derived from inverse dynamics (ID) calculations and those obtained through EMG-driven model, thereby predicting muscle forces and joint moments (Buchanan et al., 2004). This approach successfully addresses muscle redundancy while facilitating muscle-tendon property calibration. However, the accuracy of muscle force calculations in the EMG-driven model critically depends on the reliability of the collected EMG data. While surface EMG devices are widely employed in biomechanics research due to their non-invasive nature and practical applicability, they cannot capture EMG data from deep muscles that significantly contribute to joint moments. Sartori et al. developed an optimization-assisted EMG-driven model to overcome limitations in EMG data collection from specific muscle groups (Sartori et al., 2014). In this approach, activation signals for muscles with experimentally measured EMG data are fine-tuned during the optimization process, while activation signals for muscles without EMG data are estimated entirely through SO. However, this method lacks robust validation of its predictions for unmeasured muscle activation against experimental data. Moreover, the SO approach, which processes individual time frames independently, may produce unrealistic discontinuities in muscle activation patterns.
To this end, Ao et al. proposed a muscle synergy extrapolation method to estimate missing muscle activation, enhancing the EMG-driven model’s muscle force calculations (Ao et al., 2020). The muscle synergy theory originates from Bernstein’s hypothesis of a dimensionality reduction control strategy, which was later validated and developed into a theory by Bizzi et al. through frog experiments and modeling (Bernshteĭn, 1967; Bizzi et al., 1991; Bizzi et al., 2008). This theory suggests that the central nervous system (CNS) simplifies complex motor control by coordinating the activation of functionally related muscle groups. Building on this theory, Ao et al. reduced the dimensionality of experimental EMG data into time-varying synergy activation coefficients and corresponding time-invariant synergy vectors. The synergy vectors define the contribution weights of each synergy activation coefficient to individual muscle activation. By incorporating the muscle synergy structure into the prediction of muscle activation, this method not only eliminates discontinuities between adjacent time frames but also reduces the number of design variables in the optimization process. In this paper, we refer to these approaches that leverage muscle synergy extrapolation for missing EMG estimation as “synergy-assisted EMG-driven NMSK modeling.” This synergy-assisted EMG-driven NMSK muscle force computation method has been successfully applied to muscle force predictions in gait and upper limb movements (Ao et al., 2022; Li et al., 2023; Tahmid et al., 2024). However, its potential application in scenarios involving complex shoulder and elbow joint movements remains insufficiently validated.
Therefore, this study aims to develop a synergy-assisted EMG-driven NMSK model tailored for the DP technique in cross-country sit-skiing, with the objective of predicting the activation of unmeasured upper limb muscles and selecting the optimal synergy extraction method for model assistance. Specifically, kinematic, kinetic, and surface EMG data were collected from participants during a 30-s maximal effort DP test. Based on these data, an EMG-driven NMSK model was constructed, with individualized parameter calibration and validation performed to ensure accurate muscle force estimation during the DP motion. To address the limitations of experimentally collected EMG data, this study employs muscle synergy analysis to reasonably predict the EMG signals of unmeasured muscles from a limited set of measured muscles.
2 Methods
The following content will introduce a method for calculating muscle force using a synergy-assisted EMG-driven NMSK model from six aspects (Figure 1).

Figure 1. Flowchart of the synergy-assisted EMG-driven NMSK model method. The pink solid box section indicates data collection; the blue dashed box section denotes the NMSK model; the orange dashed box section is the EMG-driven NMSK model; the red solid box area represents the calibration and validation EMG-driven NMSK model; the purple dashed box part represents muscle synergy analysis; and the green solid box part illustrates the optimization process.
2.1 Participants
This study recruited three male college students from the Capital University of Physical Education and Sports as participants (age: 21 ± 2 years, weight: 76.67 ± 5.77 kg, height: 1.84 ± 0.09 m). All participants declare that they have no restrictions or pain in the shoulder, elbow, or trunk joints, indicating the absence of upper limb musculoskeletal disorders. This study was approved by the Ethics Committee of Capital University of Physical Education and Sports (Beijing, China) and conducted in accordance with relevant guidelines. All participants provided written informed consent, including consent for image publication where applicable.
2.2 Experimental data collection and processing
The participants were instructed to perform a 30-s maximum effort test on a self-developed cross-country sit-skiing smart training device (Liu et al., 2020; Liu et al., 2022) (Figure 2A). During this process, kinematic data, EMG signals, and bilateral pole forces were synchronized. Before testing, each participant was required to complete a 10-min warm-up at their preferred pace to familiarize themselves with the equipment. The rolling resistance was set to 5% of the participant’s body weight to simulate real double-poling skiing conditions. During the formal testing phase, participants were instructed to exert maximum effort throughout the entire test to ensure sufficient data collection. A minimum of 14 DP cycles was recorded for each test; tests yielding fewer than this number were discarded (Chen et al., 2023). A 10-min rest period was scheduled between each test to allow participants to recover and prevent muscle fatigue.

Figure 2. Experimental setup and simulation overview. (A) Subject performed DP on a cross-country sit-skiing smart training device. (B) Reflective marker placement and surface EMG setup were applied to the subject. (C) A simplified BUET model was used to simulate the DP, with the green arrows on the hands representing the application of pole forces. Abbreviations: double poling (DP), poling phase (PP), recovery phase (RP).
The three-dimensional trajectories were recorded using an optoelectronic camera system (Qualisys AB, Gothenburg, Sweden) sampling at 200 Hz. A total of 35 reflective markers were utilized, with 2 markers placed on each ski pole and 31 markers on the participants’ body (Figure 2B, for details refer to Chen et al. (Chen et al., 2023)). The kinematic data were filtered using a fourth-order Butterworth low-pass filter with a cutoff frequency of 10 Hz. The marker on the top of the pole was used to differentiate between the poling and recovery phases of the DP cycle. The poling phase (PP) began when the marker reached its highest position and ended at its lowest. Conversely, the recovery phase (RP) started at the marker’s lowest point and concluded when it returned to the highest position (Figure 2C).
Pole force was collected at 50 Hz using a uniaxial gauge load cell (Bengbu Zhongwan Sensor Co., Ltd, China) installed between the pole and the slider. The raw kinetic data were subsequently interpolated and smoothed using Matlab software. Two makers on each ski pole determined the direction of the pole force, with the point of force application defined as the midpoint between the two markers on the wrist, acting on the NMSK model’s hand (Figure 2C).
Surface EMG signals (Cometa Systems Co., Ltd, Italy) were collected with a sampling frequency of 2000 Hz (Figure 2B). Ten muscles (as shown in Table 1) on the right side of the body (subject’s dominant arm) were measured. EMG signals were integrated by a bandpass filter ranging from 50 to 300 Hz, then processed with full-wave rectification, linearly enveloped, and finally normalized. The normalization process involved scaling the EMG amplitudes of each muscle by the maximum EMG activation recorded for that specific muscle across all experimental trials (Bianco et al., 2018; Ao et al., 2020). Then each DP cycle was time-normalized to 100 points, with 50 points corresponding to the PP and the remaining 50 representing the RP.
2.3 Development of EMG-driven NMSK model
The EMG-driven NMSK model comprises four components: the NMSK model for calculating muscle kinematics and reference moments, activation dynamics, contraction dynamics, and forward dynamics, as detailed below.
The NMSK model was based on the OpenSim BUET model developed for cross-country sit-skiing DP propulsion (Chen et al., 2023). The BUET model includes 17 articulating rigid bodies, 35 DOFs, and 472 musculotendon actuators. The glenohumeral and sternoclavicular joints were modeled as three-degrees-of-freedom rotational joint, which simulated the movement of the joints in the coronal, sagittal, and transverse planes. The acromioclavicular joint was established as a weld joint, but it was unable to move. Based on the BUET model, individualized models were scaled for each subject, with inverse kinematics calculations first performed to obtain joint angle during DP propulsion (Supplementary Figure S2), followed by ID and muscle analysis to obtain reference moments at the shoulder and elbow, as well as kinematic data of the muscles (musculotendon length and moment arm) during DP propulsion. Finally, the obtained data underwent time normalization, dividing one DP cycle into 100 points, with 50 designated for PP and the other 50 for RP. To reduce optimization time, we have reduced the number of musculotendon actuators in the BUET model to 124 (Figure 2C).
The activation dynamics model consists of two components: the backward difference model (as shown in Equation 1) and the nonlinear model (as shown in Equation 2) (Lloyd and Besier, 2003). The preprocessed EMG signals are converted to neural excitation using a critically damped, second-order dynamic system, and muscle activation is then computed via a nonlinear transformation. The respective formulations are as follows:
where all variables and parameters are as previously defined (see (Lloyd and Besier, 2003), for details).
The calculation of muscle force in muscle contraction dynamics is based on the following key parameters: maximum isometric force (
Active muscle force-activation-length-velocity relationship as shown in in Table 2 (Thelen, 2003). A Gaussian function was employed to model the active force-length relationship (
After completing the force calculations for contraction dynamics, the joint moments are derived using the principles of forward dynamics. The resultant moment at a single joint is equal to the sum of the products of all muscle forces acting on that joint, their corresponding pennation angle cosines, and the moment arms (details are in Table 2).
Utilizing these four core components, the EMG-driven NMSK model was developed. The dynamic computation process of this model is illustrated in Figure 3, and the algorithm was implemented within the Matlab environment.

Figure 3. Dynamic calculation process of the EMG-driven model. Blue represents the muscle kinematics model; yellow represents the activation dynamics model; black represents the mechanical relationship model between muscle and tendon; green represents the active muscle force-activation-length-velocity relationship mode; orange represents the passive muscle force-length relationship model; purple represents the tendon length-tendon force relationship model; red represents the forward dynamics model.
2.4 Calibration and validation EMG-driven NMSK model
After constructing the EMG-driven NMSK model, it is necessary to calibrate individualized parameters for each subject. Referring to previous research, Kian et al. successfully calibrated the EMG-driven model for shoulder motion tasks by optimizing musculotendon parameters
The calibration process is as follows. Firstly, kinematic data, external force data, and EMG signals from 5 DP cycles were input into the EMG-driven model to calculate the joint moments of the shoulder and elbow. Subsequently, the model-predicted moments were compared with those obtained through ID. A simulated annealing algorithm was employed for iterative optimization, continuously adjusting the model parameters to minimize the root mean square error (RMSE) between the two sets of results (as shown in Equation 4; Table 2). Finally, individualized parameters for each subject were determined and output, ensuring that the model accurately captures inter-individual physiological differences. It is worth noting that for the parameters
Using 5 DP cycles data that differ from the calibration process, the individualized EMG-driven NMSK model of the subjects was validated by comparing the moment results calculated by the ID.
2.5 Synergy-assisted EMG-driven NMSK model
Based on the calibrated and validated subject-specific EMG-driven NMSK model, a synergy-assisted EMG-driven NMSK model was developed (Figure 1). In this study, the EMG data of one muscle were sequentially excluded and treated as unmeasured, while the remaining muscles’ EMG data were treated as measured. The measured EMG data were organized into a matrix
Non-negative Matrix Factorization (NMF) was applied to decompose
In this section, the synergy-assisted EMG-driven NMSK model was evaluated using data from 10 DP cycles, which were the same datasets used for calibrating and validating the EMG-driven NMSK model. The entire process was implemented in MATLAB, utilizing a simulated annealing algorithm to optimize the objective function. The predicted muscle activation values were then compared with the experimentally measured results.
2.6 Statistical analyses
To evaluate the performance of the synergy-assisted EMG-driven NMSK model in predicting unmeasured muscle activation, this study selected 10 complete DP cycles from each participant during the calibration and validation phases of the EMG-driven NMSK model for comprehensive evaluation. First,
3 Results
3.1 Assessment of EMG-driven NMSK model accuracy
Figure 4 illustrated the relationship between the moment calculations of the EMG-driven NMSK model and the experimental joint moments calculated using ID for three subjects during calibration and validation trials. The figure shows that, aside from a slight shortfall in the accuracy of predicting the arm internal-external rotation moment, the predictions for other joint moments demonstrate a high level of accuracy. Further analysis of the data in Table 3 revealed that during the calibration trials, the model performed best in the arm flexion-extension moment (%RMSE = 0.15 ± 0.03,

Figure 4. In both the calibration and validation trials, the inverse dynamics (ID) method and the EMG-driven NMSK model with all muscles were used to calculate the average joint moments for four degrees of freedom at the shoulder and elbow joints of three participants. The blue curve represents the average reference joint moments calculated by the ID method, while the red curve represents the average joint moments calculated by the EMG-driven NMSK model. The 0%–50% of the cycle corresponds to the poling phase (PP) of the double poling (DP), while the 50%–100% corresponds to the recovery phase (RP).

Table 3. The %RMSE,
3.2 Muscle activation prediction
This study employed NMF, PCA, ICA, and FA synergy-assisted EMG-driven NMSK models for muscle activation prediction, with results shown in Figures 5, 6 and Supplementary Figures S3–S9. The findings revealed that the number of synergies had no significant impact on RMSE,

Figure 5. The variation curves of the predicted average space missing muscle activation and the experimentally measured average muscle activation under different synergy number conditions were plotted. The blue curve represents the experimentally measured values, while the red curve corresponds to the predicted values calculated using the NMF synergy-assisted EMG-driven NMSK model. The 0%–50% of the cycle corresponds to the poling phase (PP) of the double poling (DP), while the 50%–100% corresponds to the recovery phase (RP). Muscle abbreviations: AD, anterior deltoid; BB, biceps brachii; BRD, brachioradialis; IF, infraspinatus; LD, latissimus dorsi; MD, middle deltoid; PD, posterior deltoid; PM, pectoralis major; TB, triceps brachii; TMj, teres major.

Figure 6. Overall RMSE,
The NMF synergy-assisted EMG-driven NMSK model demonstrates excellent performance in predicting missing muscle activation (as shown in Figure 5). The predicted muscle activation curves closely resemble the experimental measurement curves, effectively capturing the dynamic changes in muscle activation, particularly reflecting the actual trends during both the rising and falling phases of activation. In the early phase of the DP cycle (within the PP), the model’s predicted curves show a high degree of overlap with the experimental values, indicating its accuracy during the initial activation stage. The variation in the number of synergies affects the prediction of muscle activation. When the number of synergies reaches four or more, the

Figure 7. The
Further analysis using the Kruskal–Wallis test revealed the impact of the number of synergies on the RMSE,

Figure 8. The RMSE,
In addition, we evaluated the ability of the NMF synergy-assisted EMG-driven NMSK model to estimate joint moments. The results showed that changes in the number of synergies did not have a significant effect on the accuracy of joint moment estimation (see Supplementary Figures S10–S13). However, the omission of specific muscles led to significant differences in joint moment estimation (see Supplementary Tables S6–S9). Specifically, the omission of the PD significantly affected the estimation of both arm adduction–abduction and flexion–extension moments. The absence of the TB substantially impacted the estimation of arm adduction–abduction moments as well as elbow flexion–extension moments. Meanwhile, the exclusion of the IF notably influenced the estimation of arm internal–external rotation moments. Notably, the joint moment estimations obtained from the NMF synergy-assisted EMG-driven NMSK model were consistent with those obtained from the EMG-driven NMSK model with all muscles.
4 Discussion
In this study, we developed an EMG-driven NMSK model tailored for the DP technique in cross-country sit-skiing. By comparing the model predictions with experimental data, we validated its ability to accurately calculate the coupled joint moments of the shoulder and elbow during the DP motion. Furthermore, based on this model, we established a synergy-assisted EMG-driven NMSK model to predict the activation patterns of missing muscles during upper limb movements and systematically investigated the impact of synergy number on prediction performance. Regarding synergy analysis methodology, we comparatively evaluated four decomposition approaches: NMF, PCA, ICA, and FA. The results demonstrated that the NMF method significantly outperformed other methods, and when the synergy number was set to 5, the synergy-assisted EMG-driven NMSK model could accurately predict the activation patterns of any missing muscle among the 10 muscles, with low error and moderate-to-strong correlation.
By collecting surface EMG data from 10 upper limb muscles, we constructed and calibrated individualized EMG-driven NMSK models for each participant (see Figure 4). The results demonstrated that the model could accurately calculate shoulder adduction-abduction moment, shoulder flexion-extension moment, and elbow flexion-extension moment. However, the model exhibited limited accuracy in predicting shoulder internal-external rotation moment. This limitation may be attributed to several factors. First, the study did not fully capture the activation signals of all muscles involved in shoulder internal-external rotation, particularly the deep muscles (e.g., teres minor). Second, shoulder internal-external rotation involves multifunctional muscles such as the IF, TMj, PM, and LD, which may prioritize optimizing moments for other DOFs, thereby reducing their contribution to internal-external rotation moment. Additionally, the potential roles of the BB and TB in shoulder internal-external rotation moment were not considered in the model, which may further limit its prediction accuracy for this DOF.
Based on the above model, this study established a synergy-assisted EMG-driven NMSK model to evaluate upper limb muscle activation patterns during the DP motion in cross-country sit-skiing. The study specifically examined how four synergy analysis methods (NMF, PCA, ICA, FA) influence the model’s muscle activation prediction accuracy. Since synergy analysis capability is closely related to the choice of synergy number, systematically evaluating the impact of synergy number on model prediction performance is essential. The results demonstrated that the NMF method consistently achieved optimal muscle activation prediction performance across all tested synergy numbers (Figures 5, 6). This superiority originates from its unique algorithmic properties: by enforcing non-negativity constraints on both basis vectors and activation coefficients, NMF strictly adheres to the physiological characteristics of muscle activation, effectively preventing non-physiological negative values and thereby significantly enhancing the physiological plausibility and interpretability of the results. Moreover, the NMF method exhibits exceptional noise robustness and can accurately identify muscle synergy patterns without imposing strict distributional assumptions (Tresch et al., 2006; Ebied et al., 2018; Zhao et al., 2022). This approach has been widely adopted in sports and upper limb movement analysis, demonstrating reliable predictive performance (Baifa et al., 2021; Wang et al., 2021). Notably, although Tahmid et al.'s study did not conduct systematic comparisons with other synergy methods, their application of NMF for upper limb muscle prediction also yielded excellent results (Tahmid et al., 2024).
Building upon these findings, this study conducted a systematic discussion on the predictive performance of the NMF synergy-assisted EMG-NMSK model. The results showed that when the synergy number was set to 5, the model could predict the activation of any missing muscle among the 10 muscles with low error and moderate-to-strong correlation (Figure 8; Supplementary Table S4). This finding is consistent with the results of Ao et al., who reported optimal prediction performance with a synergy number of 5 or 6 in lower limb gait studies (Ao et al., 2020). The model exhibited high accuracy in predicting shoulder muscle activation patterns, which typically have single-peaked and less complex activation profiles. However, the prediction accuracy for elbow muscles was relatively lower. For example, the BB had an
When the synergy number was set to 5, the model achieved an optimal balance between reconstruction capability and generalization ability, effectively capturing the main features of muscle activation while avoiding overfitting and performance degradation. The reasons for this phenomenon are as follows: First, as the synergy number increases, the model’s ability to reconstruct the original muscle activation matrix improves. However, when the synergy number exceeds a certain threshold, the reconstruction capability of NMF saturates (Santuz et al., 2020). At the same time, the dimensionality reduction capability of synergy analysis diminishes, leading to a significant increase in the number of unknown parameters during optimization. This increase in parameters makes the model more prone to fitting noise or irrelevant details, thereby reducing its generalization ability. This phenomenon explains why RMSE reached its lowest value when the synergy number was 5 and increased when the synergy number was further increased to 6 or 7. Second, the synergy-assisted method uses the dimensionality-reduced activation matrix
Our study utilized a NMF synergy-assisted EMG-driven NMSK model to evaluate joint moment estimation capabilities when predicting the missing of individual muscles. The results indicated that variations in the number of synergies did not significantly affect moment estimation accuracy. However, among the muscles examined, the PD, TB, and IF demonstrated significantly greater impacts on the precision of moment predictions compared to other muscles (Supplementary Tables S3–S6). Specifically, PD, with its substantial moment arms for shoulder extension and abduction, plays a critical role in generating arm adduction-abduction and flexion-extension moments (Hik and Ackland, 2019). TB, characterized by a large maximum isometric force and an important moment arm for elbow extension, significantly influences elbow flexion-extension moments (Itoi et al., 2008). IF, a key component of the rotator cuff, is essential for shoulder stability and internal-external rotation (Hik and Ackland, 2019). These findings align with prior studies by Rabbi et al. and Steel et al., which emphasized that in scenarios involving fewer muscle combinations, dominant synergy muscles or those with greater maximum isometric force substantially affect predictive performance (Steele et al., 2013; Rabbi et al., 2022). Notably, while Ao et al. and Rabbi et al. focused on lower-limb movement studies, their conclusion that selecting muscles with significant joint moment contributions optimizes the accuracy of synergy-based extrapolation for predicting remaining muscle activations and joint forces using limited muscle data aligns with our findings (Ao and Fregly, 2024; Rabbi et al., 2024). Our results provide a foundation for future investigations into predicting the activation patterns of remaining muscles during DP movements with reduced muscle sets, highlighting the critical roles of PD, TB, and IF.
Despite the promising results, this study has certain limitations. First, the model did not incorporate the activation patterns of deep muscles in the shoulder and elbow joints, which may have contributed to moment calculation errors. Second, the BUET model used in this study constrained the sternoclavicular and acromioclavicular joints, without considering their potential influence on the results. Third, the small sample size (n = 3) limits the generalizability of the findings. To address these limitations, future studies will: (1) optimize the model structure to include deep muscle activations; (2) relax joint constraints in the BUET model to investigate their biomechanical effects; (3) explore the feasibility of predicting a larger number of missing muscle activations using fewer known activations; and (4) validate the model with a larger cohort to establish its generalizability.
5 Conclusion
This study developed an EMG-driven NMSK model tailored for the DP technique in cross-country sit-skiing, capable of accurately calculating the joint moments of the shoulder and elbow during coupled movements. Furthermore, to address the issue of missing muscle data in the collected dataset, a synergy-assisted EMG-driven NMSK model was further developed. The results suggest that NMF synergy-assisted method with the synergy number to 5 allows the model to reasonably predict the activation patterns of any missing muscles among the 10 muscles. However, given the limited number of participants and the observed accuracy levels, these findings should be interpreted with caution.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving humans were approved by Ethics Committee of Capital University of Physical Education and Sports. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
XC: Conceptualization, Data curation, Formal Analysis, Methodology, Project administration, Software, Validation, Writing – original draft. ZY: Formal Analysis, Methodology, Writing – original draft. XG: Methodology, Writing – original draft. YZ: Conceptualization, Writing – original draft. CL: Conceptualization, Funding acquisition, Writing – review and editing. BH: Conceptualization, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was funded by the National Natural Science Foundation of China (2102277 [CLL]) and the National Key R&D Program of China (2018YFF0300800 [BH]).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2025.1585127/full#supplementary-material
References
Ackermann, M., and van den Bogert, A. J. (2010). Optimality principles for model-based prediction of human gait. J. Biomechanics 43 (6), 1055–1060. doi:10.1016/j.jbiomech.2009.12.012
Anderson, F. C., and Pandy, M. G. (2001). Static and dynamic optimization solutions for gait are practically equivalent. J. Biomechanics 34 (2), 153–161. doi:10.1016/S0021-9290(00)00155-X
Ao, D., and Fregly, B. J. (2024). Comparison of synergy extrapolation and static optimization for estimating multiple unmeasured muscle activations during walking. J. Neuroengineering Rehabilitation 21 (1), 194. doi:10.1186/s12984-024-01490-y
Ao, D., Shourijeh, M. S., Patten, C., and Fregly, B. J. (2020). Evaluation of synergy extrapolation for predicting unmeasured muscle excitations from measured muscle synergies. Front. Comput. Neurosci. 14, 588943. doi:10.3389/fncom.2020.588943
Ao, D., Vega, M. M., Shourijeh, M. S., Patten, C., and Fregly, B. J. (2022). EMG-driven musculoskeletal model calibration with estimation of unmeasured muscle excitations via synergy extrapolation. Front. Bioeng. Biotechnol. 10, 962959. doi:10.3389/fbioe.2022.962959
Baifa, Z., Xinglong, Z., and Dongmei, L. (2021). Muscle coordination during archery shooting: a comparison of archers with different skill levels. Eur. J. Sport Sci. 23 (1), 54–61. doi:10.1080/17461391.2021.2014573
Beidokhti, H. N., Janssen, D., van de Groes, S., Hazrati, J., Van den Boogaard, T., and Verdonschot, N. (2017). The influence of ligament modelling strategies on the predictive capability of finite element models of the human knee joint. J. Biomechanics 65, 1–11. doi:10.1016/j.jbiomech.2017.08.030
Bianco, N. A., Patten, C., and Fregly, B. J. (2018). Can measured synergy excitations accurately construct unmeasured muscle excitations? J. Biomech. Eng. 140 (1). doi:10.1115/1.4038199
Bizzi, E., Mussa-Ivaldi, F. A., and Giszter, S. (1991). Computations underlying the execution of movement: a biological perspective. Science 253 (5017), 287–291. doi:10.1126/science.1857964
Bizzi, E., Cheung, V. C., d'Avella, A., Saltiel, P., and Tresch, M. (2008). Combining modules for movement. Brain Res. Rev. 57 (1), 125–133. doi:10.1016/j.brainresrev.2007.08.004
Blache, Y., Michaud, B., Rogowski, I., Monteil, K., and Begon, M. (2019). Sensitivity of shoulder musculoskeletal model predictions to muscle-tendon properties. IEEE Trans. Biomed. Eng. 66 (5), 1309–1317. doi:10.1109/TBME.2018.2872833
Buchanan, T. S., Lloyd, D. G., Manal, K., and Besier, T. F. (2004). Neuromusculoskeletal modeling: estimation of muscle forces and joint moments and movements from measurements of neural command. J. Appl. Biomechanics 20 (4), 367–395. doi:10.1123/jab.20.4.367
Chen, X., Huang, Y., Jiang, L., Sun, Q., Tian, Y., Zhou, Z., et al. (2023). Bilateral upper extremity trunk model for cross-country sit-skiing double poling propulsion: model development and validation. Med. & Biol. Eng. & Comput. 61 (2), 445–455. doi:10.1007/s11517-022-02724-8
Ebied, A., Kinney-Lang, E., Spyrou, L., and Escudero, J. (2018). Evaluation of matrix factorisation approaches for muscle synergy extraction. Med. Eng. & Phys. 57, 51–60. doi:10.1016/j.medengphy.2018.04.003
FIS (2024). The international ski competition rules (ICR)[S]. Switzerland: International Ski and Snowboard Federation. Available online at: https://www.fis-ski.com/para-nordic/documents.
Gastaldi, L., Pastorelli, S., and Frassinelli, S. (2012). A biomechanical approach to paralympic cross-country sit-ski racing. Clin. J. Sport Med. 22 (1), 58–64. doi:10.1097/JSM.0b013e31824202d3
Gastaldi, L., Mauro, S., and Pastorelli, S. (2016). Analysis of the pushing phase in paralympic cross-country sit-skiers - class LW10. J. Adv. Res. 7 (6), 971–978. doi:10.1016/j.jare.2016.10.003
Hik, F., and Ackland, D. C. (2019). The moment arms of the muscles spanning the glenohumeral joint: a systematic review. J. Anat. 234 (1), 1–15. doi:10.1111/joa.12903
Itoi, E., Lee, S.-B., Berglund, L. J., Schultz, F. M., Neale, P. G., and An, K.-N. J. J. o.M. R. (2008). Moment arms of the arm muscles at the glenohumeral joint using the tendon excursion method. J. Musculoskelet. Res. 11 (2), 45–53. doi:10.1142/s0218957708001973
Kian, A., Pizzolato, C., Halaki, M., Ginn, K., Lloyd, D., Reed, D., et al. (2021). The effectiveness of EMG-driven neuromusculoskeletal model calibration is task dependent. J. Biomechamics 129, 110698. doi:10.1016/j.jbiomech.2021.110698
Li, J., Wang, R., and Pan, L. (2023). An enhanced EMG-driven musculoskeletal model based on non-negative matrix factorization. Biomed. Signal Process. Control 79, 104178. doi:10.1016/j.bspc.2022.104178
Liu, C., Zhou, Z., Yin, J., Zhou, L., Sun, G., and Wang, Y. (2020). An intelligent training system for simulating cross-country sit-skiing. China.
Liu, C., Tian, Y., Zhou, L., Tian, Z., Sun, G., Yin, J., et al. (2022). Upper limb isokinetic muscle strength predicts the performance in cross-country sit-skiing. Sci. Rep. 12 (1), 6093. doi:10.1038/s41598-022-10103-4
Lloyd, D. G., and Besier, T. F. (2003). An EMG-driven musculoskeletal model to estimate muscle forces and knee joint moments in vivo. J. Biomechanics 36 (6), 765–776. doi:10.1016/s0021-9290(03)00010-1
Ma, Y., Xie, S., and Zhang, Y. (2016). A patient-specific EMG-driven neuromuscular model for the potential use of human-inspired gait rehabilitation robots. Comput. Biol. Med. 70, 88–98. doi:10.1016/j.compbiomed.2016.01.001
Manal, K., and Buchanan, T. S. (2003). A one-parameter neural activation to muscle activation model: estimating isometric joint moments from electromyograms. J. Biomechanics 36 (8), 1197–1202. doi:10.1016/s0021-9290(03)00152-0
Ohlsson, M. L., and Laaksonen, M. S. (2017). Sitting position affects performance in cross-country sit-skiing. Eur. J. Appl. Physiology 117 (6), 1095–1106. doi:10.1007/s00421-017-3596-y
Rabbi, M. F., Diamond, L. E., Carty, C. P., Lloyd, D. G., Davico, G., and Pizzolato, C. (2022). A muscle synergy-based method to estimate muscle activation patterns of children with cerebral palsy using data collected from typically developing children. Sci. Rep. 12 (1), 3599. doi:10.1038/s41598-022-07541-5
Rabbi, M. F., Davico, G., Lloyd, D. G., Carty, C. P., Diamond, L. E., and Pizzolato, C. (2024). Muscle synergy-informed neuromusculoskeletal modelling to estimate knee contact forces in children with cerebral palsy. Biomechanics Model. Mechanobiol. 23 (3), 1077–1090. doi:10.1007/s10237-024-01825-7
Rosso, V., Gastaldi, L., Rapp, W., Lindinger, S., Vanlandewijck, Y., and Linnamo, V. (2017). Biomechanics of simulated versus natural cross-country sit skiing. J. Electromyogr. Kinesiol. 32, 15–21. doi:10.1016/j.jelekin.2016.11.002
Santuz, A., Ekizos, A., Kunimasa, Y., Kijima, K., Ishikawa, M., and Arampatzis, A. (2020). Lower complexity of motor primitives ensures robust control of high-speed human locomotion. Heliyon 6 (10), e05377. doi:10.1016/j.heliyon.2020.e05377
Sartori, M., Farina, D., and Lloyd, D. G. (2014). Hybrid neuromusculoskeletal modeling to best track joint moments using a balance between muscle excitations derived from electromyograms and optimization. J. Biomechanics 47 (15), 3613–3621. doi:10.1016/j.jbiomech.2014.10.009
Steele, K. M., Tresch, M. C., and Perreault, E. J. (2013). The number and choice of muscles impact the results of muscle synergy analyses. Front. Comput. Neurosci. 7, 105. doi:10.3389/fncom.2013.00105
Tahmid, S., Font-Llagunes, J. M., and Yang, J. (2024). Upper extremity muscle activation pattern prediction through synergy extrapolation and electromyography-driven modeling. J. Biomechanical Eng. 146 (1), 011005. doi:10.1115/1.4063899
Thelen, D. G. (2003). Adjustment of muscle mechanics model parameters to simulate dynamic contractions in older adults. J. Biomechanical Eng. 125 (1), 70–77. doi:10.1115/1.1531112
Tian, Y., Chen, X., Liu, Y., Sun, G., Zhou, Z., Liu, C., et al. (2023). Effects of poling camber angle on the biomechanics of cross-country sit-skiing. Sci. Rep. 13 (1), 20893. doi:10.1038/s41598-023-48359-z
Trepczynski, A., Kutzner, I., Schwachmeyer, V., Heller, M. O., Pfitzner, T., and Duda, G. N. (2018). Impact of antagonistic muscle co-contraction on in vivo knee contact forces. J. NeuroEngineering Rehabilitation 15 (1), 101. doi:10.1186/s12984-018-0434-3
Tresch, M. C., Cheung, V. C. K., and d'Avella, A. (2006). Matrix factorization algorithms for the identification of muscle synergies: evaluation on simulated and experimental data sets. J. Neurophysiology 95 (4), 2199–2212. doi:10.1152/jn.00222.2005
Wang, C., Zhang, S., Hu, J., Huang, Z., and Shi, C. (2021). Upper-limb muscle synergy features in human-robot interaction with circle-drawing movements. Appl. Bionics Biomechanics 2021 (1), 1–11. doi:10.1155/2021/8850785
Zajac, F. E. (1989). Muscle and tendon: properties, models, scaling, and application to biomechanics and motor control. Crit. Rev. Biomed. Eng. 17 (4), 359–411. Available online at: https://hal.science/hal-04849267.
Keywords: EMG-driven model, muscle synergy, cross-country sit-skiing, muscle activation, upper limb
Citation: Chen X, Yuan Z, Gao X, Zhang Y, Liu C and Huo B (2025) Assessment of synergy-assisted EMG-driven NMSK model for upper limb muscle activation prediction in cross-country sit-skiing double poling. Front. Bioeng. Biotechnol. 13:1585127. doi: 10.3389/fbioe.2025.1585127
Received: 28 February 2025; Accepted: 31 July 2025;
Published: 18 August 2025.
Edited by:
Wencai Liu, Shanghai Jiao Tong University, ChinaReviewed by:
Kuan Tao, Beijing Sport University, ChinaMohammad Fazle Rabbi, Bond University, Australia
Copyright © 2025 Chen, Yuan, Gao, Zhang, Liu and Huo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Chenglin Liu, bGl1Y2hlbmdsaW5AY3VwZXMuZWR1LmNu; Bo Huo, aHVvYm9AY3VwZXMuZWR1LmNu