Spatiotemporal and Ground-Reaction Force Characteristics as Risk Factors for Running-Related Injury: A Secondary Analysis of a Randomized Trial Including 800+ Recreational Runners

Running stiffness and spatiotemporal parameters are similar between non-runners and runners with different experience levels

Spatiotemporal parameters and leg and joint stiffness are variables that represent the fundamental dynamics of running. Therefore, these variables may effectively differentiate between less-experienced and more-experienced runners' gait, possibly addressing differing injury rates between populations. We compared stiffness and spatiotemporal parameters between runners with different experience levels, including a group with no previous running experience. Healthy physically active participants (22.1 ± 3.6y) were divided into three groups, according to experience: experienced (running >1-year, 14-48 km/week; n = 23, 9F), novice (running <1-year, 5-21 km/week; n = 15, 4F) and non-runners (no running for the past 5 years; n = 17, 7F). Three-dimensional motion capture and force plates measured gait mechanics during overground running at 3.35 m·s-1. Knee, ankle and three-dimensional leg stiffness, contact time, flight time and step length were compared between groups using independent-measures ANCOVA (covariate = sex). No biomechanical variable was significantly different between the groups (leg: p = 0.652, Hedges' g = 0.09-0.17; ankle: p = 0.439, g = 0.07-0.19; knee: p = 0.153, g = 0.13-0.29; contact time: p = 0.592, g = 0.06-0.24; flight time: p = 0.513, g = 0.03-0.40; step length: p = 0.107, g = 0.26-0.61). Stiffness and spatiotemporal parameters were not different between runners with greater than 1-year of experience when compared to runners with less than 1-year experience and non-runners. Therefore, running gait may not differentially affect injury rates between experience levels.

Effect of double- density foot orthoses on ground reaction forces and lower limb muscle activities during running in adults with and without pronated feet

BACKGROUND

The analysis of ground reaction forces and muscle activities during walking or running can help clinicians decide on the usage of foot orthoses, particularly in individuals with pronated feet. Here, we aimed to investigate the effects of double- density foot orthoses on running kinetics and lower limb muscle activities in adults with and without pronated feet.

METHODS

Forty male adults with pronated feet (PF: n = 20, age = 25.4 ± 0.3 years, body-mass-index [BMI] = 23.3 ± 1.2 kg/m2) and without pronated feet (WPF: n = 20, age = 26.4 ± 1.0 years, BMI = 24.0 ± 0.7 kg/m2) volunteered to participate in this study. The study was registered with the Iranian Registry of Clinical Trials (IRCT20220129053865N1). Ground reaction forces (Fx, Fy, Fz) and lower limb muscle activities (e.g., m. gastrocnemius) were recorded using surface electromyography (EMGs) during running at a constant speed of 3.2 m/s over an 18-m walkway with an embedded force plate. EMGs were normalized to maximum voluntary isometric contractions.

RESULTS

Test-retest reliability for running speed data was excellent for PF and WPF groups and for the entire study cohort with intraclass correlation coefficients > 0.95. The 2-way ANOVA revealed lower peak Fz (p = 0.011; d = 1.226), lower time-to-peak for Fx (p = 0.023, d = 1.068), Fy (p = 0.025, d = 1.056), and Fz (p = 0.045, d = 0.931) during running with foot orthoses in PF individuals. During the loading phase, PF and WPF exhibited lower gastrocnemius (WPF: p = 0.005, d = 1.608; PF: p = 0.001, d = 2.430 ) and vastus medialis (WPF: p < 0.001, d = 2.532; PF: p < 0.001, d = 2.503) activity when running with foot orthoses.

CONCLUSIONS

Although double- density foot orthoses resulted in some beneficial biomechanical effects such as lower muscle activation (e.g., m. vastus medialis) in individuals with PF, foot orthoses constructions need further modifications to achieve even better running mechanics to enhance performance and lower limbinjury occurrence.

TRIAL REGISTRATION

IRCT20220129053865N1 (Registration date 19/08/2024).

Preferred temporal-spatial parameters, physical maturation, and sex are related to vertical and braking forces in adolescent long-distance runners

For adults, increasing cadence reduces ground reaction forces, but a lower preferred cadence does not predispose adults to experience higher ground reaction forces. Pubertal growth and motor control changes influence running mechanics, but it is unknown if preferred cadence or step length are associated with ground reaction forces for pre-adolescent and adolescent runners. Pre-adolescent and adolescent runners underwent an overground running analysis at a self-selected speed. Mixed model multiple linear regressions investigated the associations of preferred cadence, step length, physical maturation, and sex on ground reaction forces, while accounting for running speed and leg length. Running with a lower preferred cadence or longer preferred step length was associated with larger peak braking and vertical forces (p ≤ .01), being less physically mature was associated with larger vertical impact peak force and vertical loading rate (p ≤ .01), and being a male was associated with larger loading rates (p ≤ .01). A lower preferred cadence or longer preferred step length were associated with higher braking and vertical forces and being less physically mature or a male were associated with higher loading rates. An intervention to increase cadence/decrease step length could be considered if ground reaction forces are a concern for an adolescent runner.

Effects of Foot-Strike Pattern on Neuromuscular Function During a Prolonged Graded Run

PURPOSE

To study whether, during typical-level running, non-rear-foot strikers (non-RFS) or rear-foot strikers (RFS) presented a similar or different extent of neuromuscular fatigue after a prolonged graded run.

METHODS

Sixteen experienced male trail runners (8 non-RFS and 8 RFS) performed a 2.5-hour treadmill graded running exercise. Before and after exercise, neuromuscular tests were performed to assess neuromuscular fatigue of the plantar flexors. Biomechanical gait parameters were acquired with an instrumented treadmill, and electromyographic activity of the lower-limb muscles was collected as an index of muscle activation.

RESULTS

There were no significant time × foot strike interactions for neuromuscular (all P ≥ .742), muscle activation (all P ≥ .157), or biomechanical (all P ≥ .096) variables.

CONCLUSIONS

A dominant level running foot-strike pattern did not directly affect the extent of neuromuscular fatigue during a prolonged graded run. This suggests that no ideal running foot-strike pattern exists to minimize neuromuscular fatigue during prolonged-duration races wherein cumulative uphill and downhill segments are high, such as in trail running.

Longitudinal changes in Achilles tendon and triceps surae muscle architecture during a 156-km mountain ultramarathon

