Foot arch deformation and plantar fascia loading during running with rearfoot strike and forefoot strike: A dynamic finite element analysis

Forefoot strike becomes popular among runners because it facilitates better impact attenuation. However, forefoot strike may overload the plantar fascia and impose risk of plantar fasciitis. This study aimed to examine and compare the foot arch deformation and plantar fascia tension between different foot strike techniques in running using a computational modelling approach. A three-dimensional finite element foot model was reconstructed from the MRI of a healthy runner. The foot model included twenty bones, bulk soft tissue, ligaments, tendons, and plantar fascia. The time-series data of segmental kinematics, foot muscle force, and ankle joint reaction force were derived from a musculoskeletal model of the same participant based on the motion capture analysis and input as the boundary conditions for the finite element analysis. Rearfoot strike and forefoot strike running were simulated using a dynamic explicit solver. The results showed that, compared to rearfoot strike, forefoot strike reduced the foot arch height by 9.12% and increased the medial longitudinal arch angle by 2.06%. Forefoot strike also increased the plantar connective tissues stress by 18.28–200.11% and increased the plantar fascia tensile force by 18.71–109.10%. Although it is currently difficult to estimate the threshold value of stress or force that results in injury, forefoot strike runners appeared to be more vulnerable to plantar fasciitis.     

Foot strike pattern,step rate,and trunk posture combined gait modifications to reduce impact loading during running

Elevated impact loading can be detrimental to runners as it has been linked to the increased risk of tibial stress fracture and plantar fasciitis. The objective of this study was to investigate the combined effects of foot strike pattern, step rate, and anterior trunk lean gait modifications on impact loading in runners. Nineteen healthy runners performed 12 separate gait modification trials involving: three foot strike patterns (rearfoot, midfoot, and forefoot strike), two step rates (natural and 10% increased), and two anterior trunk lean postures (natural and 10-degree increased flexion). Overall, forefoot strike combined with increased step rate led to the lowest impact loading rates, and rearfoot strike combined with anterior trunk lean led to the highest impact loading rates. In addition, there were interaction effects between foot strike pattern and step rate on awkwardness and effort, such that it was both more natural and easier to transition to a combined gait modification involving forefoot strike and increased step rate than to an isolated gait modification involving either forefoot strike or increased step rate. These findings could help to inform gait modifications for runners to reduce impact loading and associated injury risks.     

Joint contact loading in forefoot and rearfoot strike patterns during running

Research concerning forefoot strike pattern (FFS) versus rearfoot strike pattern (RFS) running has focused on the ground reaction force even though internal joint contact forces are a more direct measure of the loads responsible for injury. The main purpose of this study was to determine the internal loading of the joints for each strike pattern. A secondary purpose was to determine if converted FFS and RFS runners can adequately represent habitual runners with regards to the internal joint loading. Using inverse dynamics to calculate the net joint moments and reaction forces and optimization techniques to estimate muscle forces, we determined the axial compressive loading at the ankle, knee, and hip. Subjects consisted of 15 habitual FFS and 15 habitual RFS competitive runners. Each subject ran at a preferred running velocity with their habitual strike pattern and then converted to the opposite strike pattern. Plantar flexor muscle forces and net ankle joint moments were greater in the FFS running compared to the RFS running during the first half of the stance phase. The average contact forces during this period increased by 41.7% at the ankle and 14.4% at the knee joint during FFS running. Peak ankle joint contact force was 1.5 body weights greater during FFS running (p<0.05). There was no evidence to support a difference between habitual and converted running for joint contact forces. The increased loading at the ankle joint for FFS is an area of concern for individuals considering altering their foot strike pattern.     

A simple field method to identify foot strike pattern during running

