1,692
Views
1
CrossRef citations to date
0
Altmetric
Research Article

Marker location and knee joint constraint affect the reporting of overhead squat kinematics in elite youth football players

ORCID Icon, ORCID Icon, ORCID Icon, ORCID Icon & ORCID Icon
Pages 740-757 | Received 15 May 2020, Accepted 10 Feb 2021, Published online: 05 Mar 2021

ABSTRACT

Motion capture systems are used in the analysis and interpretation of athlete movement patterns for a variety of reasons, but data integrity remains critical regardless. The extent to which marker location or constraining degrees of freedom (DOF) in the biomechanical model impacts on this integrity lacks consensus. Ten elite academy footballers performed bilateral overhead squats using a marker-based motion capture system. Kinematic data were calculated using four different marker sets with 3DOF and 6DOF configurations for the three joint rotations of the right knee. Root mean squared error differences between marker sets ranged in the sagittal plane between 1.02 and 4.19 degrees to larger values in the frontal (1.30–6.39 degrees) and transverse planes (1.33 and 7.97 degrees). The cross-correlation function of the knee kinematic time series for all eight marker-sets ranged from excellent for sagittal plane motion (>0.99) but reduced for both coronal and transverse planes (<0.9). Two-way ANOVA repeated measures calculated at peak knee flexion revealed significant differences between marker sets for frontal and transverse planes (p < 0.05). Pairwise comparisons showed significant differences between some marker sets. Marker location and constraining DOF while measuring relatively large ranges of motion in this population are important considerations for data integrity.

Introduction

The assessment of injury risk is a key strategy for professional football teams in an attempt to reduce injury occurrence in their players (McCall et al., Citation2015, Citation2017). Furthermore, improving movement quality of footballer players has been reported as important for reducing injury risk within this population (Bagwell et al., Citation2016; Smale et al., Citation2017). The ability to accurately assess movement quality can support athletic development programmes (Bergeron et al., Citation2015; Bishop, Edwards et al., Citation2016; Marques et al., Citation2017; McCall et al., Citation2017; Scibek et al., Citation2020) and improve return to play protocols when an injury occurs (Ardern et al., Citation2018; Leporace et al., Citation2013; Whittaker et al., Citation2017). A commonly evaluated movement pattern used within practice and investigative studies is the bilateral overhead squat which encompasses a large range of knee flexion (Bishop, Edwards et al., Citation2016; Donohue et al., Citation2015; Schoenfeld, Citation2010). Commonly, this assessment of movement quality relies on the experience of relevant practitioners and can be somewhat subjective, or reliant on the ability to make a judgement based on less reliable criterion (Onate et al., Citation2012). Therefore, the use of motion capture systems could enhance movement analysis and provide practitioners with more reliable quantitative data for use in subsequent athlete development programmes, rather than qualitatively ranking and rating motion as is the current practice in football.

However, despite reported sub-millimetre accuracy, marker-based motion capture systems are subject to the limitation of soft tissue artefact and muscle contraction that affect the resultant knee kinematics (Fiorentino et al., Citation2017; Schulz & Kimmel, Citation2010). Soft tissue artefact (STA), being site, participant and movement task-specific, will impact upon the selection of marker locations and the subsequent biomechanical model utilised (Benoit et al., Citation2006; Cockcroft et al., Citation2016; McFadden et al., Citation2020; Schulz & Kimmel, Citation2010). The use of bone pins or advanced imaging validation methods to optimise kinematic data (Fiorentino et al., Citation2017; Mentiplay & Clark, Citation2018; Potvin et al., Citation2017) is impractical in sports settings. While there is evidence of the impact of marker location and biomechanical model building within literature (Cockcroft et al., Citation2016; Mentiplay & Clark, Citation2018; Robinson et al., Citation2013; Slater et al., Citation2018) to our knowledge the consideration of best practice in the selection of marker location and model configuration for assessment of the overhead squat has not yet been investigated in professional footballers. Furthermore, previous research on lower limb kinematics and injury risk reduction often lack detail on both marker set and model configuration. Previously reported and validated marker locations include the Conventional Gait Model [CGM] and variations; the Plug-in-Gait model; Helen Hayes marker sets (Baker et al., Citation2017; Duffell et al., Citation2014; Schulz & Kimmel, Citation2010); Cluster-based model (Mentiplay & Clark, Citation2018) and 6DOF gait models (Buczek et al., Citation2010; Collins et al., Citation2009; Schmitz et al., Citation2016; Żuk & Pezowicz, Citation2015). However, marker locations used in clinical gait analysis are often not validated for other movement patterns (Mentiplay & Clark, Citation2018; Schulz & Kimmel, Citation2010).

Attempts to improve the fidelity of segment and segmental interaction data, particularly for the lower limbs include the use of rigid cluster marker sets placed on the thigh and shank (Buczek et al., Citation2010; Collins et al., Citation2009; Schache et al., Citation2008) and constraining the degrees of freedom within the model for anatomical joints in an attempt to reduce the influence of STA (Duprey et al., Citation2010; Gasparutto et al., Citation2015; Potvin et al., Citation2017; Richard et al., Citation2017). However, there remains some ambiguity around the influence of these decisions on data reliability which may be in part attributable to the nature of the study participants and the type of movement patterns under scrutiny. Additionally, while some marker location and DOF decisions have been reported as reliable in sagittal and frontal plane rotations (Slater et al., Citation2018), others reported agreement only in the sagittal plane, with no improvement in data reliability when rigid clusters of markers were used (Mantovani & Lamontagne, Citation2017; Schulz & Kimmel, Citation2010). In addition, constraining the knee joint to 3DOF has been previously reported not to influence the accuracy of the kinematic outcomes and was reported to have contributed to increased error (Andersen et al., Citation2010; Fiorentino et al., Citation2017; Mentiplay & Clark, Citation2018; Potvin et al., Citation2017; Richard et al., Citation2017; Wen et al., Citation2018). In contrast, marker set test retest reliability of 3DOF and 6DOF models was deemed good to excellent for all except the transverse plane motion, with slightly higher agreement of the 6DOF marker set models reported (Mentiplay & Clark, Citation2018).