This study aimed to assess the longitudinal changes in triceps surae muscle-tendon architecture during a mountain ultramarathon. Experienced trail runners [n = 55, 78% men, age: 45.2 (13.5) yr] participated in a 156-km trail run (6,000 m climbing) consisting of six 26-km laps. The resting architectural properties of triceps surae muscle-tendon were measured using ultrasound imaging for Achilles tendon cross-sectional area (AT CSA), medial gastrocnemius muscle pennation angle, thickness, length, and fiber length. Measurements were performed the day before the race (baseline), at 52 km (T1), at 104 km (T2), at 156 km (T3), and 12 h after the race (H12). Among finishers (n = 41), there was a significant biphasic change in AT CSA during the race (P = 0.001). First, a significant decrease in AT CSA occurred between baseline and T1 (P = 0.006), with a greater decrease for participants averaging speed >8 km/h (P = 0.014). Second, there was a significant increase in AT CSA especially between T2 and T3 (P = 0.006) that was correlated with a decrease in average speed (P = 0.001) and alteration of spaciotemporal running parameters (P < 0.05). Changes in muscle-tendon architecture were not significantly different between finishers (n = 41) and nonfinishers (n = 14). In 47 participants (85.5%) who completed the follow-up, AT CSA at H12 was greater compared with baseline (P = 0.010). The main finding is the significant and biphasic modification of the AT CSA during a 156-km mountain ultramarathon with an initial decrease corresponding to mechanical stress followed by a secondary increase suggesting adaptive mechanotransduction persisting after 12 h.NEW & NOTEWORTHY Achilles tendon cross-sectional area (AT CSA) demonstrated significant adaptive modifications during a 156 km mountain ultramarathon in trained athletes. Initially, a decrease in AT CSA, especially at higher running speeds, is consecutive to the biomechanical stress on the plantar flexor muscle-tendon unit (MTU). Subsequently, there is a significant increase in AT CSA persisting up to 12 h after the race, which likely corresponds to an adaptive process to limit the compressive and tensile load on the tendon.

Effect of manual manipulation on mechanical gait parameters

CONTEXT

A variety of manual manipulation techniques are utilized in clinical practice to alleviate pain and improve musculoskeletal function. Many manual practitioners analyze gait patterns and asymmetries in their assessment of the patient, and an increasing number of gait motion capture studies are taking place with recent improvements in motion capture technology. This study is the first systematic review of whether these manual modalities have been shown to produce an objectively measurable change in gait mechanics.

OBJECTIVES

This study was designed to perform a systematic review of the literature to assess the impact of manual medicine modalities on biomechanical parameters of gait.

METHODS

A master search term composed of keywords and Medical Subject Headings (MeSH) search terms from an initial scan of relevant articles was utilized to search six databases. We screened the titles and abstracts of the resulting papers for relevance and then assessed their quality with the Cochrane Risk of Bias Tool. Clinical trials that featured both a manual manipulation intervention and multiple mechanical gait parameters were included. Case reports and other studies that only measured gait speed or other subjective measures of mobility were excluded.

RESULTS

We included 20 studies in our final analysis. They utilize manipulation techniques primarily from osteopathic, chiropractic, massage, and physiotherapy backgrounds. The conditions studied primarily included problems with the back, knee, and ankle, as well as healthy patients and Parkinson's patients. Control groups were highly variable, if not absent. Most studies measured their gait parameters utilizing either multicamera motion capture systems or force platforms.

CONCLUSIONS

Twelve of 20 papers included in the final analysis demonstrated a significant effect of manipulation on gait variables, many of which included either step length, walking speed, or sagittal range of motion (ROM) in joints of the lower extremity. However, the results and study design are too heterogeneous to draw robust conclusions from these studies as a whole. While there are initial indications that certain modalities may yield a change in certain gait parameters, the quality of evidence is low and there is insufficient evidence to conclude that manual therapies induce changes in biomechanical gait parameters. Studies are heterogeneous with respect to the populations studied and the interventions performed. Comparators were variable or absent across the studies, as were the outcome variables measured. More could be learned in the future with consistent methodology around blinding and sham treatment, and if the gait parameters measured were standardized and of a more robust clinical significance.

Habitual Running Style Matters: Duty Factor, and Not Stride Frequency, Relates to Loading Magnitude

Running style is temporally defined by a duty factor and stride frequency and believed to be related to the loading experienced during ever step. However, the exact relationship between both temporal variables and loading magnitude is still unknown. We aimed to identify the relationship between a duty factor and stride frequency with external load measures, joint reaction forces and joint moments. Thirty-one healthy female recreational runners ran across a 25-m runway at a speed of 2.30 ± 0.05 m·s-1. Ground reaction forces and motion capture data were used to determine the maximal vertical ground reaction force, the vertical instantaneous loading rate, peak braking force, peak joint extension moments and peak joint reaction forces at the knee and the ankle. The habitual duty factor and stride frequency of runners did not correlate with each other. The duty factor was found to be a significant predictor of maximal vertical ground reaction force (R2 = 0.585), peak braking force (R2 = 0.153), peak knee extension moment (R2 = 0.149), ankle plantar flexion moment (R2 = 0.225) and peak joint reaction forces at the knee (R2 = 0.591) and the ankle (R2 = 0.592), but not of the vertical instantaneous loading rate. Stride frequency had no significant predictive value. In conclusion, the maximal loading and potential injury risk of female recreational runners running with high duty factors are lower compared to those of peers running with lower duty factors.

The Dogma of Running Injuries: Perceptions of Adolescent and Adult Runners

CONTEXT

Adults perceive certain factors to increase or decrease the risk of sustaining running-related injuries, but many of their perceptions are not supported by research. Little is known about the perceptions that adolescent runners hold. Investigating perceptions for adolescent runners is needed to assist in the development of future injury educational materials, as these resources may need to be tailored differently for adolescents and adults.

OBJECTIVE

To identify factors that adolescent runners perceive as risk or protective factors for running-related injuries and to compare these perceptions with those of adult runners.

DESIGN

Cross-sectional study.

SETTING

Online survey.

PATIENTS OR OTHER PARTICIPANTS

We surveyed 302 adolescent (164 females, 138 males; age = 16.0 ± 1.4 years [range, 12-19 years]) and 357 adult runners (197 women, 160 men; age = 40.7 ± 11.8 years [range, 20-77 years]).

MAIN OUTCOME MEASURE(S)