Identifying foot strike patterns in running is an important issue for sport clinicians, coaches and footwear industrials. Current methods allow the monitoring of either many steps in laboratory conditions or only a few steps in the field. Because measuring running biomechanics during actual practice is critical, our purpose is to validate a method aiming at identifying foot strike patterns during continuous field measurements. Based on heel and metatarsal accelerations, this method requires two uniaxial accelerometers. The time between heel and metatarsal acceleration peaks (THM) was compared to the foot strike angle in the sagittal plane (α_foot) obtained by 2D video analysis for various conditions of speed, slope, footwear, foot strike and state of fatigue. Acceleration and kinematic measurements were performed at 1000 Hz and 120 Hz, respectively, during 2-min treadmill running bouts. Significant correlations were observed between THM and α_foot for 14 out of 15 conditions. The overall correlation coefficient was r=0.916 (P<0.0001, n=288). The THM method is thus highly reliable for a wide range of speeds and slopes, and for all types of foot strike except for extreme forefoot strike during which the heel rarely or never strikes the ground, and for different footwears and states of fatigue. We proposed a classification based on THM: FFS<−5.49 ms<MFS<15.2 ms<RFS. With only a few precautions being necessary to ensure appropriate use of this method, it is reliable for distinguishing rearfoot and non-rearfoot strikers in situ.     

Foot internal stress distribution during impact in barefoot running as function of the strike pattern

The aim of the present study is to examine the impact absorption mechanism of the foot for different strike patterns (rearfoot, midfoot and forefoot) using a continuum mechanics approach. A three-dimensional finite element model of the foot was employed to estimate the stress distribution in the foot at the moment of impact during barefoot running. The effects of stress attenuating factors such as the landing angle and the surface stiffness were also analyzed. We characterized rear and forefoot plantar sole behavior in an experimental test, which allowed for refined modeling of plantar pressures for the different strike patterns. Modeling results on the internal stress distributions allow predictions of the susceptibility to injury for particular anatomical structures in the foot.     

Acute changes in foot strike pattern and cadence affect running parameters associated with tibial stress fractures

Tibial stress fractures are a common and debilitating injury that occur in distance runners. Runners may be able to decrease tibial stress fracture risk by adopting a running pattern that reduces biomechanical parameters associated with a history of tibial stress fracture. The purpose of this study was to test the hypothesis that converting to a forefoot striking pattern or increasing cadence without focusing on changing foot strike type would reduce injury risk parameters in recreational runners. Running kinematics, ground reaction forces and tibial accelerations were recorded from seventeen healthy, habitual rearfoot striking runners while running in their natural running pattern and after two acute retraining conditions: (1) converting to forefoot striking without focusing on cadence and (2) increasing cadence without focusing on foot strike. We found that converting to forefoot striking decreased two risk factors for tibial stress fracture: average and peak loading rates. Increasing cadence decreased one risk factor: peak hip adduction angle. Our results demonstrate that acute adaptation to forefoot striking reduces different injury risk parameters than acute adaptation to increased cadence and suggest that both modifications may reduce the risk of tibial stress fractures.     

Foot joint coupling variability differences between habitual rearfoot and forefoot runners prior to and following an exhaustive run

As joint coupling variability has been associated with running-related lower extremity injury, the purpose of this study was to identify how variability within the foot may be different between forefoot (FFS) and rearfoot strike (RFS) runners. Identifying typical variability in uninjured runners may contribute to understanding of ideal coordination associated with running foot strike patterns.Fifteen FFS and 15 RFS runners performed a maximal-effort 5 km treadmill run. A 7-segment foot model identified 6 functional articulations (rearfoot, medial and lateral midfoot and forefoot, and 1st metatarsophalangeal) for analysis. Beginning and end of the run motion capture data were analyzed. Vector coding was used to calculate 6 joint couples. Standard deviations of the coupling angles were used to identify variability within subphases of stance (loading, mid-stance, terminal, and pre-swing). Mixed between-within subjects ANOVAs compared differences between the foot strikes, pre and post run.Increased variability was identified within medial foot coupling for FFS and within lateral foot coupling for RFS during loading and mid-stance. The exhaustive run increased variability during mid-stance for both groups.Interpretation. Joint coupling variability profiles for FFS and RFS runners suggest different foot regions have varying coordination needs which should be considered when comparing the strike patterns.     