In practice, the ability to accurately report normative 6DOF knee kinematics that includes joint translation is important within professional football, but a challenge for motion capture systems (Lu & O’ Connor, Citation1999; Mentiplay & Clark, Citation2018; Ojeda et al., Citation2014; Richard et al., Citation2017; Smale et al., Citation2017). This joint translation may be important to quantify in developing athletic populations who may present with underlying hypermobility influenced by growth rates and maturation stage resulting in an increased risk of injury (Ryan et al., Citation2018; Smale et al., Citation2017). These inconsistencies within the literature regarding appropriate marker locations and model constraints, particularly when measuring large-amplitude athletic movement tasks prompted this study. The purpose of this study was to compare knee joint kinematics derived using four different lower limb marker sets and two degrees of freedom models (6DOF and 3DOF), while tracking the overhead squat in an elite youth football population. The null hypothesis was that there would be no difference between the eight configurations in the determination of knee kinematics during the overhead squat.

Methods

Ethics approval was obtained from the National University of Ireland, Galway, Medical Ethics Committee and individual informed consent gained from each participant before testing. All participants were informed of the purpose of testing and advised that they could withdraw at any stage. All data was stored in accordance with GDPR guidelines, including ensuring anonymity of all participants.

Participants

Ten Academy Elite football players (n = 10) of mean age (SD): 18.5 (±1.3) years, height 1.83 (±0.04) metres and weight 79.2 (±6.2) kg, volunteered to participate in this study. All participants self-reported their right leg as their dominant kicking leg. Participants wore tight-fitting clothing to reduce extraneous marker movement and the same footwear to standardise testing protocols.

Inclusion criteria

All participants were physically able to perform the overhead squat motion without any restriction. Participants were full-time Premier League academy football players for a minimum of two years and were engaged in supervised strength training. All participants had achieved full maturation status or 100% of Peak Adult Height at the time of testing (Khamis & Roche, Citation1994; Malina et al., Citation2015).

Exclusion criteria

Participants were excluded if they reported any musculoskeletal injury or illness that hindered participation in full training or games within the preceding 3 months.

Overhead squat

All participants were instructed using the same verbal and visual demonstration of the overhead squat (). Participants were requested to stand with their feet shoulder-width apart, holding a 120-cm wooden dowel pressed overhead with extended arms at the initiation of the trail and then to asked to complete the squat to maximum depth with good trunk control (Bishop, Edwards et al., Citation2016). The squat depth was self-determined by participants as per standardised functional movement assessment protocols (Cook et al., Citation2006; Scibek et al., Citation2020). Motion timing was controlled to two seconds down, and two seconds to return to the start position (Bishop, Villiere et al., Citation2016). All participants familiarised themselves with the movement prior to data collection under the supervision of a Chartered Physiotherapist.

Figure 1. (a) Full lower limb marker sets detailing the cluster with medial and lateral knee joint markers. The SCS for the knee joint axes was aligned with (X) positive axis in mediolateral direction, positive (Y) axis anterior/posterior and positive (Z) axis vertical directions respectively. (b) Full thigh four marker cluster detailing the Anterior partial cluster (grey), Posterior partial cluster (white) makers that were used with medial and lateral knee joint markers.

Figure 1. (a) Full lower limb marker sets detailing the cluster with medial and lateral knee joint markers. The SCS for the knee joint axes was aligned with (X) positive axis in mediolateral direction, positive (Y) axis anterior/posterior and positive (Z) axis vertical directions respectively. (b) Full thigh four marker cluster detailing the Anterior partial cluster (grey), Posterior partial cluster (white) makers that were used with medial and lateral knee joint markers.

Data collection

Kinematic data were collected using eight infrared cameras (Miqus, Qualisys Medical Ltd., Sweden) operating at 100 Hz surrounding the capture space; a dedicated area adjacent to the strength and conditioning facility of the Academy. The capture volume was calibrated according to the manufacturers’ instructions with the camera system showing a maximum calibration residual of 1 mm for each camera.

Qualisys Track Manager ™ (QTM, Version 2.16 Qualisys Medical Ltd., Sweden) was used to reconstruct the three-dimensional coordinates of each of the 19 mm spherical reflective markers placed on the following sites: right and left Anterior Superior Iliac Spine (ASIS), right and left posterior superior iliac spine (PSIS) right and left medial and lateral femoral condyles, right and left medial malleoli ankle and right and left lateral malleoli ankle. Rigid clusters were placed on the right and left lateral thighs and shanks. The markers on the lateral thigh and shank consisted of a rigid plate (131 x 80 mm) consisting of four markers with the top two markers width 7 cm apart and lengthwise 9 cm (). The rigid plate was attached with adhesive tape and velcro strapping to the lateral thigh approximately midway between the greater trochanter and the lateral knee markers. Markers were placed on all participants by the same musculoskeletal physiotherapist with over 20 years of experience, under the guidance of a biomechanist with over 15 years of motion capture experience. All reflective markers were visible in all repetitions of the movement to ensure concurrent data collection. A static trial was collected in the anatomical position for subsequent segment definitions before completing the four overhead squat movement trials .

Figure 2. Bilateral overhead squat movement test to compare marker set and models.

Figure 2. Bilateral overhead squat movement test to compare marker set and models.

Data processing

Kinematic models for each participant were created using Visual 3D™ (version 6.01.16, C-motion, Germantown, MD, USA) as follows:

Segmental and joint centre definitions

Using the right-hand orthogonal rule, the segment coordinate system (SCS) axes were aligned with the (X) positive axis in the mediolateral direction, positive (Y) axis in the anterior/posterior direction and positive (Z) axis in the vertical direction ()). The hip and knee joint centres were determined using a functional joint methods approach and individualised for all participants (Schwartz & Rozumalski, Citation2005). The squat motion served as the calibration method for the determination of segmental interaction to define the mediolateral knee joint axis (Philp et al., Citation2019), and the hip joint centre. The longitudinal axis of the femur was defined as the line between the functional hip joint centre and the midpoint between the medial and lateral femoral condyle markers projected onto the functional knee axis. All eight marker set model configurations used the same functional hip centres and functional knee axes to define segment interaction according to (Schwartz & Rozumalski, Citation2005). The shank and pelvis segments were tracked in the same way throughout the overhead squat. Thus, the primary difference between the eight models was how they track the segments articulating at the knee joint. A Cardan rotational sequence of X, Y, Z (Lees et al., Citation2010) was applied to the model.

Segment Tracking