Participants completed a survey with questions about whether factors related to training habits, footwear, biomechanics, strength, stretching, or nutrition influence the risk of sustaining a running-related injury. If ≥75% of adolescents indicated that a factor increases or decreases the risk of sustaining an injury, we considered that factor to be a perceived risk or protective factor, respectively. We also performed Fisher's exact test to compare the proportion of adolescent and adult runners who responded with "increase," "decrease," "neither increase or decrease," or "I don't know" to each question.

RESULTS

Adolescent runners perceived training habits, footwear, biomechanics, strength, stretching, and nutrition to increase or decrease the risk of sustaining a running-related injury. A larger proportion of adolescents than adults perceived that more footwear cushioning and stretching decrease injury risk, but a smaller proportion perceived that overtraining increases injury risk and strength decreases injury risk.

CONCLUSIONS

Differences in perceptions exist between adolescent and adult runners, and future educational materials and research questions may need to be tailored for different running populations.

Evaluation of the optimal number of steps to obtain reliable running spatio-temporal parameters and their variability

BACKGROUND

Spatio-temporal running parameters and their variability help to determine a runner's running style. However, determining whether a change is due to the measurement or to a specific condition such as an injury is a matter of debate, as no recommendation on the number of steps required to obtain reliable assessments exists.

RESEARCH QUESTION

What is the optimal number of steps required to measure different spatio-temporal parameters and study their variability at different running speeds?

METHODS

Twenty-five runners performed three experimental sessions of three bouts of treadmill running at 8, 10 and 12 km/h separated by 24 h. We measured cadence, stride, step, contact and flight time. We calculated the duty factor and the leg stiffness index (Kleg). Mean spatio-temporal parameters and linear (coefficient of variation, standard deviation) and non-linear (Higuchi fractal index, α1 coefficient of detrended fluctuation analysis) analyses were computed for different numbers of steps. Relative reliability was determined using the intraclass coefficient correlation. The minimal number of steps which present a good reliability level was considered as the optimal number of steps for measurement. Absolute reliability was assessed by calculating minimal detectable change.

RESULTS

To assess the mean values of spatio-temporal running parameters, between 16 and 150 steps were required. We were unable to obtain an optimal number of steps for cadence, stride and step-time variabilities for all speeds. For the linear analyses, we deduced the optimal number of steps for Kleg and the contact time (around 350 steps). Non-linear analyses measurements required between 350 and 540 steps, depending on the parameter.

SIGNIFICANCE

Researchers and clinicians should optimize experimental conditions (number of steps and running speed) depending on the parameter or the variability analysis targeted. Future studies must use absolute reliability metrics to report changes in response to a specific condition with no bias due to measurement error.

Predicting overstriding with wearable IMUs during treadmill and overground running

Running injuries are prevalent, but their exact mechanisms remain unknown largely due to limited real-world biomechanical analysis. Reducing overstriding, the horizontal distance that the foot lands ahead of the body, may be relevant to reducing injury risk. Here, we leverage the geometric relationship between overstriding and lower extremity sagittal segment angles to demonstrate that wearable inertial measurement units (IMUs) can predict overstriding during treadmill and overground running in the laboratory. Ten recreational runners matched their strides to a metronome to systematically vary overstriding during constant-speed treadmill running and showed similar overstriding variation during comfortable-speed overground running. Linear mixed models were used to analyze repeated measures of overstriding and sagittal segment angles measured with motion capture and IMUs. Sagittal segment angles measured with IMUs explained 95% and 98% of the variance in overstriding during treadmill and overground running, respectively. We also found that sagittal segment angles measured with IMUs correlated with peak braking force and explained 88% and 80% of the variance during treadmill and overground running, respectively. This study highlights the potential for IMUs to provide insights into landing and loading patterns over time in real-world running environments, and motivates future research on feedback to modify form and prevent injury.

Risk Factors for Running-Related Injury in High School and Collegiate Cross-country Runners: A Systematic Review

OBJECTIVE: To summarize and describe risk factors for running-related injuries (RRIs) among high school and collegiate cross-country runners. DESIGN: Descriptive systematic review. LITERATURE SEARCH: Four databases (Scopus, SPORTDiscus, CINAHL, Cochrane) were searched from inception to August 2023. STUDY SELECTION CRITERIA: Studies assessing RRI risk factors in high school or collegiate runners using a prospective design with at least 1 season of follow-up were included. DATA SYNTHESIS: Results across each study for a given risk factor were summarized and described. The NOS and GRADE frameworks were used to evaluate quality of each study and certainty of evidence for each risk factor. RESULTS: Twenty-four studies were included. Overall, study quality and certainty of evidence were low to moderate. Females or runners with prior RRI or increased RED-S (relative energy deficiency in sport) risk factors were most at risk for RRI, as were runners with a quadriceps angle of >20° and lower step rates. Runners with weaker thigh muscle groups had increased risk of anterior knee pain. Certainty of evidence regarding training, sleep, and specialization was low, but suggests that changes in training volume, poorer sleep, and increased specialization may increase RRI risk. CONCLUSION: The strongest predictors of RRI in high school and collegiate cross-country runners were sex and RRI history, which are nonmodifiable. There was moderate certainty that increased RED-S risk factors increased RRI risk, particularly bone stress injuries. There was limited evidence that changes in training and sleep quality influenced RRI risk, but these are modifiable factors that should be studied further in this population. J Orthop Sports Phys Ther 2024;54(2):1-13. Epub 16 November 2023. doi:10.2519/jospt.2023.11550.

Running-Related Injuries Among More Than 7000 Runners in 87 Different Countries: The Garmin-RUNSAFE Running Health Study