Foot strike patterns of runners at the 15-km point during an elite-level half marathon

There are various recommendations by many coaches regarding foot landing techniques in distance running that are meant to improve running performance and prevent injuries. Several studies have investigated the kinematic and kinetic differences between rearfoot strike (RFS), midfoot strike (MFS), and forefoot strike (FFS) patterns at foot landing and their effects on running efficiency on a treadmill and over ground conditions. However, little is known about the actual condition of the foot strike pattern during an actual road race at the elite level of competition. The purpose of the present study was to document actual foot strike patterns during a half marathon in which elite international level runners, including Olympians, compete. Four hundred fifteen runners were filmed by 2 120-Hz video cameras in the height of 0.15 m placed at the 15.0-km point and obtained sagittal foot landing and taking off images for 283 runners. Rearfoot strike was observed in 74.9% of all analyzed runners, MFS in 23.7%, and FFS in 1.4%. The percentage of MFS was higher in the faster runners group, when all runners were ranked and divided into 50 runner groups at the 15.0-km point of the competition. In the top 50, which included up to the 69th place runner in actual order who passed the 15-km point at 45 minutes, 53 second (this speed represents 5.45 m x s(-1), or 15 minutes, 17 seconds per 5 km), RFS, MFS, and FFS were 62.0, 36.0, and 2.0%, respectively. Contact time (CT) clearly increased for the slower runners, or the placement order increased (r = 0.71, p < or = 0.05). The CT for RFS + FFS for every 50 runners group significantly increased with increase of the placement order. The CT for RFS was significantly longer than MFS + FFS (200.0 +/- 21.3 vs. 183.0 +/- 16 millisecond). Apparent inversion (INV) of the foot at the foot strike was observed in 42% of all runners. The percentage of INV for MFS was higher than for RFS and FFS (62.5, 32.0, and 50%, respectively). The CT with INV for MFS + FFS was significantly shorter than the CT with and without INV for RFS. Furthermore, the CT with INV was significantly shorter than push-off time without INV for RFS. The findings of this study indicate that foot strike patterns are related to running speed. The percentage of RFS increases with the decreasing of the running speed; conversely, the percentage of MFS increases as the running speed increases. A shorter contact time and a higher frequency of inversion at the foot contact might contribute to higher running economy.     

Multi-segment foot kinematics and ground reaction forces during gait of individuals with plantar fasciitis

Background

Clinically, plantar fasciitis (PF) is believed to be a result and/or prolonged by overpronation and excessive loading, but there is little biomechanical data to support this assertion. The purpose of this study was to determine the differences between healthy individuals and those with PF in (1) rearfoot motion, (2) medial forefoot motion, (3) first metatarsal phalangeal joint (FMPJ) motion, and (4) ground reaction forces (GRF).

Methods

We recruited healthy (n=22) and chronic PF individuals (n=22, symptomatic over three months) of similar age, height, weight, and foot shape (p>0.05). Retro-reflective skin markers were fixed according to a multi-segment foot and shank model. Ground reaction forces and three dimensional kinematics of the shank, rearfoot, medial forefoot, and hallux segment were captured as individuals walked at 1.35 ms⁻¹.

Results

Despite similarities in foot anthropometrics, when compared to healthy individuals, individuals with PF exhibited significantly (p<0.05) (1) greater total rearfoot eversion, (2) greater forefoot plantar flexion at initial contact, (3) greater total sagittal plane forefoot motion, (4) greater maximum FMPJ dorsiflexion, and (5) decreased vertical GRF during propulsion.

Conclusion