Eight different segment reconstructions were created in Visual 3D™ using the 4 marker location sets and two joint constraint methods, 3DOF and 6DOF. For the 3DOF models, DOF was constrained to the three rotations only.

  1. Anatomical (Ana3DOF, Ana6DOF): this model used the functional hip joint centre and the medial and lateral knee markers to track the thigh (Schulz & Kimmel, Citation2010; Slater et al., Citation2018).

  2. Anterior Partial Cluster Knee (Ant3DOF, Ant6DOF): this model used two anterior thigh cluster markers and the medial and lateral knee markers to track the thigh.

  3. Cluster (Clust3DOF, Clust6DOF): all four markers on the rigid thigh cluster marker was used to track the thigh (Mentiplay & Clark, Citation2018; Schulz & Kimmel, Citation2010).

  4. Posterior Partial Cluster Knee (Post3DOF, Post6DOF): this model used the two posterior thigh cluster markers and the medial and lateral knee markers to track the thigh.

For the purpose of data extraction, the start of the overhead squat was defined as the instant when the pelvis vertical downward velocity reached a threshold of 0.1 m·s−1, and the end at the instant the participant returned to the original upright stance position. The data was time normalised to 101 data points between these two events. A 4th order low-pass Butterworth filter of 7 Hz was applied to the data (Cortes et al., Citation2014).

Data analysis

The time series of all three knee joint rotations during the overhead squat was evaluated for each marker and DOF configuration. The post-processed knee joint kinematics were examined using summary statistics, and the kinematic difference analysed. The range of differences in measurement between all eight models was expected to fall between 3 to 5 degrees (McGinley et al., Citation2009; Slater et al., Citation2018).

Root mean squared error (RMSE) was used to determine the kinematic differences between all marker sets and model configurations over the whole time series. RMSE values within a five-degree kinematic threshold were considered of ‘high fidelity’ (Slater et al., Citation2018). Cross-Correlation Function (Baxter et al., Citation2016) compared the kinematic agreement between all marker set and model combinations over the time series. A correlation of less than 0.9 was deemed substantially different (Slater et al., Citation2018). RMSE and cross-correlation values were calculated using R (R CoreTeam, Citation2019).

Peak knee joint flexion was identified as the maximum knee flexion angle at the deepest point of the overhead squat motion. The corresponding knee frontal and transverse kinematics were determined and compared at this event, as peak knee flexion kinematics reflect the risk of knee injury (Hewett & Bates, Citation2017; McLean et al., Citation2005; Myer et al., Citation2015). The effect of marker sets and model constraint on the derived knee joint kinematic data at peak knee flexion were examined using two-way repeated measures ANOVAs. The assumption of sphericity was verified following the method of Girden (Citation1992). Paired t-tests with Bonferroni correction with 95% CI were carried out to identify the mean differences between the four marker sets. One-way within-subject ANOVAs with Bonferroni correction examined any significant interaction. Partial eta2 was used to determine the effect size. The assumption of normality was checked by visual inspection of the histogram, q-q plot and the box plot of the data within the groups. Z-values of skewness and kurtosis and a Shapiro-Wilks test were also performed on the data. There were no violations of these assumptions. Statistical analyses were completed using the statistical software package SPSS (version 25, SPSS Inc., USA) with a significance level of p < 0.05.

Results

The results for the mean of the time series for the right knee kinematics of all participants using the four marker sets and the 3 and 6DOF models are illustrated in . The time series for the non-sagittal plane marker set and models demonstrate variation.

Figure 3. Group mean (n = 10) of the full time-series for the overhead squat in the sagittal plane using the four marker sets and both model configurations (3DOF, 6DOF). The right knee flexion angles are negative.

Figure 3. Group mean (n = 10) of the full time-series for the overhead squat in the sagittal plane using the four marker sets and both model configurations (3DOF, 6DOF). The right knee flexion angles are negative.

Figure 4. Group mean of the full time series for the overhead squat in the frontal plane for the four marker sets and both DOF model configurations (3DOF, 6DOF). Abduction (negative) and adduction (positive) angles.

Figure 4. Group mean of the full time series for the overhead squat in the frontal plane for the four marker sets and both DOF model configurations (3DOF, 6DOF). Abduction (negative) and adduction (positive) angles.

Figure 5. Group mean of the full time series for the overhead squat in the transverse plane for all four marker sets and both model configurations (3DOF, 6DOF). Internal rotation (positive) and external rotation (negative) angles.

Figure 5. Group mean of the full time series for the overhead squat in the transverse plane for all four marker sets and both model configurations (3DOF, 6DOF). Internal rotation (positive) and external rotation (negative) angles.

Sagittal plane

All RMSE values were below the suggested 5-degree threshold for error. The lowest RMSE between models was found between the Anterior Partial Cluster-3DOF (Ant-3DOF) and the Posterior Partial Cluster-3DOF models (Post-3DOF) (1.02 SD 1.08 degrees) with a cross correlation function of r = 0.99 (SD 0.00). The largest RMSE was noted between the Posterior Partial Cluster knee-6DOF (Post-6DOF) (4.19 SD 1.41 degrees) and the Cluster-3DOF model (Clust-3DOF) with a cross correlation function of r = 0.99 (SD 0.00). The calculated kinematic data for all the configurations in this plane correlated well (0.99) (See supplementary material).

At mean peak knee flexion, each of the marker set and DOF model combinations showed similar values (). The 2-way ANOVA for repeated measurements in this plane revealed a significant main effect for marker set, (F (1.77,15.92) = 6.00, p = 0.014, ηp2 = 0.40). Pairwise comparisons revealed a significant difference only between the Anterior Partial Cluster and Cluster marker set of 2.43 degrees (95% CI [0.24, 4.62]). No significant main effects for DOF or the interaction between marker set and DOF were observed.

Figure 6. Mean Peak flexion angles in the sagittal plane for the overhead squat using all four marker sets and both model configurations (6DOF, 3DOF). Flexion (negative) and extension (positive) angles. Examination of the significant main effect for marker set in the pairwise comparisons showed only a significant difference between the Anterior Partial Cluster (AntClusterKnee) and Cluster marker sets, averaged over DOF models.

Figure 6. Mean Peak flexion angles in the sagittal plane for the overhead squat using all four marker sets and both model configurations (6DOF, 3DOF). Flexion (negative) and extension (positive) angles. Examination of the significant main effect for marker set in the pairwise comparisons showed only a significant difference between the Anterior Partial Cluster (AntClusterKnee) and Cluster marker sets, averaged over DOF models.