OBJECTIVE: To describe the cumulative injury proportion after 1000 and 2000 km of running among runners from 87 countries worldwide using wearable devices. Secondly, examine if the cumulative injury proportion differed between runners from different countries. DESIGN: Cohort study with an 18-month follow-up. METHODS: Runners aged ≥18 years who were familiar with the English language, and who were using a Garmin sports watch that supported tracking of running were eligible for inclusion. The exposure was residential country; self-reported running-related injury was the primary outcome. A generalized linear model was used to estimate the cumulative injury proportion for each country and the cumulative risk difference between the countries (country with the lowest risk used as reference). Data were analyzed at 1000 and 2000 km. RESULTS: The proportions of injured runners among the 7605 included runners from 87 different countries were 57.6% [95% CI: 56.9%, 59.0%] at 1000 km and 69.8% [95% CI: 68.3%, 71.4%] at 2000 km. Runners from the Czech Republic (40.3% [95% CI: 28.7%, 51.9%]), Austria (41.1% [95% CI: 25.9%, 52.2%]), and Germany (41.9% [95% CI: 36.0%, 47.9%]) had the lowest cumulative injury proportions at 1000 km, whereas Ireland (75.4% [95% CI: 60.4%, 90.4%]), Great Britain and Northern Ireland (73.2% [95% CI: 69.3%, 77.1%]), and Finland (67.5% [95% CI: 47.2%, 87.7%]) had the highest proportions. At 2000 km, Poland (47.7% [95% CI: 36.0%, 59.4%]), Slovenia (52.2% [95% CI: 28.5%, 75.8%]), and Croatia (54.2% [95% CI: 35.6%, 72.7%]) had the lowest proportions of injured runners. The highest cumulative injury proportions were reported in Great Britain and Northern Ireland (83.6% [95% CI: 79.6%, 87.6%]) and the Netherlands (78.3% [95% CI: 70.6%, 85.9%]). CONCLUSION: More than half of the population of adult runners from 87 countries using wearable devices sustained a running-related injury during follow-up. There were considerable between-country differences in injury proportions. J Orthop Sports Phys Ther 2024;54(2):1-9. Epub 16 November 2023. doi:10.2519/jospt.2023.11959.

Gait asymmetry in spatiotemporal and kinetic variables does not increase running-related injury risk in lower limbs: a secondary analysis of a randomised trial including 800+ recreational runners

Objective

To investigate asymmetry in spatiotemporal and kinetic variables in 800+ recreational runners, identify determinants of asymmetry, investigate if asymmetry is related to greater running injury risk and compare spatiotemporal and kinetic variables between the involved and uninvolved limb at baseline in runners having sustained an injury during follow-up.

Methods

836 healthy recreational runners (38.6% women) were tested on an instrumented treadmill at their preferred running speed at baseline and followed up for 6 months. From ground reaction force recordings, spatiotemporal and kinetic variables were derived for each lower limb. The Symmetry Index was computed for each variable. Correlations and multiple regression analyses were performed to identify potential determinants of asymmetry. Cox regression analyses investigated the association between asymmetry and injury risk. Analysis of variance for repeated measures was used to compare the involved and uninvolved limbs in runners who had sustained injuries during follow-up.

Results

107 participants reported at least one running-related injury. Leg length discrepancy and fat mass were the most common determinants of asymmetry, but all correlation coefficients were negligible (0.01-0.13) and explained variance was very low (multivariable-adjusted R2<0.01-0.03). Greater asymmetry for flight time and peak breaking force was associated with lower injury risk (HR (95% CI): 0.80 (0.64 to 0.99) and 0.96 (0.93 to 0.98), respectively). No between-limb differences were observed in runners having sustained an injury.

Conclusion

Gait asymmetry was not associated with higher injury risk for investigated spatiotemporal and kinetic variables.

Trial registration number

NCT03115437.

Changes in running economy and running technique following 6 months of running with and without wearable-based real-time feedback

BACKGROUND

An increasing number of commercially available wearables provide real-time feedback on running biomechanics with the aim to reduce injury risk or improve performance.

OBJECTIVE

Investigate whether real-time feedback by wearable insoles (ARION) alters running biomechanics and improves running economy more as compared to unsupervised running training. We also explored the correlation between changes in running biomechanics and running economy.

METHODS

Forty recreational runners were randomized to an intervention and control group and performed ~6 months of in-field training with or without wearable-based real-time feedback on running technique and speed. Running economy and running biomechanics were measured in lab conditions without feedback pre and post intervention at four speeds.

RESULTS

Twenty-two individuals (13 control, 9 intervention) completed both tests. Both groups significantly reduced their energetic cost by an average of -6.1% and -7.7% for the control and intervention groups, respectively. The reduction in energy cost did not significantly differ between groups overall (-0.07 ± 0.14 J∙kg∙m-1 , -1.5%, p = 0.63). There were significant changes in spatiotemporal metrics, but their magnitude was minor and did not differ between the groups. There were no significant changes in running kinematics within or between groups. However, alterations in running biomechanics beyond typical session-to-session variation were observed during some in-field sessions for individuals that received real-time feedback.

CONCLUSION

Alterations in running biomechanics as observed during some in-field sessions for individuals receiving wearable-based real-time feedback did not result in significant differences in running economy or running biomechanics when measured in controlled lab conditions without feedback.

Preinjury Knee and Ankle Mechanics during Running Are Reduced among Collegiate Runners Who Develop Achilles Tendinopathy

INTRODUCTION

Achilles tendinopathies (AT) are common in runners, but prospective data assessing running mechanics associated with developing AT are limited. Asymmetry in running mechanics is also considered a risk factor for injury, although it is unknown if the problematic mechanics occur on the injured limb only or are present bilaterally.

PURPOSE

This study aimed to prospectively identify differences in preinjury running biomechanics in collegiate runners who did and did not develop AT and determine if between-limb asymmetries were associated with which limb developed AT.

METHODS

Running gait data were obtained preseason on healthy collegiate cross-country runners, and AT incidence was prospectively recorded each year. Spatiotemporal, ground reaction forces, and joint kinematics and kinetics were analyzed. Linear mixed-effects models assessed differences in biomechanics between those who did and did not develop AT during the subsequent year. Generalized linear mixed-effects models determined if the asymmetry direction was associated with which limb developed an AT, with odds ratios (OR) and 95% confidence intervals (95% CI) reported.

RESULTS

Data from 106 runners were analyzed and 15 developed AT. Preinjury biomechanics of runners who developed AT showed less peak knee flexion (noninjured: 45.9° (45.2°-46.6°), injured: 43.2° (41.5°-44.9°), P < 0.001), ankle dorsiflexion (noninjured: 28.7° (28.0°-30.2°), injured: 26.0° (23.8°-28.3°), P = 0.01), and knee extensor moment (noninjured: -2.18 (N·m)·kg -1 (-2.24 to -2.12 (N·m)·kg -1 ), injured: -2.00 (N·m)·kg -1 (-2.17 to -1.84 (N·m)·kg -1 ), P = 0.02). The limb demonstrating less peak knee flexion had greater odds of sustaining an AT (OR, 1.29 (1.00-1.65), P = 0.05).

CONCLUSIONS

Knee and ankle kinematics, in addition to knee kinetics, were associated with developing an AT. Monitoring these mechanics may be useful for prospectively identifying runners at risk of developing AT.