These data suggest that compared to healthy individuals, individuals with PF exhibit significant differences in foot kinematics and kinetics. Consistent with the theoretical injury mechanisms of PF, we found these individuals to have greater total rearfoot eversion and peak FMPJ dorsiflexion, which may put undue loads on the plantar fascia. Meanwhile, increased medial forefoot plantar flexion at initial contact and decreased propulsive GRF are suggestive of compensatory responses, perhaps to manage pain.     

Dynamic Patterns of Forces and Loading Rate in Runners with Unilateral Plantar Fasciitis: A Cross-Sectional Study

Aim/Hypothesis

The etiology of plantar fasciitis (PF) has been related to several risk factors, but the magnitude of the plantar load is the most commonly described factor. Although PF is the third most-common injury in runners, only two studies have investigated this factor in runners, and their results are still inconclusive regarding the injury stage.

Objective

Analyze and compare the plantar loads and vertical loading rate during running of runners in the acute stage of PF to those in the chronic stage of the injury in relation to healthy runners.

Methods

Forty-five runners with unilateral PF (30 acute and 15 chronic) and 30 healthy control runners were evaluated while running at 12 km/h for 40 meters wearing standardized running shoes and Pedar-X insoles. The contact area and time, maximum force, and force-time integral over the rearfoot, midfoot, and forefoot were recorded and the loading rate (20–80% of the first vertical peak) was calculated. Groups were compared by ANOVAs (p<0.05).

Results

Maximum force and force-time integral over the rearfoot and the loading rate was higher in runners with PF (acute and chronic) compared with controls (p<0.01). Runners with PF in the acute stage showed lower loading rate and maximum force over the rearfoot compared to runners in the chronic stage (p<0.01).

Conclusion

Runners with PF showed different dynamic patterns of plantar loads during running over the rearfoot area depending on the injury stage (acute or chronic). In the acute stage of PF, runners presented lower loading rate and forces over the rearfoot, possibly due to dynamic mechanisms related to pain protection of the calcaneal area.     

Midtarsal locking,the windlass mechanism,and running strike pattern: A kinematic and kinetic assessment

Changes in running strike pattern affect ankle and knee mechanics, but little is known about the influence of strike pattern on the joints distal to the ankle. The purpose of this study was to explore the effects of forefoot strike (FFS) and rearfoot strike (RFS) running patterns on foot kinematics and kinetics, from the perspectives of the midtarsal locking theory and the windlass mechanism. Per the midtarsal locking theory, we hypothesized that the ankle would be more inverted in early stance when using a FFS, resulting in decreased midtarsal joint excursions and increased dynamic stiffness. Associated with a more engaged windlass mechanism, we hypothesized that a FFS would elicit increased metatarsophalangeal joint excursions and negative work in late stance. Eighteen healthy female runners ran overground with both FFS and RFS patterns. Instrumented motion capture and a validated multi-segment foot model were used to analyze midtarsal and metatarsophalangeal joint kinematics and kinetics. During early stance in FFS the ankle was more inverted, with concurrently decreased midtarsal eversion (p < 0.001) and abduction excursions (p = 0.003) but increased dorsiflexion excursion (p = 0.005). Dynamic midtarsal stiffness did not differ (p = 0.761). During late stance in FFS, metatarsophalangeal extension was increased (p = 0.009), with concurrently increased negative work (p < 0.001). In addition, there was simultaneously increased midtarsal positive work (p < 0.001), suggesting enhanced power transfer in FFS. Clear evidence for the presence of midtarsal locking was not observed in either strike pattern during running. However, the windlass mechanism appeared to be engaged to a greater extent during FFS.     

Quantifying rearfoot-forefoot coordination in human walking