Frontal plane

Lowest RMSE differences for the frontal plane were revealed between the Anterior Partial Cluster knee-6DOF (Ant-6DOF) and the Posterior Partial Cluster-6DOF models (Post-6DOF) (1.30 SD 0.93 degrees) with a cross-correlation function of r = 0.92 (SD 0.09). Highest RMSE values were found between the Posterior-6DOF (Post-6DOF) and the Cluster-3DOF (6.39 SD 2.77 degrees) with a cross correlation function of r = 0.42 (SD 0.36) (See supplementary material). The majority of configurations in this plane demonstrated lower agreement and differences > 5 degrees. Frontal plane kinematics at the moment of peak knee flexion are presented in . A significant main effect was again noted for marker set, F (2.06, 18.6) = 11.37, p = 0.001, ηp2 = 0.56. Post hoc pairwise comparisons for the main effect of marker set showed that when comparing frontal plane angles at maximal knee flexion, significant mean differences were noted between the Anatomical and Anterior Partial Cluster marker sets (−2.22 degrees, 95%CI [−3.42, −1.01), Anatomical and Posterior Partial Cluster (−3.89 degrees, 95% CI [−6.31, −1.47]) and Cluster and Posterior Partial Cluster (−3.83 degrees, 95% CI [−7.38, −0.28]) respectively. The Posterior partial cluster and both marker set models revealed mean rotations in opposite directions. No significant effects for DOF or interaction between marker set and DOF were observed.

Figure 7. Mean frontal plane angles at the peak knee flexion event during the overhead squat, using four marker sets and both models (6DOF and 3DOF). Abduction (negative) and adduction (positive) angles. Examination of the significant main effect for marker set, revealed significant differences between the anatomical, anterior partial cluster and the posterior partial cluster respectively and cluster and posterior cluster marker sets. No significant effect was noted for DOF. The posterior partial cluster marker sets by DOF model revealed rotations in opposing directions.

Figure 7. Mean frontal plane angles at the peak knee flexion event during the overhead squat, using four marker sets and both models (6DOF and 3DOF). Abduction (negative) and adduction (positive) angles. Examination of the significant main effect for marker set, revealed significant differences between the anatomical, anterior partial cluster and the posterior partial cluster respectively and cluster and posterior cluster marker sets. No significant effect was noted for DOF. The posterior partial cluster marker sets by DOF model revealed rotations in opposing directions.

Transverse plane

The RMSE was smallest between the Anterior Partial Cluster-6DOF (Ant-6DOF) and the Posterior partial Cluster-3DOF (Post-3DOF) (1.33 SD 0.56 degrees), with a cross correlation function of r = 0.92 (SD 0.08). The largest difference was found between the Posterior partial Cluster-6DOF (Post-6DOF) and the Cluster 4–3DOF models (Clust-3DOF) (7.97 RMSE SD 2.66 degrees), with a cross correlation function of r = 0.65 (SD 0.31). Many of the models in this plane revealed lower agreement and differences greater than 5 degrees (See supplementary material).

A significant interaction was observed between both marker set and DOF in this plane (F (1.33,12.00) = 13.05, p = 0.002, ηp2 = 0.59) (see ). Further examination using paired t-tests and one-way repeated measures ANOVAs revealed significantly higher mean differences for the 6DOF versus 3DOF models for the Anatomical (3.59, 95%CI [0.49, 6.69] degrees) and Cluster marker sets (3.03, 95%CI [1.06, 4.99] degrees) respectively. While the Posterior Partial Cluster (6DOF) revealed significantly lower values compared with the 3DOF model (mean difference −3.67, 95%CI [−5.83, −1.51] degrees). A one-way repeated measures ANOVA comparing all marker set (6DOF) models was also significant, F (1.65,14.87) = 93.46, p < 0.001, ηp2 = 0.91. Pairwise comparisons of the 6DOF models highlighted significant differences between all marker sets except between the Anterior Partial Cluster and Posterior Partial Cluster marker sets. These 6DOF significant differences ranged from the lowest mean difference between the Anatomical and Anterior Partial Cluster of 2.01 degrees, 95% CI [0.78, 3.24], to the highest mean difference of 7.31 degrees, 95% CI [5.08, 9.54] between the Anterior Cluster and Cluster 6DOF pairs. For the 3DOF models there was also a significant main effect for marker set, F (1.40,12.61) = 5.75, p = 0.024, ηp2 = 0.39. However, a significant difference was only observed between the Anterior Partial Cluster (Ant3DOF) and Cluster (Clust3DOF) (3.85 degrees, 95%CI [1.79, 5.91) models. A large SD was evident for all marker sets in this plane ().

Figure 8. Mean transverse plane angles at the peak knee flexion event during the overhead squat for all four marker set and both models (6DOF and 3DOF). Internal rotation (positive) and external rotation (negative) angles. The P-values are indicated for the significant comparisons based on the significant interaction between marker set and DOF model.

Significant differences between the 6DOF and 3DOF models were noted between all marker sets except the Anterior Cluster (AntClusterKnee). For all 3DOF models (lighter coloured bar chart) only significant differences were noted between the AntClusterKnee and Cluster marker sets. Large SD were revealed for all models in this plane.
Figure 8. Mean transverse plane angles at the peak knee flexion event during the overhead squat for all four marker set and both models (6DOF and 3DOF). Internal rotation (positive) and external rotation (negative) angles. The P-values are indicated for the significant comparisons based on the significant interaction between marker set and DOF model.

Discussion and implications

This study explored the influence of four different marker sets using both 3DOF and 6DOF models on the derived knee kinematic data during the bilateral overhead squat in elite youth footballers. The sagittal plane results at peak knee flexion showed similar mean values below the suggested kinematic fidelity threshold of 5-degrees for motion capture error (McGinley et al., Citation2009; Slater et al., Citation2018) with low RMSE values and excellent correlation. However, significant differences were revealed between the Anterior Partial Cluster and Cluster marker sets (2.43 degrees).

In this study, the difficulty in evaluating the quality of movement during the overhead squat in the frontal and transverse plane due to marker set and model configuration is highlighted by the differences at peak knee flexion being greater than the suggested 5-degree threshold of agreement, with higher RMSE and lower cross correlation function values. The frontal plane kinematics demonstrated significant differences between the Anatomical and both the Anterior Partial Cluster and posterior partial cluster respectively, and between the Cluster and Posterior partial cluster marker sets. The findings in the transverse plane, revealed many large differences, a large SD variation and a significant interaction for marker set and DOF. Our results reflect other studies where the combination of smaller angles and difficulty tracking both these planar motions was reported (Andersen et al., Citation2010; Duffell et al., Citation2014; Richard et al., Citation2017), in particular during the quasi-squat (Clément et al., Citation2018, Citation2015). Moreover, the frontal and transverse plane results at peak knee flexion revealed physiologically implausible kinematic data specifically for the Anatomical (6DOF), Cluster, and Posterior partial-cluster marker sets, when considered using the benchmark of gold standard validated motion (Bennett et al., Citation2018; Benoit et al., Citation2006; Donohue et al., Citation2015; Philp et al., Citation2019; Potvin et al., Citation2017; Richard et al., Citation2017; Sangeux et al., Citation2017; Sauret et al., Citation2016; Smale et al., Citation2017). Based on these results, these marker sets do not appear to provide valid calculated knee kinematic data for tracking a larger range of motion task in an athletic population.

The cluster marker set findings conflict with the existing recommendations, which have suggested that the use of clusters helps to reduce the kinematic error (Buczek et al., Citation2010; Collins et al., Citation2009). The cluster marker sets using both 3DOF and 6DOF performed poorly when tracking the frontal and transverse plane motion, with high RMSE values (frontal plane; 4–5 degrees, transverse plane; RMSE > 5 degrees) and poorer cross correlation values, in agreement with (Mentiplay & Clark, Citation2018). Further, recent results mapping soft-tissue artefact distribution suggest that marker clusters do not compensate for the kinematic error, particularly in the thigh (Barré et al., Citation2015). However, the Anatomical marker set which excluded the thigh cluster markers did not compare well with other models in this study, in contrast to previous findings (Mantovani & Lamontagne, Citation2017; Schulz & Kimmel, Citation2010; Slater et al., Citation2018). Using a four marker cluster with no additional tracking markers, or markers which are defined solely by anatomical structures did not appear to be an optimal solution for assessing the overhead squat movement in an athletic population.

Reaching consensus on the optimal modelling solution during the overhead squat, particularly in the transverse plane, is difficult. Constraining the knee joint motion to 3DOF revealed significant differences between only one pair of marker sets and suggests that this may be the solution to improve data fidelity during the overhead squat supported by some previous findings (Andersen et al., Citation2010; Duprey et al., Citation2010). The lower mean differences between marker sets in 6DOF found in this study suggested that the 6DOF models may be superior. However, the 6DOF marker set comparison revealed all three marker sets were significantly different in contrast to some previous findings (Buczek et al., Citation2010; Fiorentino et al., Citation2017; Mentiplay & Clark, Citation2018; Potvin et al., Citation2017). The transverse plane DOF, pairwise comparisons between 3DOF and 6DOF models, revealed significant differences for all marker sets except the anterior partial cluster models, which is of interest given the larger thigh muscle mass in trained athletic populations (Milsom et al., Citation2015; Peek et al., Citation2018). These results highlight the importance of DOF constraints, particularly in the transverse plane for future research and corroborate the findings of (Schellenberg et al., Citation2017). The findings suggest that future investigation of marker location and biomechanical modelling during other commonly assessed athletic movement tasks, including more dynamic movement tests, is required in order to improve the fidelity and usefulness of derived data using marker-based 3D motion capture. Further, the nature of participants could be a determining factor in the selection of marker locations and joint constraints.

There are limitations to this study, including the small sample size and the reliance upon consistent marker placement for accurate joint definitions between participants. Marker set crosstalk was not controlled for in this study (Smale et al., Citation2017). The absence of direct skeletal measurement using the gold standard, despite these methods underestimating abnormal kinematics, is an additional limitation (Akbarshahi et al., Citation2010; Li et al., Citation2012). The overhead squat in this study was performed at a relatively slow velocity, which may not be representative of more dynamic tasks and which may result in different outcomes.

Conclusion

Marker location and model configuration are important considerations when assessing athletic movement tasks. The significant differences between marker sets and model combinations during the overhead squat at peak knee flexion, particularly in the frontal and transverse planes highlight these inconsistencies. Therefore, caution is needed when comparing studies using different marker sets, models, motion tasks or populations. While it was out with the scope of this study to determine an optimal configuration for this population, primarily due to the inability to compare with reliable alternative tools, the Anterior Partial-Cluster using either the 3DOF or 6DOF models, appeared to be the most robust method for analysing knee rotational kinematics during the overhead squat in an athletic population. The benefit of not constraining DOF may also be useful in practice and research where joint translations are of interest. Further, these results suggest further research is necessary concerning marker set and model when assessing more dynamic athletic movement tasks, particularly for research into the mitigation of lower limb injury risk.

Supplemental material

Supplemental Material

Download PDF (397.9 KB)

Acknowledgments

The authors of this manuscript would like to thank Arsenal Academy Sports Science and Sports Medicine Department and Dr. Alan McCall for assistance with this paper.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Supplementary material

Supplemental data for this article can be accessed here.

Additional information

Funding

This work was supported by Science Foundation Ireland under grant number SFI/12/RC/2289 (Insight).

References

  • Akbarshahi, M., Schache, A. G., Fernandez, J. W., Baker, R., Banks, S., & Pandy, M. G. (2010). Non-invasive assessment of soft-tissue artifact and its effect on knee joint kinematics during functional activity. Journal of Biomechanics, 43(7), 1292–1301. https://doi.org/10.1016/j.jbiomech.2010.01.002
  • Andersen, M. S., Benoit, D. L., Damsgaard, M., Ramsey, D. K., & Rasmussen, J. (2010). Do kinematic models reduce the effects of soft tissue artefacts in skin marker-based motion analysis? An in vivo study of knee kinematics. Journal of Biomechanics, 43(2), 268–273. https://doi.org/10.1016/j.jbiomech.2009.08.034
  • Ardern, C. L., Ekås, G., Grindem, H., Moksnes, H., Anderson, A., Chotel, F., Cohen, M., Forssblad, M., Ganley, T. J., Feller, J. A., Karlsson, J., Kocher, M. S., LaPrade, R. F., McNamee, M., Mandelbaum, B., Micheli, L., Mohtadi, N., Reider, B., Roe, J., Seil, R., & Engebretsen, L. (2018). 2018 International Olympic Committee consensus statement on prevention, diagnosis and management of paediatric anterior cruciate ligament (ACL) injuries. Knee Surgery, Sports Traumatology, Arthroscopy: Official Journal of the ESSKA, 26(4), 989–1010. https://doi.org/10.1007/s00167-018-4865-y
  • Bagwell, J. J., Snibbe, J., Gerhardt, M., & Powers, C. M. (2016). Hip kinematics and kinetics in persons with and without cam femoroacetabular impingement during a deep squat task. Clinical Biomechanics, 31, 87–92. https://doi.org/10.1016/j.clinbiomech.2015.09.016
  • Baker, R., Leboeuf, F., Reay, J., & Sangeux, M. (2017). The conventional gait model - Success and limitations. In B. Müller, S. I. Wolf, G.-P. Brueggemann, Z. Deng, A. McIntosh, F. Miller, & W. S. Selbie (Eds.), Handbook of human motion (pp. 1–19). Springer International Publishing.
  • Barre, A., Jolles, B. M., Theumann, N., & Aminian, K. (2015). Soft tissue artifact distribution on lower limbs during treadmill gait: Influence of skin markers’ location on cluster design. Journal of Biomechanics, 48(10), 1965–1971. https://doi.org/10.1016/j.jbiomech.2015.04.007
  • Baxter, J. R., Sturnick, D. R., Demetracopoulos, C. A., Ellis, S. J., & Deland, J. T. (2016). Cadaveric gait simulation reproduces foot and ankle kinematics from population‐specific inputs. Journal of Orthopaedic Research, 34(9), 1663–1668. https://doi.org/10.1002/jor.23169
  • Bennett, H. J., Fleenor, K., & Weinhandl, J. T. (2018). A normative database of hip and knee joint biomechanics during dynamic tasks using anatomical regression prediction methods. Journal of Biomechanics, 81, 122–131. https://doi.org/10.1016/j.jbiomech.2018.10.003
  • Benoit, D. L., Ramsey, D. K., Lamontagne, M., Xu, L., Wretenberg, P., & Renström, P. (2006). Effect of skin movement artifact on knee kinematics during gait and cutting motions measured in vivo. Gait & Posture, 24(2), 152–164. https://doi.org/10.1016/j.gaitpost.2005.04.012
  • Bergeron, M. F., Mountjoy, M., Armstrong, N., Chia, M., Côté, J., Emery, C. A., Faigenbaum, A., Hall, G., Kriemler, S., Léglise, M., & Malina, R. M. (2015). International Olympic Committee consensus statement on youth athletic development. British Journal of Sports Medicine, 49(13), 843–851. http://dx.doi.org/10.1136/bjsports-2015-094962
  • Bishop, C., Edwards, M., & Turner, A. N. (2016). Screening movement dysfunctions using the overhead squat. Professional Strength & Conditioning, 42, 22–30.
  • Bishop, C., Villiere, A., & Turner, A. N. (2016). Addressing movement patterns by using the overhead squat. Professional Strength & Conditioning Journal, 40, 7–12.
  • Buczek, F. L., Rainbow, M. J., Cooney, K. M., Walker, M. R., & Sanders, J. O. (2010). Implications of using hierarchical and six degree-of-freedom models for normal gait analyses. Gait & Posture, 31(1), 57–63. https://doi.org/10.1016/j.gaitpost.2009.08.245
  • Clément, J., Dumas, R., Hagemeister, N., & De Guise, J. A. (2015). Soft tissue artifact compensation in knee kinematics by multi-body optimization: Performance of subject-specific knee joint models. Journal of Biomechanics, 48(14), 3796–3802. https://doi.org/10.1016/j.jbiomech.2015.09.040
  • Clément, J., De Guise, J. A., Fuentes, A., & Hagemeister, N. (2018). Comparison of soft tissue artifact and its effects on knee kinematics between non-obese and obese subjects performing a squatting activity recorded using an exoskeleton. Gait & Posture, 61, 197–203. https://doi.org/10.1016/j.gaitpost.2018.01.009
  • Cockcroft, J., Louw, Q., & Baker, R. (2016). Proximal placement of lateral thigh skin markers reduces soft tissue artefact during normal gait using the Conventional Gait Model. Computer Methods Biomechical Biomedical Engineering, 19(14), 1497–1504. https://doi.org/10.1080/10255842.2016.1157865
  • Collins, T. D., Ghoussayni, S. N., Ewins, D. J., & Kent, J. A. (2009). A six degrees-of-freedom marker set for gait analysis: Repeatability and comparison with a modified Helen Hayes set. Gait & Posture, 30(2), 173–180. https://doi.org/10.1016/j.gaitpost.2009.04.004
  • Cook, G., Burton, L., & Hoogenboom, B. (2006). Pre-participation screening: The use of fundamental movements as an assessment of function-part 1. North American Journal of Sports Physical Therapy, 1(2), 62–72.
  • Cortes, N., Greska, E., Ambegaonkar, J. P., Kollock, R. O., Caswell, S. V., & Onate, J. A. (2014). Knee kinematics is altered post-fatigue while performing a crossover task. Knee Surgery, Sports Traumatology, Arthroscopy, 22(9), 2202–2208. https://doi.org/10.1007/s00167-013-2673-y
  • Donohue, M. R., Ellis, S. M., Heinbaugh, E. M., Stephenson, M. L., Zhu, Q., & Dai, B. (2015). Differences and correlations in knee and hip mechanics during single-leg landing, single-leg squat, double-leg landing, and double-leg squat tasks. Research in Sports Medicine, 23(4), 394–411. https://doi.org/10.1080/15438627.2015.1076413
  • Duffell, L. D., Hope, N., & McGregor, A. H. (2014). Comparison of kinematic and kinetic parameters calculated using a cluster-based model and Vicon’s plug-in gait. Proceedings of the Institution of Mechanical Engineers. Part H, Journal of Engineering in Medicine, 228(2), 206–210. https://doi.org/10.1177/0954411913518747
  • Duprey, S., Cheze, L., & Dumas, R. (2010). Influence of joint constraints on lower limb kinematics estimation from skin markers using global optimization. Journal of Biomechanics, 43(14), 2858–2862. https://doi.org/10.1016/j.jbiomech.2010.06.010
  • Fiorentino, N. M., Atkins, P. R., Kutschke, M. J., Goebel, J. M., Foreman, K. B., & Anderson, A. E. (2017). Soft tissue artifact causes significant errors in the calculation of joint angles and range of motion at the hip. Gait & Posture, 55, 184–190. https://doi.org/10.1016/j.gaitpost.2017.03.033
  • Gasparutto, X., Sancisi, N., Jacquelin, E., Parenti-Castelli, V., & Dumas, R. (2015). Validation of a multi-body optimization with knee kinematic models including ligament constraints. Journal of Biomechanics, 48(6), 1141–1146. https://doi.org/10.1016/j.jbiomech.2015.01.010
  • Girden, E. R. (1992). ANOVA: Repeated measures. Newbury Park, CA: Sage Publications.
  • Hewett, T. E., & Bates, N. A. (2017). Preventive biomechanics: A paradigm shift with a translational approach to injury prevention. The American Journal of Sports Medicine, 45(11), 2654–2664. https://doi.org/10.1177/0363546516686080
  • Khamis, H. J., & Roche, A. F. (1994). Predicting adult stature without using skeletal age: The Khamis-Roche method. Pediatrics, 94(4), 504–507.
  • Lees, A., Asai, T., Andersen, H., Nunome, H., & Stetrzing, T. (2010). The biomechanics of kicking in soccer: A review. Journal of Sports Sciences, 28(8), 805–817. https://doi.org/10.1080/02640414.2010.481305
  • Leporace, G., Praxedes, J., Pereira, G. R., Pinto, S. M., Chagas, D., Metsavaht, L., Chame, F., & Batista, L. A. (2013). Influence of a preventive training program on lower limb kinematics and vertical jump height of male volleyball athletes. Physical Therapy in Sport, 14(1), 35–43. https://doi.org/10.1016/j.ptsp.2012.02.005
  • Li, K., Zheng, L., Tashman, S., & Zhang, X. (2012). The inaccuracy of surface-measured model-derived tibiofemoral kinematics. Journal of Biomechanics, 45(15), 2719–2723. https://doi.org/10.1016/j.jbiomech.2012.08.007
  • Lu, T. W., & O’ Connor, J. J. (1999). Bone position estimation from skin marker co-ordinates using global optimisation with joint constraints. Journal of Biomechanics, 32(2), 129–134. https://doi.org/10.1016/S0021-9290(98)00158-4
  • Malina, R. M., Rogol, A. D., Cumming, S. P., Coelho E Silva, M. J., & Figueiredo, A. J. (2015). Biological maturation of youth athletes: Assessment and implications. British Journal of Sports Medicine, 49(13), 852–859. https://doi.org/10.1136/bjsports-2015-094623
  • Mantovani, G., & Lamontagne, M. (2017). How different marker sets affect joint angles in inverse kinematics framework. Journal of Biomechanical Engineering, 139(4), 044503. https://doi.org/10.1115/1.4034708
  • Marques, V. B., Medeiros, T. M., De Souza Stigger, F., Nakamura, F. Y., & Baroni, B. M. (2017). The Functional Movement Screen (FMS™) in elite young soccer players between 14 and 20 years: Composite scores, individual-test scores and asymmetries. International Journal of Sports Physical Therapy, 12(6), 977–985. https://doi.org/10.16603/ijspt20170977
  • McCall, A., Carling, C., Davison, M., Nedelec, M., Le Gall, F., Berthoin, S., & Dupont, G. (2015). Injury risk factors, screening tests and preventative strategies: A systematic review of the evidence that underpins the perceptions and practices of 44 football (soccer) teams from various premier leagues. British Journal of Sports Medicine, 49(9), 583–589. https://doi.org/10.1136/bjsports-2014-094104
  • McCall, A., Lewin, C., O’Driscoll, G., Witvrouw, E., & Ardern, C. (2017). Return to play: The challenge of balancing research and practice. British Journal of Sports Medicine, 51(9), 702–703. https://doi.org/10.1136/bjsports-2016-096752
  • McFadden, C., Daniels, K., & Strike, S. (2020). The sensitivity of joint kinematics and kinetics to marker placement during a change of direction task. Journal of Biomechanics, 101. https://doi.org/10.1016/j.jbiomech.2020.109635
  • McGinley, J. L., Baker, R., Wolfe, R., & Morris, M. E. (2009). The reliability of three-dimensional kinematic gait measurements: A systematic review. Gait & Posture, 29(3), 360–369. https://doi.org/10.1016/j.gaitpost.2008.09.003
  • McLean, S. G., Walker, K., Ford, K., Myer, G., Hewett, T., & Van den Bogert, A. J. (2005). Evaluation of a two dimensional analysis method as a screening and evaluation tool for anterior cruciate ligament injury. British Journal of Sports Medicine, 39(6), 355–362. https://doi.org/10.1136/bjsm.2005.018598
  • Mentiplay, B. F., & Clark, R. A. (2018). Modified conventional gait model versus cluster tracking: Test-retest reliability, agreement and impact of inverse kinematics with joint constraints on kinematic and kinetic data. Gait & Posture, 64, 75–83. https://doi.org/10.1016/j.gaitpost.2018.05.033
  • Milsom, J., Naughton, R., O’Boyle, A., Iqbal, Z., Morgans, R., Drust, B., & Morton, J. P. (2015). Body composition assessment of English Premier League soccer players: A comparative DXA analysis of first team, U21 and U18 squads. Journal of Sports Sciences, 33(17), 1799–1806. https://doi.org/10.1080/02640414.2015.1012101
  • Myer, G. D., Ford, K. R., Di Stasi, S. L., Foss, K. D. B., Micheli, L. J., & Hewett, T. E. (2015). High knee abduction moments are common risk factors for patellofemoral pain (PFP) and anterior cruciate ligament (ACL) injury in girls: Is PFP itself a predictor for subsequent ACL injury? British Journal of Sports Medicine, 49(2), 118–122. https://doi.org/10.1136/bjsports-2013-092536
  • Ojeda, J., Martínez-Reina, J., & Mayo, J. (2014). A method to evaluate human skeletal models using marker residuals and global optimization. Mechanism and Machine Theory, 73, 259–272. https://doi.org/10.1016/j.mechmachtheory.2013.11.003
  • Onate, J. A., Dewey, T., Kollock, R. O., Thomas, K. S., Van Lunen, B. L., DeMaio, M., & Ringleb, S. I. (2012). Real-time intersession and interrater reliability of the functional movement screen. The Journal of Strength & Conditioning Research, 26(2), 408–415. https://doi.org/10.1519/JSC.0b013e318220e6fa
  • Peek, K., Gatherer, D., Bennett, K. J. M., Fransen, J., & Watsford, M. (2018). Muscle strength characteristics of the hamstrings and quadriceps in players from a high-level youth football (soccer) Academy. Research in Sports Medicine, 26(3), 276–288. https://doi.org/10.1080/15438627.2018.1447475
  • Philp, F., Leboeuf, F., Pandyan, A., & Stewart, C. (2019). “Dynamic knee valgus” - Are we measuring what we think we’re measuring? An evaluation of static and functional knee calibration methods for application in gait and clinical screening tests of the overhead squat and hurdle step. Gait & Posture, 70, 298–304. https://doi.org/10.1016/j.gaitpost.2019.03.006
  • Potvin, B. M., Shourijeh, M. S., Smale, K. B., & Benoit, D. L. (2017). A practical solution to reduce soft tissue artifact error at the knee using adaptive kinematic constraints. Journal of Biomechanics, 62, 124–131. https://doi.org/10.1016/j.jbiomech.2017.02.006
  • R CoreTeam. (2019). R: A language and environment for statistical computing. R Foundation for Statistical Computing. https://www.R-project.org/
  • Richard, V., Cappozzo, A., & Dumas, R. (2017). Comparative assessment of knee joint models used in multi-body kinematics optimisation for soft tissue artefact compensation. Journal of Biomechanics, 62, 95–101. https://doi.org/10.1016/j.jbiomech.2017.01.030
  • Robinson, M. A., Donnelly, C. J., Tsao, J., & Vanrenterghem, J. (2013). Impact of knee modeling approach on indicators and classification of ACL injury risk. Medicine and Science in Sports and Exercise, 46(7), 1269–1276. https://doi.org/10.1249/MSS.0000000000000236
  • Ryan, D., Lewin, C., Forsythe, S., & McCall, A. (2018). Developing world-class soccer players: An example of the academy physical development program from an English premier league team. Strength & Conditioning Journal, 40(3), 2–11. https://doi.org/10.1519/ssc.0000000000000340
  • Sangeux, M., Barré, A., & Aminian, K. (2017). Evaluation of knee functional calibration with and without the effect of soft tissue artefact. Journal of Biomechanics, 62, 53–59. https://doi.org/10.1016/j.jbiomech.2016.10.049
  • Sauret, C., Pillet, H., Skalli, W., & Sangeux, M. (2016). On the use of knee functional calibration to determine the medio-lateral axis of the femur in gait analysis: Comparison with EOS biplanar radiographs as reference. Gait & Posture, 50, 180–184. https://doi.org/10.1016/j.gaitpost.2016.09.008
  • Schache, A. G., Baker, R., & Lamoreux, L. W. (2008). Influence of thigh cluster configuration on the estimation of hip axial rotation. Gait & Posture, 27(1), 60–69. https://doi.org/10.1016/j.gaitpost.2007.01.002
  • Schellenberg, F., Taylor, W. R., Jonkers, I., & Lorenzetti, S. (2017). Robustness of kinematic weighting and scaling concepts for musculoskeletal simulation. Computer Methods in Biomechanics and Biomedical Engineering, 20(7), 720–729. https://doi.org/10.1080/10255842.2017.1295305
  • Schmitz, A., Buczek, F. L., Bruening, D., Rainbow, M. J., Cooney, K., & Thelen, D. (2016). Comparison of hierarchical and six degrees-of-freedom marker sets in analyzing gait kinematics. Computer Methods in Biomechanics and Biomedical Engineering, 19(2), 199–207. https://doi.org/10.1080/10255842.2015.1006208
  • Schoenfeld, B. J. (2010). Squatting kinematics and kinetics and their application to exercise performance. The Journal of Strength & Conditioning Research, 24(12), 3497–3506. https://doi.org/10.1519/JSC.0b013e3181bac2d7
  • Schulz, B. W., & Kimmel, W. L. (2010). Can hip and knee kinematics be improved by eliminating thigh markers? Clinical Biomechanics, 25(7), 687–692. https://doi.org/10.1016/j.clinbiomech.2010.04.002
  • Schwartz, M. H., & Rozumalski, A. (2005). A new method for estimating joint parameters from motion data. Journal of Biomechanics, 38(1), 107–116. https://doi.org/10.1016/j.jbiomech.2004.03.009
  • Scibek, E. P., Moran, M. F., & Edmond, S. L. (2020). Accuracy of functional movement screen deep squat scoring and the influence of optimized scoring criteria: A 3-dimensional kinematic approach. Journal of Sport Rehabilitation, 1(aop). https://doi.org/10.1123/jsr.2019-0041
  • Slater, A. A., Hullfish, T. J., & Baxter, J. R. (2018). The impact of thigh and shank marker quantity on lower extremity kinematics using A constrained model. BMC Musculoskeletal Disorders, 19(399). https://doi.org/10.1186/s12891-018-2329-7
  • Smale, K. B., Potvin, B. M., Shourijeh, M. S., & Benoit, D. L. (2017). Knee joint kinematics and kinetics during the hop and cut after soft tissue artifact suppression: Time to reconsider ACL injury mechanisms? Journal of Biomechanics, 62, 132–139. https://doi.org/10.1016/j.jbiomech.2017.06.049
  • Wen, Y., Huang, H., Yu, Y., Zhang, S., Yang, J., Ao, Y., & Xia, S. (2018). Effect of tibia marker placement on knee joint kinematic analysis. Gait & Posture, 60, 99–103. https://doi.org/10.1016/j.gaitpost.2017.11.020
  • Whittaker, J. L., Booysen, N., De La Motte, S., Dennett, L., Lewis, C. L., Wilson, D., McKay, C., Warner, M., Padua, D., Emery, C. A., & Stokes, M. (2017). Predicting sport and occupational lower extremity injury risk through movement quality screening: A systematic review. British Journal of Sports Medicine, 51(7), 580–585. https://doi.org/10.1136/bjsports-2016-096760
  • Zuk, M., & Pezowicz, C. (2015). Kinematic analysis of a six-degrees-of-freedom model based on ISB recommendation: A repeatability analysis and comparison with conventional gait model. Applied Bionics and Biomechanics, 2015. https://doi.org/10.1155/2015/503713