Exploring running styles in the field through cadence and duty factor modulation

According to the dual-axis model, running styles can be defined by cadence and duty factor, variables that have been associated with running performance, economy and injury risk. To guide runners in exploring different running styles, effective instructions to modulate cadence and duty factor are needed. Such instructions have been established for treadmill running, but not for overground running, during which speed can be varied. In this study, five participants completed eight field training sessions over a 4-week training period with acoustic instructions to modulate cadence, duty factor, and, in combination, running style. Instructions were provided via audio files. Running data were collected with sports watches. Participants' experiences with guided-exploration training were evaluated with the user experience questionnaire. Data analysis revealed acoustic pacing and verbal instructions to be effective in respectively modulating cadence and duty factor, albeit with co-varying effects on speed and the non-targeted variable (i.e. duty factor or cadence). Combining acoustic pacing and verbal instructions mitigated these co-varying effects considerably, allowing for running-style modulations in intended directions (particularly towards the styles with increased cadence and increased duty factor). User experience of this form of guided-exploration training was overall positive, but could be improved in terms of autonomy (dependability). In conclusion, combining acoustic pacing and verbal instructions for running-style modulation is effective in overground running.

Preseason Vertical Center of Mass Displacement During Running and Bone Mineral Density Z-Score Are Risk Factors for Bone Stress Injury Risk in Collegiate Cross-country Runners

OBJECTIVES: To (1) assess relationships between running biomechanics, bone health, and bone stress injuries (BSIs), and (2) determine which variables constitute the most parsimonious BSI risk model among collegiate cross-country runners. DESIGN: Prospective, observational cohort study. METHODS: Running gait and bone mineral density (BMD) data from healthy collegiate cross-country runners were collected at preseason over 6 seasons. A generalized estimating equation model with backward selection was used to develop the most parsimonious model for estimating BSI risk, controlling for sex, running speed, and prior BSI. The variables assessed were spatiotemporal, ground reaction force, and joint kinematics, based on previous literature. Quasi-likelihood under the independence model criterion values and R2 values were used to select the best-fitting model. RESULTS: Data from 103 runners were included in the analysis. The best-fitting model included vertical center of mass (COM) displacement and BMD z-score. Injury risk increased with greater vertical COM displacement (unit = 0.5 cm; relative risk [RR] = 1.14; 95% confidence interval [CI]: 1.01, 1.29; P = .04) and decreased with greater BMD z-score (unit = 0.5; RR = 0.83; 95% CI: 0.72, 0.95; P = .007). The model performed similarly when step rate was included instead of vertical COM displacement. CONCLUSION: Vertical COM displacement and BMD z-score contributed to the best model for estimating risk the risk of bone stress injury in cross-country runners. Step rate was also an important variable for assessing injury risk. J Orthop Sports Phys Ther 2023;53(12):1-8. Epub 20 October 2023. doi:10.2519/jospt.2023.11860.

Towards functionally individualised designed footwear recommendation for overuse injury prevention: a scoping review

Injury prevention is essential in running due to the risk of overuse injury development. Tailoring running shoes to individual needs may be a promising strategy to reduce this risk. Novel manufacturing processes allow the production of individualised running shoes that incorporate features that meet individual biomechanical and experiential needs. However, specific ways to individualise footwear to reduce injury risk are poorly understood. Therefore, this scoping review provides an overview of (1) footwear design features that have the potential for individualisation; and (2) the literature on the differential responses to footwear design features between selected groups of individuals. These purposes focus exclusively on reducing the risk of overuse injuries. We included studies in the English language on adults that analysed: (1) potential interaction effects between footwear design features and subgroups of runners or covariates (e.g., age, sex) for running-related biomechanical risk factors or injury incidences; (2) footwear comfort perception for a systematically modified footwear design feature. Most of the included articles (n = 107) analysed male runners. Female runners may be more susceptible to footwear-induced changes and overuse injury development; future research should target more heterogonous sampling. Several footwear design features (e.g., midsole characteristics, upper, outsole profile) show potential for individualisation. However, the literature addressing individualised footwear solutions and the potential to reduce biomechanical risk factors is limited. Future studies should leverage more extensive data collections considering relevant covariates and subgroups while systematically modifying isolated footwear design features to inform footwear individualisation.

Acceleration-Based Estimation of Vertical Ground Reaction Forces during Running: A Comparison of Methods across Running Speeds, Surfaces, and Foot Strike Patterns

Twenty-seven methods of estimating vertical ground reaction force first peak, loading rate, second peak, average, and/or time series from a single wearable accelerometer worn on the shank or approximate center of mass during running were compared. Force estimation errors were quantified for 74 participants across different running surfaces, speeds, and foot strike angles and biases, repeatability coefficients, and limits of agreement were modeled with linear mixed effects to quantify the accuracy, reliability, and precision. Several methods accurately and reliably estimated the first peak and loading rate, however, none could do so precisely (the limits of agreement exceeded ±65% of target values). Thus, we do not recommend first peak or loading rate estimation from accelerometers with the methods currently available. In contrast, the second peak, average, and time series could all be estimated accurately, reliably, and precisely with several different methods. Of these, we recommend the 'Pogson' methods due to their accuracy, reliability, and precision as well as their stability across surfaces, speeds, and foot strike angles.

Reference Values and Determinants of Spatiotemporal and Kinetic Variables in Recreational Runners

Background

Identifying atypical lower limb biomechanics may help prevent the occurrence or recurrence of running-related injuries. No reference values for spatiotemporal or kinetic variables in healthy recreational runners are available in the scientific literature to support clinical management.

Purpose

To (1) present speed- and sex-stratified reference values for spatiotemporal and kinetic variables in healthy adult recreational runners; (2) identify the determinants of these biomechanical variables; and (3) develop reference regression equations that can be used as a guide in a clinical context.

Study Design

Descriptive laboratory study.

Methods

This study involved 860 healthy recreational runners (age, 19-65 years [38.5% women]) tested on an instrumented treadmill at their preferred running speed in randomly allocated, standardized running shoes with either hard or soft cushioning. Twelve common spatiotemporal and kinetic variables-including contact time, flight time, duty factor, vertical oscillation, step cadence, step length, vertical impact peak (VIP), time to VIP, vertical average loading rate, vertical stiffness, peak vertical ground-reaction force (GRF), and peak braking force-were derived from GRF recordings. Reference values for each biomechanical variable were calculated using descriptive statistics and stratified by sex and running speed category (≤7, 8, 9, 10, 11, 12, 13, 14, and ≥15 km/h). Correlations and multiple regression analyses were performed to identify potential determinants independently associated with each biomechanical variable and generate reference equations.