A method is proposed to facilitate the quantification and interpretation of inter-joint/-segment coordination. This technique is illustrated using rearfoot-forefoot kinematic data. We expand existing vector coding techniques and introduce a set of operational terms through which the coordinative patterns between the rearfoot segment and the forefoot segment are summarized: in-phase, anti-phase, rearfoot phase and forefoot phase. The literature on foot mechanics has characterized the stable foot at pushoff by a decreasing medial longitudinal arch angle in the sagittal plane, which is accompanied by forefoot pronation and concurrent rearfoot supination-in other words, anti-phase motion. Nine skin markers were placed on the rearfoot and forefoot segments according to a multi-segment foot model. Three healthy subjects performed standing calibration and walking trials (1.35ms(-1)), while a three-dimensional motion capture system acquired their kinematics. Rearfoot-forefoot joint angles were derived and the arch angle was inferred from the sagittal plane. Coupling angles of rearfoot and forefoot segments were derived and categorized into one of the four coordination patterns. Arch kinematics were consistent with the literature; in stance, the arch angle reached peak dorsiflexion, and then decreased rapidly. However, anti-phase coordination was not the predominant pattern during mid- or late stance. These preliminary data suggest that the coordinative interactions between the rearfoot and the forefoot are more complicated than previously described. The technique offers a new perspective on coordination and may provide insight into deformations of underlying tissues, such as the plantar fascia.     

Effects of a foot orthosis inspired by the concept of a twisted osteoligamentous plate on the kinematics of foot-ankle complex during walking: A proof of concept

It has been suggested that the foot acts as a twisted osteoligamentous plate to control pronation and facilitate supination during walking. The aim of this study was to investigate the effect of an orthosis inspired by the concept of a foot’s twisted osteoligamentous plate on the kinematics of foot-ankle complex. Thirty-five subjects underwent a kinematic assessment of the foot-ankle complex during walking using three different orthoses: (1) Twisted Plate Spring (TPS) orthosis: inspired by the concept of a twisted osteoligamentous plate shape and made with a spring-like material (carbon fiber); (2) Flat orthosis: control orthosis made of a non-elastic material with a non-inclined surface; and (3) Rigid orthosis: control orthosis made of a non-elastic material, with the same shape of the TPS. Repeated measures analyses of variance demonstrated that the TPS reduced the duration and magnitude of rearfoot eversion (p ≤ 0.03), increased rearfoot inversion relative to shank (p < 0.01), increased forefoot eversion relative to rearfoot (p < 0.01), and increased peak of plantar flexion of forefoot relative to rearfoot during the propulsive phase (p = 0.01) compared to Flat orthosis. The effects of the TPS were different from the Rigid orthosis, demonstrating that, alongside shape, material properties were a determinant factor for the obtained results. The findings of this study help clarify the role of a mechanism similar to a twisted osteoligamentous plate on controlling foot pronation and facilitating supination during the stance phase of walking.     

Effect of the Fatigue Induced by a 110-km Ultramarathon on Tibial Impact Acceleration and Lower Leg Kinematics

Ultramarathon runners are exposed to a high number of impact shocks and to severe neuromuscular fatigue. Runners may manage mechanical stress and muscle fatigue by changing their running kinematics. Our purposes were to study (i) the effects of a 110-km mountain ultramarathon (MUM) on tibial shock acceleration and lower limb kinematics, and (ii) whether kinematic changes are modulated according to the severity of neuromuscular fatigue. Twenty-three runners participated in the study. Pre- and post-MUM, neuromuscular tests were performed to assess knee extensor (KE) and plantar flexor (PF) central and peripheral fatigue, and a treadmill running bouts was completed during which step frequency, peak acceleration, median frequency and impact frequency content were measured from tibial acceleration, as well as foot-to-treadmill, tibia-to-treadmill, and ankle flexion angles at initial contact, and ankle range of motion using video analysis. Large neuromuscular fatigue, including peripheral changes and deficits in voluntary activation, was observed in KE and PF. MVC decrements of ~35% for KE and of ~28% for PF were noted. Among biomechanical variables, step frequency increased by ~2.7% and the ankle range of motion decreased by ~4.1% post-MUM. Runners adopting a non rearfoot strike pre-MUM adopted a less plantarflexed foot strike pattern post-MUM while those adopting a rearfoot strike pre-MUM tended to adopt a less dorsiflexed foot strike pattern post-MUM. Positive correlations were observed between percent changes in peripheral PF fatigue and the ankle range of motion. Peripheral PF fatigue was also significantly correlated to both percent changes in step frequency and the ankle angle at contact. This study suggests that in a fatigued state, ultratrail runners use compensatory/protective adjustments leading to a flatter foot landing and this is done in a fatigue dose-dependent manner. This strategy may aim at minimizing the overall load applied to the musculoskeletal system, including impact shock and muscle stretch.     