Results

The mean running speed was 10.5 ± 1.3 km/h and 9 ± 1.1 km/h in men and women, respectively. While all potential predictors were significantly correlated with many of the 12 biomechanical variables, only running speed showed high correlations (r > 0.7). The adjusted R2 of the multiple regression equations ranged from 0.19 to 0.88.

Conclusion

This study provides reference values and equations that may guide clinicians and researchers in interpreting spatiotemporal and kinetic variables in recreational runners.

Clinical Relevance

The reference values can be used as targets for clinicians working with recreational runners in cases where there is a clinical suspicion of a causal relationship between atypical biomechanics and running-related injury.

The biomechanical effects of insoles with different cushioning on the knee joints of people with different body mass index grades

Background: Enhancing knee protection for individuals who are overweight and obese is crucial. Cushioning insoles may improve knee biomechanics and play a significant protective role. However, the impact of insoles with varying cushioning properties on knee joints in individuals with different body mass index (BMI) categories remains unknown. Our aim was to investigate the biomechanical effects of insoles with different cushioning properties on knee joints across different BMI grades. Methods: Gravity-driven impact tests were used to characterize the cushioning properties of three types of Artificial Cartilage Foam (ACF18, 28, and 38) and ethylene-vinyl acetate (EVA) insoles. Knee joint sagittal, coronal, and vertical axis angles and moments were collected from healthy-weight (BMI 18.5-23.9 kg/m2, n = 15), overweight (BMI 24.0-27.9 kg/m2, n = 16), and obese (BMI ≥28.0 kg/m2, n = 15) individuals randomly assigned four different insoles during a drop jump. The Kruskal-Wallis test and mixed model repeated measures analysis of variance were used to compare differences among cushioning and biomechanical data across various insoles, respectively. Results: ACF showed higher cushioning than EVA, and ACF38 was the highest among the three types of ACF (all p < 0.001). During the drop jump, the knee flexion angles and moments of the ACF insoles were lower than those of the EVA insoles, the knee adduction angles of the ACF18 and ACF28 insoles were lower than those of the EVA insoles, and ACF18 insoles increased the first cushion time (all p < 0.05) for all participants in whom biomechanical variables demonstrated no interactions between insoles and BMI. Regarding the BMI-dependent biomechanical variables, compared with the EVA insoles, ACF28 insoles decreased the knee flexion angle and ACF38 insoles decreased the knee adduction and rotation moment in the healthy-weight group; ACF18 insoles decreased the knee flexion angle and ACF38 insoles decreased the knee moment in the overweight group; ACF28 insoles decreased the knee flexion and adduction moment, and ACF38 insoles decreased the knee flexion angle and rotation moment in the obese group (all p < 0.05). Conclusion: Insoles with higher cushioning properties could improve knee biomechanics and provide better knee joint protection in people across different BMI ranges.

The accuracy of commercially available instrumented insoles (ARION) for measuring spatiotemporal running metrics

Spatiotemporal metrics such as step frequency have been associated with running injuries in some studies. Wearables can measure these metrics and provide real-time feedback in-field, but are often not validated. This study assessed the validity of commercially available wireless instrumented insoles (ARION) for quantifying spatiotemporal metrics during level running at different speeds (2.78-5.0 m s-1 ,) and slopes (3° and 6° up/downhill) to an instrumented treadmill. Mean raw, percentage and absolute percentage error, and limits of agreement (LoA) were calculated. Agreement was statistically quantified using four thresholds: excellent, <5%; good, <10%; acceptable, <15%; and poor, >15% error. Excellent agreement (<5% error) was achieved for stride time across all conditions, and for step frequency across all but one condition with good agreement. Contact time and swing time generally showed at least good agreement. The mean difference across all conditions was -0.95% for contact time, 0.11% for stride time, 0.6% for swing time, -0.11% for step frequency, and -0.09% when averaged across all outcomes and conditions. The accuracy at an individual level was generally good to excellent, being <10% for all but two conditions, with these conditions being <15%. Additional experiments among four runners showed that step length could also be measured with an accuracy of 1.76% across different speeds with an updated version of the insoles. These findings suggests that the ARION wearable may not only be useful for large-scale in-field studies investigating group differences, but also to quantify spatiotemporal metrics with generally good to excellent accuracy for individual runners.

The Association of Joint Power Kinetic Variables with Running Injuries: A Case-Control Study

Background

There is conflicting data on which kinetic variables are important to consider with running injuries. Furthermore, less is understood regarding differences in these variables when considering demographics such as age, sex, weight, and running speed. The primary question was what joint power kinetic variables were different between non-injured and injured runners.

Purpose

The purpose of this study was to identify if there were differences in joint power kinetic variables between non-injured runners and injured runners.

Study Design

Case-Control Study.

Methods

Kinetic data were collected on 122 runners (26 non-injured and 96 injured) over three years with a Bertec force plated treadmill and Qualisys 3D motion capture. The subjects were considered eligible if they self-identified themselves as runners or had running as a key component of their activity. The subjects ran at a comfortable, self-selected pace while two 10-second trials of recordings were used to calculate the means of peak power generated at the hips, knees, and ankles of each gait cycle. Foot strike was categorized by kinematic data. Two sample T-tests were used to compare peak power variables at the hips, knees, and ankles between non-injured and injured runners. Logistic regression analyses examined how a combination of demographics and peak power variables were associated with injuries.

Results

No peak power variable at the hip, knee, or ankle was significantly different between injured and non-injured runners (p=0.07-0.87). However, higher hip power absorbed was found to be protective against injuries (odds ratio, .16; 95% CI .025-.88) when considering demographics using a logistic regression model including sex, foot strike, BMI, speed, age, and power variables from the hip, knee, and ankle. The area under the ROC curve was .74, which is acceptable discrimination.

Conclusion

When controlling for age, sex, BMI, foot strike, and speed; higher hip power absorbed was found to be protective against injury. This could be due to the hip muscles' unique role in absorbing force during early stance phase.

Level of Evidence

3b©The Author(s).

Adolescents running in conventional running shoes have lower vertical instantaneous loading rates but greater asymmetry than running barefoot or in partial-minimal shoes

Footwear may moderate the transiently heightened asymmetry in lower limb loading associated with peak growth in adolescence during running. This repeated-measures study compared the magnitude and symmetry of peak vertical ground reaction force and instantaneous loading rates (VILRs) in adolescents during barefoot and shod running. Ten adolescents (age, 10.6 ± 1.7 years) ran at self-selected speed (1.7 ± 0.3 m/s) on an instrumented treadmill under three counter-balanced conditions; barefoot and shod with partial-minimal and conventional running shoes. All participants were within one year of their estimated peak height velocity based on sex-specific regression equations. Foot-strike patterns, peak vertical ground reaction force and VILRs were recorded during 20 seconds of steady-state running. Symmetry of ground reaction forces was assessed using the symmetry index. Repeated-measures ANOVAs were used to compare conditions (α=.05). Adolescents used a rearfoot foot-strike pattern during barefoot and shod running. Use of conventional shoes resulted in a lower VILR (P < .05, dz = 0.9), but higher VILR asymmetry (P < .05) than running barefoot (dz = 1.5) or in partial-minimal shoes (dz = 1.6). Conventional running shoes result in a lower VILR than running unshod or in partial-minimal shoes but may have the unintended consequence of increasing VILR asymmetry. The findings may have implications for performance, musculoskeletal development and injury in adolescents.

Biomechanical adaptations during exhaustive runs at 90 to 120% of peak aerobic speed

The aim of this study was to examine how running biomechanics (spatiotemporal and kinetic variables) adapt with exhaustion during treadmill runs at 90, 100, 110, and 120% of the peak aerobic speed (PS) of a maximal incremental aerobic test. Thirteen male runners performed a maximal incremental aerobic test on an instrumented treadmill to determine their PS. Biomechanical variables were evaluated at the start, mid, and end of each run until volitional exhaustion. The change of running biomechanics with fatigue was similar among the four tested speeds. Duty factor and contact and propulsion times increased with exhaustion (P ≤ 0.004; F ≥ 10.32) while flight time decreased (P = 0.02; F = 6.67) and stride frequency stayed unchanged (P = 0.97; F = 0.00). A decrease in vertical and propulsive peak forces were obtained with exhaustion (P ≤ 0.002; F ≥ 11.52). There was no change in the impact peak with exhaustion (P = 0.41; F = 1.05). For runners showing impact peaks, the number of impact peaks increased (P ≤ 0.04; [Formula: see text] ≥ 6.40) together with the vertical loading rate (P = 0.005; F = 9.61). No changes in total, external, and internal positive mechanical work was reported with exhaustion (P ≥ 0.12; F ≤ 2.32). Results suggest a tendency towards a "smoother" vertical and horizontal running pattern with exhaustion. A smoother running pattern refers to the development of protective adjustments, leading to a reduction of the load applied to the musculoskeletal system at each running step. This transition seemed continuous between the start and end of the running trials and could be adopted by the runners to decrease the muscle force level during the propulsion phase. Despite these changes with exhaustion, there were no changes in either gesture speed (no alteration of stride frequency) or positive mechanical work, advocating that runners unconsciously organize themselves to maintain a constant whole-body mechanical work output.

Time course of muscle activation, energetics and mechanics of running in minimalist and traditional cushioned shoes during level running

The study aimed to compare the ankle muscles activation, biomechanics and energetics of running in male runners during submaximal level run using minimalist (MinRS) and traditional cushioned (TrdRS) running shoes. During 45-min running in MinRS and TrdRS, the ankle muscles pre- and co-activation, biomechanics, and energetics of running of 16 male endurance runners (25.5 ± 3.5 yr) were assessed using surface electromyography (tibialis anterior and gastrocnemius lateralis), instrumented treadmill and indirect calorimetry, respectively. The net energy cost of running (C_r) was similar for both conditions (P = 0.25) with a significant increase over time (P < 0.0001). Step frequency (P < 0.001), and total mechanical work (P = 0.001) were significantly higher in MinRS than in TrdRS with no evolution over time (P = 0.28 and P = 0.85, respectively). The ankle muscles pre- and co-activation during the contact phase did not differ between the two shoe conditions (P ≥ 0.33) or over time (P ≥ 0.15). In conclusion, during 45-min running, Cr and muscle pre- and co-activation were not significantly different between MinRS and TrdRS with significantly higher step frequency and total mechanical work noted in the former than in the latter. Moreover, C_r significantly increased during the 45-min trial in both shoe conditions along with no significant change over time in muscle activation and biomechanical variables.

Differences in kinetic variables between injured and uninjured rearfoot runners: A hierarchical cluster analysis

Running is a popular form of physical activity with a high incidence of running-related injuries. However, the etiology of running-related injuries remains elusive, possibly due to the heterogeneity of movement patterns. The purpose of this study was to investigate whether different clusters existed within a large group of injured and uninjured runners based on their kinetic gait patterns. A sample of 134 injured and uninjured runners were acquired from an existing database and 12 discrete kinetic and spatiotemporal variables which are commonly associated with running injuries were extracted from the ground reaction force waveforms. A principal components analysis followed by an unsupervised hierarchical cluster analysis was performed. The results revealed two distinct clusters of runners which were not associated with injury status (OR = 1.14 [0.57, 2.30], χ²  = 0.143, p = 0.706) or sex (OR = 1.72 [0.85, 3.49], χ²  = 2.3258, p = 0.127). These results suggest that while there appeared to be evidence for two distinct clusters within a large sample of injured and uninjured runners, there is no association between the kinetic variables and running related injuries.

Comparison of different machine learning models to enhance sacral acceleration-based estimations of running stride temporal variables and peak vertical ground reaction force

Machine learning (ML) was used to predict contact (tc) and flight (tf) time, duty factor (DF) and peak vertical force (Fv,max) from IMU-based estimations. One hundred runners ran on an instrumented treadmill (9-13 km/h) while wearing a sacral-mounted IMU. Linear regression (LR), support vector regression and two-layer neural-network were trained (80 participants) using IMU-based estimations, running speed, stride frequency and body mass. Predictions (remaining 20 participants) were compared to gold standard (kinetic data collected using the force plate) by calculating the mean absolute percentage error (MAPE). MAPEs of Fv,max did not significantly differ among its estimation and predictions (P = 0.37), while prediction MAPEs for tc, tf and DF were significantly smaller than corresponding estimation MAPEs (P ≤ 0.003). There were no significant differences among prediction MAPEs obtained from the three ML models (P ≥ 0.80). Errors of the ML models were equal to or smaller than (≤32%) the smallest real difference for the four variables, while errors of the estimations were not (15-45%), indicating that ML models were sufficiently accurate to detect a clinically important difference. The simplest ML model (LR) should be used to improve the accuracy of the IMU-based estimations. These improvements may be beneficial when monitoring running-related injury risk factors in real-world settings.