Effect of peak pressure and pressure gradient on subsurface shear stresses in the neuropathic foot

The pressure distribution on the plantar surface of the foot may provide insights into the stresses within the subsurface tissues of patients with diabetes mellitus and peripheral neuropathy (PN) who are at risk for skin breakdown. The purposes of this study were to (1) estimate the stress distribution in the subsurface soft tissue from a measured surface pressure distribution and determine any differences between values in the forefoot and rearfoot, and (2) determine the relationship between maximum shear stress (MSS) (magnitude and depth) and characteristics of the pressure distribution. The measured in-shoe pressure distributions during walking characterized by the peak plantar pressure and maximum pressure gradient on the plantar surface of the feet for 20 subjects with diabetes, PN and history of a mid foot or forefoot plantar ulcer were analyzed. The effects of peak pressure and maximum pressure gradient at the peak pressure location on the stress components in the subsurface soft tissue were studied using a potential function method to estimate the subsurface tissue stress. The calculated MSSs are larger in magnitude and located closer to the surface in the forefoot, where most skin breakdown occurs, compared to the rearfoot. In addition, the MSS (magnitude and depth) is highly correlated with the pressure gradient (r=-0.77 & 0.61) and the peak pressure (r=-0.61 & 0.91). The peak pressure and the maximum pressure gradient obtained from the surface pressure distribution appear to be important variables to identify where MSSs are located in the subsurface tissues on the plantar foot that may lead to skin break down.     

Predicting timing of foot strike during running,independent of striking technique,using principal component analysis of joint angles

As 3-dimensional (3D) motion-capture for clinical gait analysis continues to evolve, new methods must be developed to improve the detection of gait cycle events based on kinematic data. Recently, the application of principal component analysis (PCA) to gait data has shown promise in detecting important biomechanical features. Therefore, the purpose of this study was to define a new foot strike detection method for a continuum of striking techniques, by applying PCA to joint angle waveforms. In accordance with Newtonian mechanics, it was hypothesized that transient features in the sagittal-plane accelerations of the lower extremity would be linked with the impulsive application of force to the foot at foot strike. Kinematic and kinetic data from treadmill running were selected for 154 subjects, from a database of gait biomechanics. Ankle, knee and hip sagittal plane angular acceleration kinematic curves were chained together to form a row input to a PCA matrix. A linear polynomial was calculated based on PCA scores, and a 10-fold cross-validation was performed to evaluate prediction accuracy against gold-standard foot strike as determined by a 10 N rise in the vertical ground reaction force. Results show 89–94% of all predicted foot strikes were within 4 frames (20 ms) of the gold standard with the largest error being 28 ms. It is concluded that this new foot strike detection is an improvement on existing methods and can be applied regardless of whether the runner exhibits a rearfoot, midfoot, or forefoot strike pattern.     

Differences in muscle activity between natural forefoot and rearfoot strikers during running

Running research has focused on reducing injuries by changing running technique. One proposed method is to change from rearfoot striking (RFS) to forefoot striking (FFS) because FFS is thought to be a more natural running pattern that may reduce loading and injury risk. Muscle activity affects loading and influences running patterns; however, the differences in muscle activity between natural FFS runners and natural RFS runners are unknown. The purpose of this study was to measure muscle activity in natural FFS runners and natural RFS runners. We tested the hypotheses that tibialis anterior activity would be significantly lower while activity of the plantarflexors would be significantly greater in FFS runners, compared to RFS runners, during late swing phase and early stance phase. Gait kinematics, ground reaction forces and electromyographic patterns of ten muscles were collected from twelve natural RFS runners and ten natural FFS runners. The root mean square (RMS) of each muscle?s activity was calculated during terminal swing phase and early stance phase. We found significantly lower RMS activity in the tibialis anterior in FFS runners during terminal swing phase, compared to RFS runners. In contrast, the medial and lateral gastrocnemius showed significantly greater RMS activity in terminal swing phase in FFS runners. No significant differences were found during early stance phase for the tibialis anterior or the plantarflexors. Recognizing the differences in muscle activity between FFS and RFS runners is an important step toward understanding how foot strike patterns may contribute to different types of injury.     

The extensibility of the plantar fascia influences the windlass mechanism during human running

The arch of the human foot is unique among hominins as it is compliant at ground contact but sufficiently stiff to enable push-off. These behaviours are partly facilitated by the ligamentous plantar fascia whose role is central to two mechanisms. The ideal windlass mechanism assumes that the plantar fascia has a nearly constant length to directly couple toe dorsiflexion with a change in arch shape. However, the plantar fascia also stretches and then shortens throughout gait as the arch-spring stores and releases elastic energy. We aimed to understand how the extensible plantar fascia could behave as an ideal windlass when it has been shown to strain throughout gait, potentially compromising the one-to-one coupling between toe arc length and arch length. We measured foot bone motion and plantar fascia elongation using high-speed X-ray during running. We discovered that toe plantarflexion delays plantar fascia stretching at foot strike, which probably modifies the distribution of the load through other arch tissues. Through a pure windlass effect in propulsion, a quasi-isometric plantar fascia''s shortening is delayed to later in stance. The plantar fascia then shortens concurrently to the windlass mechanism, likely enhancing arch recoil at push-off.     

Initial contact and toe off event identification for rearfoot and non-rearfoot strike pattern treadmill running at different speeds

The purpose of this study was to determine the validity of kinematic based initial contact (IC) and toe-off (TO) identification algorithms for rearfoot and non-rearfoot runners across a broad range of treadmill running speeds. 14 healthy active participants completed six 20–60 s treadmill running trials at 6 speeds: 2.24, 2.68, 3.13, 3.58, 4.02, and 4.48 ms⁻¹. 3D kinematic data were collected for the last 20 s of each trial. Force plates (FP) were used as the gold standard to determine ICFP and TOFP for each step. Three algorithms for finding IC, IC_Milner, IC_Alvim, IC_Alvim-mod, and one algorithm for finding toe off, TO_Fellin, were chosen for analysis. Root mean square errors (RMSE) and difference scores with 95% confidence intervals were computed for IC, TO and stance time (ST). IC_Alvim RMSE ranged from 0.175 to 0.219 s. ST_Alvim RMSE ranged from 0.168 to 0.216 s. IC_Alvim-mod RMSE ranged from 0.105 to 0.131 s. ST_Alvim-mod RMSE ranged from 0.108 to 0.129 s. IC_Milner RMSE ranged 0.012 to 0.015 s. ST_Milner RMSE ranged 0.019 to 0.024 s. IC_Milner accuracy was inversely related to speed. IC_Milner corrected with a linear regression equation reduced differences to- 0.006 ± 0.012 s with 86% of foot strikes identified within 20 ms and 58% with 10 ms. TO_Fellin RMSE ranged from 0.012 to 0.016 s. IC_Milner adjusted for speed and TO_Fellin can be used to predict IC and TO within a broad range of treadmill running speeds (2.24–4.48 ms⁻¹) and for rearfoot and non-rearfoot strikers.