Increased oxygen uptake in well-trained runners during uphill high intensity running intervals: A randomized crossover testing

The time spent above 90% of maximal oxygen uptake ( V ˙ O₂max) during high-intensity interval training (HIIT) sessions is intended to be maximized to improve V ˙ O₂max. Since uphill running serves as a promising means to increase metabolic cost, we compared even and moderately inclined running in terms of time ≥90% V ˙ O₂max and its corresponding physiological surrogates. Seventeen well-trained runners (8 females & 9 males; 25.8 ± 6.8yrs; 1.75 ± 0.08m; 63.2 ± 8.4kg; V ˙ O₂max: 63.3 ± 4.2 ml/min/kg) randomly completed both a horizontal (1% incline) and uphill (8% incline) HIIT protocol (4-times 5min, with 90s rest). Mean oxygen uptake ( V ˙ O₂mean), peak oxygen uptake ( V ˙ O₂peak), lactate, heart rate (HR), and perceived exertion (RPE) were measured. Uphill HIIT revealed higher (p ≤ 0.012; partial eta-squared (pes) ≥ 0.351) V ˙ O₂mean (uphill: 3.3 ± 0.6 vs. horizontal: 3.2 ± 0.5 L/min; standardized mean difference (SMD) = 0.15), V ˙ O₂peak (uphill: 4.0 ± 0.7 vs. horizontal: 3.8 ± 0.7 L/min; SMD = 0.19), and accumulated time ≥90% V ˙ O₂max (uphill: 9.1 ± 4.6 vs. horizontal: 6.4 ± 4.0 min; SMD = 0.62) compared to even HIIT. Lactate, HR, and RPE responses did not show mode*time rANOVA interaction effects (p ≥ 0.097; pes ≤0.14). Compared to horizontal HIIT, moderate uphill HIIT revealed higher fractions of V ˙ O₂max at comparable perceived efforts, heartrate and lactate response. Therefore, moderate uphill HiiT notably increased time spent above 90% V ˙ O₂max.

Correspondence Between Values of Vertical Loading Rate and Oxygen Consumption During Inclined Running

Purpose

The aim of this study was to provide a theoretical model to predict the vertical loading rate (VLR) at different slopes and speeds during incline running.

Methods

Twenty-nine healthy subjects running at least once a week performed in a randomized order 4-min running trials on an instrumented treadmill at various speeds (8, 10, 12, and 14 km h−1) and slopes (− 20%, − 10%, − 5%, 0%, + 5%, + 10%, + 15%, + 20%). Heart rate, gas exchanges and ground reaction forces were recorded. The VLR was then calculated as the slope of the vertical force between 20 and 80% of the duration from initial foot contact to the impact peak.

Results

There was no difference in VLR between the four different uphill conditions at given running speeds, but it was reduced by 27% at 5% slope and by 54% at 10% slope for the same metabolic demand (similar $${\dot{\text{V}}\text{O}}_{{2}}$$ V ˙ O 2 ), when compared to level running. The average VLR measured at maximal aerobic intensity during level running would be decreased by 52.7% at + 5%, by 63.0% at + 10%, and by 73.3% at + 15% slope. Moreover, VLR was dependent on the slope in downhill conditions.

Conclusion

This study highlights the possibility to use uphill running to minimize rate of mechanical load (i.e., osteoarticular load) from foot impact on the ground and as a time-efficient exercise routine (i.e., same energy expenditure than in level running in less time).

Running-style modulation: Effects of stance-time and flight-time instructions on duty factor and cadence

BACKGROUND

The duty factor (reflecting the ratio of stance to flight time) is an important variable related to running performance, economy, and injury risk. According to the dual-axis model, the duty factor and the cadence are sufficient to describe an individual's running style at a certain speed. To test this model, one should be able to modulate both variables independently. While acoustic pacing is an established method for cadence modulation, no such method is available for duty-factor modulation.

RESEARCH QUESTIONS

Can people modulate their duty factor based on verbal instructions to change either their stance or flight time without changing their cadence? And, if so, which instruction is most effective?

METHODS

Twelve participants ran on an instrumented treadmill and completed four training blocks starting with a baseline trial and ending with a performance trial in which they followed verbal instructions to both increase and decrease their stance and flight time. Acoustic pacing at their preferred cadence was present during the first part of each trial. We calculated the duty factor and cadence for paced and non-paced parts of each trial, assessed the effectiveness of the instructions aimed at changing the duty factor, and examined the effects of instructions and acoustic pacing on cadence using Bayesian statistics.

RESULTS

The duty factor changed in intended directions with verbal instructions to increase and decrease the stance and flight time (18.04 ≤ BF₁₀ ≤ 4954.42), without differences between the instructions or during and after acoustic pacing. The instructions and acoustic pacing did not result in a consistent change in cadence (0.40 ≤ BF₁₀ ≤ 2.59).

SIGNIFICANCE

Runners can change their duty factor through verbal instructions pertaining to stance or flight time, without clear concomitant effects on cadence. Running styles can thus be altered with verbal instructions to change stance or flight time for duty-factor modulation, optionally combined with acoustic pacing to prescribe cadence.

Examination of running pattern consistency across speeds

Duty factor (DF) and step frequency (SF) are key running pattern determinants. However, running patterns may change with speed if DF and SF changes are inconsistent across speeds. We examined whether the relative positioning of runners was consistent: 1) across five running speeds (10-18 km/h) for four temporal variables [DF, SF, and their subcomponents: contact (t c ) and flight (t f ) time]; and 2) across these four temporal variables at these five speeds. Three-dimensional whole-body kinematics were acquired from 52 runners, and deviations from the median for each variable (normalised to minimum-maximum values) were extracted. Across speeds for all variables, correlations on the relative positioning of individuals were high to very high for 2-4 km/h speed differences, and moderate to high for 6-8 km/h differences. Across variables for all speeds, correlations were low between DF-SF, very high between DF-t f , and low to high between DF-t c , SF-t c , and SF-t f . Hence, the consistency in running patterns decreased as speed differences increased, suggesting that running patterns be assessed using a range of speeds. Consistency in running patterns at a given speed was low between DF and SF, corroborating suggestions that using both variables can encapsulate the full running pattern spectrum.