Two commonly proposed mechanical explanations for the walk-to-run transition (WRT) include the prevention of muscular over-exertion (effort) and the minimization of peak musculoskeletal loads and thus injury risk. The purpose of this study was to address these hypotheses at a joint level by analysing the effect of speed on discrete lower-limb joint kinetic parameters in humans across a wide range of walking and running speeds including walking above and running below the WRT speed. Joint work, peak instantaneous joint power, and peak joint moments in the sagittal and frontal plane of the ankle, knee and hip from eight participants were collected for 10 walking speeds (30–120% of their WRT) and 10 running speeds (80–170% of their WRT) on a force plate instrumented treadmill. Of the parameters analysed, three satisfied our statistical criteria of the ‘effort-load’ hypothesis of the WRT. Mechanical parameters that provide an acute signal (peak moment and peak power) were more strongly associated with the gait transition than parameters that reflect the mechanical function across a portion of the stride. We found that both the ankle (peak instantaneous joint power during swing) and hip mechanics (peak instantaneous joint power and peak joint moments in stance) can influence the transition from walking to running in human locomotion and may represent a cascade of mechanical events beginning at the ankle and leading to an unfavourable compensation at the hip. Both the ankle and hip mechanisms may contribute to gait transition by lowering the muscular effort of running compared with walking at the WRT speed. Although few of the examined joint variables satisfied our hypothesis of the WRT, most showed a general marked increase when switching from walking to running across all speeds where both walking and running are possible, highlighting the fundamental differences in the mechanics of walking and running. While not eliciting the WRT per se, these variables may initiate the transition between stable walking and running patterns. Those variables that were invariant of gait were predominantly found in the swing phase.

Why do we switch from walking to running to move faster? This seemingly simple question has received considerable attention, yet, whilst several factors have been proposed as triggers for the walk-to-run transition (WRT), their mechanisms remain debated. For example, it has been argued that humans (Margaria et al., 1963; Mercier et al., 1994) and other animals (Hoyt and Taylor, 1981; Rubenson et al., 2004; Watson et al., 2011) switch between walking and running to minimize metabolic energy use. However, some studies in humans have reported that the WRT occurs at a speed when it remains metabolically advantageous to walk, thus putting into question whether energy cost is the underlying determinant of the WRT (Hreljac, 1993a; Tseh et al., 2002). For example, dynamic systems theory acknowledges the importance of minimizing energy expenditure but views locomotor economy as a consequence of dynamic stability and not the cause of the transition per se. From a dynamical systems perspective, the WRT is influenced by non-specific control parameters that move the system from one stable coordination pattern to another, but does not prescribe these states. As such, the WRT may represent an abrupt shift from walking to running to avoid dynamic instability, initiated by one or more discrete mechanical variables, for example the increase in peak ankle power at the onset of running (Farris and Sawicki, 2012a).

Alternatively, discrete mechanical variables may themselves be directly associated with gait transition, rather than only reflecting the system stability or energetics. Two such commonly proposed explanations for the WRT based on discrete mechanical variables include (1) the prevention of muscular over-exertion (effort) (Hreljac, 1995; Hreljac et al., 2001; Bartlett and Kram, 2008) and (2) the minimization of peak musculoskeletal loads and therefore the chance of injury (Biewener and Taylor, 1986; Farley and Taylor, 1991; Hreljac et al., 2008). This can occur if variables associated with muscular effort and/or load become elevated when walking above the WRT but are lowered upon switching to running (referred to here as the ‘effort–load’ hypothesis of the WRT). In humans, the clearest example of these factors influencing gait transition have been identified at the ankle (Farris and Sawicki, 2012b; Hreljac, 1995a; Hreljac et al., 2008; Prilutsky and Gregor, 2001; MacLeod et al., 2014). However, the prevalence of other joints [e.g. the hip (Minetti et al., 1994)] and associated muscle groups affecting the WRT in humans may be obscured because of a lack of comprehensive information on the simultaneous individual mechanics of all the major lower leg joints across the gait cycle at speeds where both walking and running are possible. Furthermore, which discrete joint mechanical parameters are most closely associated with gait transition is poorly understood, and, in particular, whether joint loading, work or power have varying influence on the WRT remains untested.

Although these forms of mechanical analyses at the level of the joints do not assess individual muscles directly, they do provide a means to comprehensively assess the end-effect of the lower-limb muscle functions. They also provide a powerful approach by which discrete muscular effort and mechanical loading hypotheses of gait transition can be tested broadly across different joints and muscle groups, and afford further insight into the underlying mechanisms of gait transition. The purpose of this study was, therefore, to assess how discrete joint mechanics impacts the WRT by systematically measuring the effect of speed on lower-limb joint kinetics and mechanical energetics in humans across a wide range of walking and running speeds including walking above and running below the preferred WRT speed. We asked two main questions. (1) How are mechanical work, peak power and/or peak loading (moments) at the individual joints associated with the WRT? (2) Do the individual joints (ankle, knee and hip) affect the WRT differently? Furthermore, in order to better understand the mechanics of switching between walking and running, we also asked which joint work, power and moment parameters exhibit abrupt changes between walking and running but do not influence the WRT per se. Although not the central question of the present study, these parameters can inform motor control theories, such as dynamic systems theory, as they may represent the behavioural manifestation of higher order control.

In testing these questions we established statistic-based criteria for accepting a variable as satisfying the effort–load hypothesis for the WRT. We analysed variables that can provide a physiological signal for altering gait mechanics including: joint work, peak instantaneous joint power and both peak sagittal and frontal plane joint moments, the former being linked to muscle force required for body support and the latter being linked to lateral stability and ligament loading, and thus may be more closely related to joint injury mechanisms (Besier et al., 2001a; Besier et al., 2001b). We assessed how these mechanical parameters changed with respect to the stance and swing phases of gait independently (we only included those variables that represent the major action at the joints in these phases).

Example traces of the group mean joint moment and instantaneous power curves during walking and running at the WRT speeds are presented in Fig. 1. The peak moments and peak instantaneous powers during the stance and swing phases of walking and running that are used in our analyses are identified, as are the joint power bursts that define our joint work variables. Results from the speed/gait series analyses are presented below and group mean data are provided for each variable at each speed in the supplementary material (Tables S1–S3, Figs S1–S3).

Mechanical work

The majority of joint work variables exhibited a statistically significant increase with speed, with the exception of the ankle stance negative work (see Table 1 for ANOVA statistics including speed main effects from both two-way and one-way post hoc analyses). The individual speeds identified as being significantly different between walking and running (a priori and post hoc analyses) are identified in Fig. 2.

No mechanical work parameters satisfied all the statistical requirements of the effort–load hypothesis of the WRT (see Materials and methods). The line of best fit for the positive ankle swing work (Fig. 2D) during walking was found to increase above that of running at the WRT speed but did not exhibit a statistically higher value during walking compared with running at the WRT. Table 1 details the ANOVA results and the work variables that satisfied our definition of a general change between walking and running (gait main effect) but which were not identified as satisfying the effort–load hypothesis of gait transition.

Fig. 1.

Joint powers and moments at the WRT. Grouped average moments (A–F) and instantaneous powers (G–L) for the ankle (A,B and G,H), knee (C,D and I,J) and hip (E,F and K,L) joints when walking (A,C,E and G,I,K) and running (B,D,F and H,J,L) at the walk-to-run transition (WRT). The shaded area on the joint moment graphs (A–F) represents the s.d. of the group means. The dark and light shading on the power graphs (G–L) represent positive and negative joint work, respectively, with group s.d. omitted for clarity. AM1 and AM2 (A,B) indicate the peak ankle stance-phase dorsiflexion and plantarflexion moments, respectively. The peak swing-phase moments have not been indicated on the graph because of the scale. KM1 and KM2 indicate the peak knee stance-phase extension and flexion moments, respectively, while KM3 and KM4 indicate the peak knee swing-phase extension and flexion moments, respectively (C,D). HM1 and HM2 represent the peak hip stance-phase extension and flexion moments, respectively, while HM3 and HM4 represent the peak hip swing-phase flexion and extension moments (E,F). AP1 and AP2 represent the peak positive and negative ankle stance-phase power (G,H). The swing-phase power has been not been indicated because of the scale. KP1 indicates the peak positive knee stance-phase power, and KP2 and KP3 indicate the peak negative knee stance-phase and swing-phase power, respectively (I,J). HP1 and HP2 represent the peak positive and negative hip stance-phase power, respectively, and HP3 and HP4 represent the peak positive hip swing-phase power during flexion and extension, respectively. The negative knee swing-phase flexion work (performed during net flexion moments in terminal swing) was computed from the KP3 power burst; the positive hip stance-phase extension work and negative flexion work (performed during net extension and flexion moments, respectively) were computed from the HP1 and HP2 power bursts, respectively; the positive hip swing-phase flexion work (performed during net flexion moments in early swing) was computed from the HP3 power burst.

Fig. 1.

Joint powers and moments at the WRT. Grouped average moments (A–F) and instantaneous powers (G–L) for the ankle (A,B and G,H), knee (C,D and I,J) and hip (E,F and K,L) joints when walking (A,C,E and G,I,K) and running (B,D,F and H,J,L) at the walk-to-run transition (WRT). The shaded area on the joint moment graphs (A–F) represents the s.d. of the group means. The dark and light shading on the power graphs (G–L) represent positive and negative joint work, respectively, with group s.d. omitted for clarity. AM1 and AM2 (A,B) indicate the peak ankle stance-phase dorsiflexion and plantarflexion moments, respectively. The peak swing-phase moments have not been indicated on the graph because of the scale. KM1 and KM2 indicate the peak knee stance-phase extension and flexion moments, respectively, while KM3 and KM4 indicate the peak knee swing-phase extension and flexion moments, respectively (C,D). HM1 and HM2 represent the peak hip stance-phase extension and flexion moments, respectively, while HM3 and HM4 represent the peak hip swing-phase flexion and extension moments (E,F). AP1 and AP2 represent the peak positive and negative ankle stance-phase power (G,H). The swing-phase power has been not been indicated because of the scale. KP1 indicates the peak positive knee stance-phase power, and KP2 and KP3 indicate the peak negative knee stance-phase and swing-phase power, respectively (I,J). HP1 and HP2 represent the peak positive and negative hip stance-phase power, respectively, and HP3 and HP4 represent the peak positive hip swing-phase power during flexion and extension, respectively. The negative knee swing-phase flexion work (performed during net flexion moments in terminal swing) was computed from the KP3 power burst; the positive hip stance-phase extension work and negative flexion work (performed during net extension and flexion moments, respectively) were computed from the HP1 and HP2 power bursts, respectively; the positive hip swing-phase flexion work (performed during net flexion moments in early swing) was computed from the HP3 power burst.

Table 1.

Joint work

Joint work
Joint work

Peak instantaneous joint powers

The majority of instantaneous joint power variables exhibited a main effect of speed (ANOVA, Table 2), including both power generation and absorption (Fig. 3). The individual speeds identified as being significantly different between walking and running (a priori and post hoc analyses) are identified in Fig. 3. From these statistical analyses, the peak positive hip stance power and the peak positive ankle swing power were the only variables identified that satisfied all the statistical requirements of the effort–load hypothesis of the WRT (Fig. 3C,D). The line of best fit for the peak positive hip swing (flexion) power (Fig. 3F) during walking was found to increase above that of running close to the WRT speed, but did not exhibit a statistically higher value during walking compared with running at the WRT.

The summary of ANOVA results including gait main effects and gait–speed interaction effects on peak powers, and the peak power variables that satisfied our definition of a general change between walking and running (gait main effect) but that did not satisfying the effort–load hypothesis of gait transition are outlined in Table 2.

Peak joint moments

The majority of peak joint moment variables increased with speed with the only exception being the peak knee stance abduction moments (Fig. 4G) (see Table 3 for ANOVA statistics including speed main effects from both two-way and one-way post hoc analyses). The individual speeds identified as being significantly different between walking and running (a priori and post hoc analyses) are identified on Fig. 4.

Fig. 2.

Joint work during walking and running plotted against speed. Data are means ± s.d. for walking (filled symbols) and running (open symbols). Speed is expressed as a percentage of the WRT speed; 100% is depicted by the vertical dashed line. (A–C) Positive and negative ankle (A), negative knee (B) and positive hip stance extension work and negative flexion work (performed during net extension and flexion moments, respectively; C). (D–F) Positive ankle swing dorsiflexion work (performed during net dorsiflexion moments in swing; D), negative knee swing flexion work (performed during net flexion moments in terminal swing; E) and positive hip swing flexion work (performed during net flexion moments in early swing; F). Circles and squares represent peak positive and negative work, respectively. Speeds where joint work during walking was significantly greater than that during running (a priori tests). *Significant differences between walking and running established in post hoc analyses (run only when an interaction effect between gait and speed was found). If both a priori and post hoc significant differences were found, only the asterisk is labelled.

Fig. 2.

Joint work during walking and running plotted against speed. Data are means ± s.d. for walking (filled symbols) and running (open symbols). Speed is expressed as a percentage of the WRT speed; 100% is depicted by the vertical dashed line. (A–C) Positive and negative ankle (A), negative knee (B) and positive hip stance extension work and negative flexion work (performed during net extension and flexion moments, respectively; C). (D–F) Positive ankle swing dorsiflexion work (performed during net dorsiflexion moments in swing; D), negative knee swing flexion work (performed during net flexion moments in terminal swing; E) and positive hip swing flexion work (performed during net flexion moments in early swing; F). Circles and squares represent peak positive and negative work, respectively. Speeds where joint work during walking was significantly greater than that during running (a priori tests). *Significant differences between walking and running established in post hoc analyses (run only when an interaction effect between gait and speed was found). If both a priori and post hoc significant differences were found, only the asterisk is labelled.

Table 2.

Peak instantaneous joint powers

Peak instantaneous joint powers
Peak instantaneous joint powers

The peak hip stance flexion moment satisfied the criteria for the effort–load hypothesis of the WRT (Fig. 4C, positive values; Table 3). The line of best fit for the peak hip stance extension moment (Fig. 4C, negative values; Table 3) increased above that of running at the WRT, but did not exhibit a significantly larger value for walking compared with running at the WRT after the Benjamini correction for multiple comparisons. The other peak joint moments did not satisfy the criteria for the effort–load hypothesis of the WRT. The summary of ANOVA results including gait main effects and gait–speed interaction effects on peak joint moments, and the peak joint moment variables that satisfied our definition of a general change between walking and running (gait main effect) but that did not satisfying the effort–load hypothesis of gait transition are outlined in Table 3.

Humans and other terrestrial animals spontaneously switch from a walking to running gait as they increase locomotor speed. This study aimed to determine the extent to which discrete mechanical variables at the individual joints that can affect muscular effort and musculoskeletal loads are associated with the WRT in humans. The majority of the examined joint parameters showed a general marked increase between walking and running at speeds where both walking and running are possible. However, out of the parameters examined, only three satisfied our statistical criteria of the effort–load hypothesis of the WRT. Mechanical parameters that provide an acute signal (peak moments and peak power) were more strongly associated with the gait transition than mechanical work parameters that reflect the muscular function across a portion of the stride, and were present both at the ankle and at the hip.

Fig. 3.

Peak instantaneous joint power during walking and running plotted against speed. Data are means ± s.d. for walking (filled symbols) and running (open symbols). Speed is expressed as a percentage of the WRT speed; 100% is depicted by the vertical dashed line. (A–C) Peak positive and negative ankle (A), negative knee (B) and positive and negative hip stance powers (C). (D–F) Peak positive ankle swing power (D), peak negative knee swing power (E) and peak positive hip swing power in flexion (F). Circles and squares represent peak positive and negative power, respectively. For C and D, the variable satisfied the statistical criteria of the effort–load hypothesis of the WRT. Speeds where joint power during walking was significantly greater than that during running (a priori tests). *Significant differences between walking and running established in post hoc analyses (run only when an interaction effect between gait and speed was found). If both a priori and post hoc significant differences were found, only the asterisk is labelled.

Fig. 3.

Peak instantaneous joint power during walking and running plotted against speed. Data are means ± s.d. for walking (filled symbols) and running (open symbols). Speed is expressed as a percentage of the WRT speed; 100% is depicted by the vertical dashed line. (A–C) Peak positive and negative ankle (A), negative knee (B) and positive and negative hip stance powers (C). (D–F) Peak positive ankle swing power (D), peak negative knee swing power (E) and peak positive hip swing power in flexion (F). Circles and squares represent peak positive and negative power, respectively. For C and D, the variable satisfied the statistical criteria of the effort–load hypothesis of the WRT. Speeds where joint power during walking was significantly greater than that during running (a priori tests). *Significant differences between walking and running established in post hoc analyses (run only when an interaction effect between gait and speed was found). If both a priori and post hoc significant differences were found, only the asterisk is labelled.

Fig. 4.

Peak joint moments during walking and running plotted against speed. Data are means ± s.d. for walking (filled symbols) and running (open symbols). Speed is expressed as a percentage of the WRT speed; 100% is depicted by the vertical dashed line. (A–C) Peak ankle stance-phase plantarflexion moments (A), peak knee stance-phase extension (B) and hip stance-phase flexion and extension moments (C). (D–F) Peak ankle swing-phase dorsiflexion moments (D), peak knee swing-phase flexion moments (E) and hip swing-phase flexion and extension moments (F). (G,H) Peak knee (G) and hip (H) stance-phase abduction moments. Positive values (circles) represent net joint flexion moments (dorsiflexion at the ankle); negative values (squares) represent net joint extension moments (plantarflexion at the ankle, abduction at the knee and hip). For C, the variable satisfied the statistical criteria of the effort–load hypothesis of the WRT. Speeds where walking was significantly greater than running (a priori tests). *Significant differences between walking and running established in post hoc analyses (run only when an interaction effect between gait and speed was found). If both a priori and post hoc significant differences were found, only the asterisk is labelled.

Fig. 4.

Peak joint moments during walking and running plotted against speed. Data are means ± s.d. for walking (filled symbols) and running (open symbols). Speed is expressed as a percentage of the WRT speed; 100% is depicted by the vertical dashed line. (A–C) Peak ankle stance-phase plantarflexion moments (A), peak knee stance-phase extension (B) and hip stance-phase flexion and extension moments (C). (D–F) Peak ankle swing-phase dorsiflexion moments (D), peak knee swing-phase flexion moments (E) and hip swing-phase flexion and extension moments (F). (G,H) Peak knee (G) and hip (H) stance-phase abduction moments. Positive values (circles) represent net joint flexion moments (dorsiflexion at the ankle); negative values (squares) represent net joint extension moments (plantarflexion at the ankle, abduction at the knee and hip). For C, the variable satisfied the statistical criteria of the effort–load hypothesis of the WRT. Speeds where walking was significantly greater than running (a priori tests). *Significant differences between walking and running established in post hoc analyses (run only when an interaction effect between gait and speed was found). If both a priori and post hoc significant differences were found, only the asterisk is labelled.

Which discrete joint mechanical variables might affect the WRT?

Previous studies have identified discrete joint kinematics (Hreljac, 1995; Minetti et al., 1994), joint moments (Prilutsky and Gregor, 2001) and peak joint power (Hreljac et al., 2008) that may be linked to the WRT. Because these studies either focused on a sub-set of joints or mechanical variables separately, it is difficult to deduce from them whether each parameter is equally important in the WRT. The current study's comprehensive inverse dynamic analyses indicate that various joint-level mechanical variables are linked to the WRT, including both peak joint power and peak joint moments. This might result, in part, because these parameters are each associated with the same muscular actions. For example, the required increase in peak positive hip stance power (Fig. 3C) may occur partly because of the increase in peak hip stance flexion moments (Fig. 4C, positive values), as joint power is a function of joint moment.

Table 3.

Peak joint moments

Peak joint moments
Peak joint moments

It is nevertheless interesting that all of the mechanical variables satisfying the effort–load hypothesis of the WRT identified in this study (peak joint power and peak joint moments) are those that can provide an acute physiological signal. This is consistent with the theory that gait transition occurs spontaneously, is initiated within a single stride (Segers et al., 2013) and is triggered by discrete variables that can be acutely sensed. In contrast, joint work reflects the joint's mechanical function over a portion of the stride. Similar to the overall metabolic rate, the longer time over which these variables fluctuate might preclude them from causing the transition step itself [although AMP-activated protein kinase (Winder, 2001) and/or other cellular mechanisms (Clanton et al., 2013) may provide acute sensing of the energy/work state of a muscle]. However, the increase in joint mechanical work at the ankle and hip during fast walking above that required for running offers a reasonable explanation, together with instantaneous peaks in joint power and moments, for why humans avoid walking at these speeds and instead choose to either walk slower or run faster. They might also help explain, in part, the higher cost of walking compared with running at fast walking speeds.

It is likely that the variables consistent with the effort–load hypothesis identified in this study reflect local muscle-level stimuli. These local joint-level effects were, however, not restricted to a single phase of the gait cycle. Variables linked to the WRT were identified both in the swing phase (peak positive ankle dorsiflexion power) and in the stance phase (peak hip power and flexion moments). Furthermore, the peak hip flexion moments and peak hip power occur at different times during stance. Given that gait transition is likely triggered within a short time frame, it remains unclear whether one of these variables is more strongly associated with the transition, whether each independently affects the WRT at different times, or whether they possibly act cumulatively.

Evidence for both ankle- and hip-based mechanics influencing the WRT

Several studies have identified effort (Hreljac, 1995; Hreljac et al., 2001; Hreljac et al., 2008) or fatigue (Malcolm et al., 2009; Segers et al., 2007) in the ankle dorsiflexors during the swing phase as a main factor influencing the WRT. The findings of this study support this hypothesis, providing further evidence that ankle dorsiflexion in swing is a key mechanism affecting the transition from walking to running. Indeed, this study found that the peak swing-phase dorsiflexion power (Fig. 3D) satisfies the effort–load hypothesis of the WRT transition, with peak power values during walking exceeding those of running at the WRT speed, similar to that of MacLeod et al. (MacLeod et al., 2014). Our results for ankle swing-phase dorsiflexion work (Fig. 2D) further indicate that ankle dorsiflexion mechanics become unfavourable at walking speeds beyond the WRT. It is plausible that the high peak joint power and work are the underlying mechanisms leading to dorsiflexion fatigue at fast walking speeds and that the reduction of these variables upon switching to running mitigates the fatigue and effort sensed in the muscles.

It has also been argued that stance-phase plantarflexion is a factor contributing to the WRT. Neptune and Sasaki (Neptune and Sasaki, 2005), and more recently Arnold et al. (Arnold et al., 2013), showed in simulation studies that the ability to generate force is compromised at fast walking speeds because of sub-optimal force–length–velocity characteristics of the triceps surae muscles. Switching to running allowed these muscles to function at a more favourable length and velocity, thus increasing their force capacity. More direct evidence for this has recently been found by combining ultrasound imaging and gait analysis, in which the velocity of the medial gastrocnemius fascicles has been identified as limiting force and power production during walking at the WRT speed (Farris and Sawicki, 2012b). Our speed-series joint-level analysis provides further support for this theory, whereby a plateau in peak instantaneous positive power and plantarflexion moments occurs during walking above the WRT (Fig. 3A, Fig. 4A). The switch to running leads to a marked increase in ankle joint power and plantarflexion moments (Fig. 3A, Fig. 4A) and is in agreement with the higher estimated triceps surae muscle force and power during running at the WRT speed (Neptune and Sasaki, 2005; Farris and Sawicki, 2012b; Arnold et al., 2013). We cannot rule out, however, that the combination of the limitation in plantarflexion torque and dorsiflexion fatigue may contribute to the loss of local dynamic stability at the ankle joint when walking at speeds beyond the WRT (Jordan et al., 2009), and may together reflect a dynamic systems determinant of the WRT as opposed to a discrete ankle joint determinant.

In addition to the aforementioned ankle mechanisms, our study also identified novel discrete hip-based mechanical variables that are associated with the WRT. The increase in peak positive hip stance power and peak flexion moments (Fig. 3C and Fig. 4C positive values) during walking to values above those required for running at the WRT may contribute to the transition to a running gait. It is intriguing to consider whether the sharp increase in the peak positive hip stance-phase power beyond the WRT speed (Fig. 3C) occurs because of the inability of the ankle to produce sufficient stance-phase peak power and joint moments (Fig. 3A and Fig. 4A). In this regard, the trigger for the WRT may, in part, be the result of a cascade of mechanical events beginning at the ankle and leading to an unfavourable compensation at the hip, both of which may provide the critical signal for altering gait. Unlike Prilutsky and Gregor (Prilutsky and Gregor, 2001), we did not observe a clear unfavourable effect on hip swing-phase joint moments at fast walking speeds. This may be due to the moderately smaller range of walking speeds in the present study as a result of the difficulty our subjects had in maintaining faster walking speeds. That hip swing-phase mechanics might influence the WRT was, however, supported by the trend of greater peak hip flexion joint power in walking compared with running above the WRT speed (Fig. 3F).

Together, these results suggest that the switch between walking and running may occur not only to reduce effort in ankle muscles but also in the hip musculature. Some previous support for a hip-based mechanism affecting the WRT can be found from electromyography analyses of the hip muscles (rectus femoris and biceps femoris) at walking and running speeds spanning the WRT speed (Prilutsky and Gregor, 2001). Furthermore, simulation studies in which individual muscle mechanics were predicted for walking and running at and above and below the WRT indicate favourable reductions in peak power and work in hip muscle fibres as a result of switching to running (Sasaki and Neptune, 2006a; Sasaki and Neptune, 2006b), although this has yet to be shown experimentally in humans. Interestingly, the redistribution of average mechanical power from the hip to the ankle that accompanies the reduction in hip work at the WRT has been suggested to contribute to the greater locomotor efficiency of running (Farris and Sawicki, 2012a). This may occur because of the purported greater plantarflexor efficiency during running, and might also contribute to gait selection.

Is injury avoidance a factor in human gait transition?

Whether the elevation in peak hip sagittal plane joint moments (Fig. 4C) would influence the WRT through an injury reduction mechanism [as has been suggested in horses to reduce tendon loads (Farley and Taylor, 1991)] is questionable. The increased hip flexion and extension moments are most likely accompanied by an increase in muscle force that, while potentially increasing the effort of locomotion, are at levels that are not expected to pose any significant musculoskeletal injury risk.

More likely to be linked with injury mechanisms are the frontal plane loads at the knee (Besier et al., 2001a; Besier et al., 2001b) and the hip. The present study provides among the first measurements of the response of non-sagittal loading at these joints across speed. Interestingly, while the abduction loading at the knee and hip increase at the same rate with speed during both walking and running, the load level is higher during running at all speeds where both gaits were analysed (Fig. 4). Running may therefore place the joints at larger risk of injury in general, but there is no evidence from our study that frontal plane loading per se influences the WRT in a manner to specifically reduce injury risk at the WRT speed.

General effect of gait transition

Whether joint mechanical differences between walking and running are general across all speeds where walking and running are possible, or whether they arise due to a speed-related effect has not been extensively examined. Our analysis found that most of the joint parameters in both the sagittal and frontal planes exhibited a marked increase with a shift from walking to running, both below and above the WRT. Furthermore, half of variables lacked an interaction effect between gait and speed (Tables 1, 2 and 3). These findings suggest that, overall, there is an increase in the joint mechanical variables between walking and running that are often unaffected by speed. Interestingly, those variables that were invariant of gait were predominantly found in the swing phase (Tables 1, 2 and 3). That differences in joint mechanics between walking and running are more predominant in the stance phase of gait is consistent with the body centre-of-mass paradigms of walking and running (inverted pendulum versus spring mass), which are dictated primarily by stance dynamics (Saibene and Minetti, 2003). We also found some joint variables at the hip and the ankle that exhibited speed-dependent differences between walking and running. For example, the ankle swing-phase work (Fig. 2D), and the peak positive hip stance- and swing-phase flexion power (Fig. 3C,F) and stance-phase moments (Fig. 4C) all lacked a main effect of gait but exhibited an interaction effect between gait and speed (Tables 1, 2 and 3). In these variables, either there was a difference between gaits only at faster speeds or the differences were in opposite directions above versus below the WRT. In this latter scenario, these variables might help explain both the WRT as well as the run-to-walk transition.

While the general increase in these joint parameter magnitudes alone may not explain the WRT, they may, however, reflect which joint variables are responsible for initiating the transition between stable patterns of coordination (walking/running). As such, while not affecting the WRT per se, they may represent the key variables underpinning a dynamic systems interpretation of the WRT and may be representative of the central motor plan for moving between attractor states and avoiding system instability.

Limitations

Our joint-level inverse dynamic analyses represent the net effect of all muscles and structures that span the joints. As has been outlined previously (Sasaki et al., 2009), inverse dynamic analyses do not necessarily reflect the mechanics of individual muscles. We have not taken into account co-contraction between antagonist muscles, force sharing between synergist muscles or the distribution of work and power between muscle fibres and tendon. It is also important to stress that our criteria for the effort–load hypothesis of the WRT to be satisfied are based on the statistical identification of a unique increase in a parameter during walking compared with running at the WRT. This is the case when a parameter for walking first increases above running at the WRT. Whilst these criteria are designed to detect a unique event at the WRT, we cannot rule out that a variable reaches a critical value that influences gait transition without being uniquely identifiable at the WRT speed using the above criteria. They also do not identify other mechanisms that may influence gait transition, such as a restricted capacity for stance-phase peak ankle power and moments output that were evident in our analysis of joint mechanics (Fig. 3A and Fig. 4A). Finally, our analysis of the WRT does not specifically assess whether non-specific control parameters, which may be reflected in our discrete joint measurements, are responsible for the WRT (rather than the discrete variables themselves).

Conclusion

Discrete joint-level mechanisms at both the ankle and hip that are thought to increase muscular effort have been identified as being associated with the WRT in humans. Of the examined variables, only those that provide an acute signal satisfied our effort–load hypothesis of the WRT. We hypothesize that the WRT in humans is dictated, in part, by a limitation in ankle moment and power generation that results in a compensation at the hip that increases the effort in hip muscles above that which is required during running. Finally, our analyses suggest that the differences in joint mechanics between walking and running in most joint-level parameters are consistent across different speeds where walking and running are both feasible gaits.

Subjects

Eight healthy and recreationally active subjects (N=4 men, N=4 women), with no history of major lower-limb injury were recruited for this study (age 24.8±1.8 years, height 170.0±9.4 cm, mass 69.6±13.2 kg, means ± s.d.). All testing procedures were approved by the Human Research Ethics Committee of the University of Western Australia. Written informed consent was obtained from all subjects.

Walk-to-run gait transition speed

All subjects were accustomed to running on treadmills and were initially familiarized to walk and run on a force plate instrumented split-belt treadmill (Bertec, Columbus, OH, USA). The subject's preferred WRT speed was determined following a protocol similar to those used elsewhere (Hreljac, 1993a; Hreljac, 1993b; Bartlett and Kram, 2008). The subjects walked at randomized speeds ranging from 1.5 to 2.5 m s−1 at intervals of 0.1 m s−1 for ~30 s to 1 min per speed. At each speed, the subjects were instructed to voluntarily choose between walking and running and were allowed to switch between gaits. The slowest speed where the subjects confirmed running to be their ‘most comfortable’ gait was selected as the WRT speed. This was repeated three times and the subject's range and average WRT speed were determined (group average 2.00±0.09 m s−1, ranging from 1.90 to 2.10 m s−1).

Three-dimensional (3D) gait analysis

Experiments were performed on a split belt force plate treadmill (Bertec). Subjects walked at speeds ranging from 30% to 120% of their WRT speed (0.5–2.5 m s−1) and ran at speeds ranging from 80% to 170% (1.5–3.5 m s−1) of their WRT. For both walking and running, speed increments of 10% were analysed with the subject walking/running at a set speed for at least 1 min. Five speeds spanning 80–120% of the WRT were performed for both walking and running. The walking and running speeds were randomized to prevent any order effects. Five strides per subject per speed were analysed where the subject maintained their anterior–posterior and medial–lateral position on the treadmill. Individual strides where gait was spontaneously changed between walking and running were not analysed for this study.

3D ground reaction forces (2000 Hz) and marker position data (250 Hz; 8-camera VICON MX motion-analysis system, Vicon Oxford Metrics, Oxford, UK) were captured and integrated in the Vicon Nexus data acquisition suite (version 1.8.3). Retroreflective markers were placed either on bony landmarks or as part of triad-marker clusters. Single markers were placed on the left and right anterior superior iliac spines, left and right posterior superior iliac spines, left and right head of first and fifth metatarsals, and left and right calcaneus. Triad-marker clusters were placed on the left and right thigh and leg. 3D kinematic and inverse dynamic calculations were computed in Vicon Nexus and Bodybuilder (Vicon Oxford Metrics) in accordance with procedures established previously (Besier et al., 2003), including a functional knee axis and joint centre (mean helical axis) and a functional hip centre (sphere-fit) calculation. Prior to inverse dynamic modelling, marker position data were filtered using a fourth-order, 12 Hz (walking) and 8 Hz (running) low-pass Butterworth filter, with the filter frequencies selected by performing a residual analysis of filtered versus unfiltered data (MATLAB, The MathWorks, Natick, MA, USA). Ground reaction force signals were filtered at the same frequency as the kinematic data to prevent joint moment artifacts (Kristianslund et al., 2012).

Joint moments and instantaneous joint power calculations

Joint moments were expressed in the distal segment anatomical coordinate systems as per Besier et al. (Besier et al., 2003). Only flexion/extension and adduction/abduction joint moments were used for WRT transition analyses because of the high variability in the long-axis joint loading. Joint powers were computed as the net power across all three joint planes. Joint moments and powers were computed for the right leg and both normalized to body mass (N kg−1 and W kg−1, respectively). Numerical values of the peak joint moments and powers were initially identified using MATLAB. Secondary visual inspection identified whether the peak joint moments and instantaneous peak joint powers corresponded to the physiologically relevant phases of the gait cycle (see Fig. 1).

Joint work

Positive and negative joint work () were calculated for the ankle, knee and hip from the positive and negative values of the instantaneous joint power curves (), respectively:
(1)

The work from the individual joints (Eqn 1) was computed for the stance and swing phases independently. Work was normalized to body mass (J kg−1). When clear bursts of positive and/or negative work in the gait phase (stance/swing) were produced separately under net flexion and extension joint moments, both the individual flexion and extension work were computed (Fig. 1). This included positive hip joint work produced during the swing and stance phases, negative work at the hip in stance, negative work at the knee in swing, and the positive work at the ankle in swing. We report the flexion or extension work depending on which parameter represents the primary function of the joint in the gait phase (see Fig. 1).

Testing of the effort–load hypothesis for mechanical variables

Similar to Hreljac (Hreljac, 1993b; Hreljac, 1995), we defined a set of criteria for a mechanical variable to satisfy the effort–load hypothesis of the WRT. Our analysis is based on the hypothesis that the WRT occurs when an ‘undesirable’ discrete variable surpasses a certain threshold value with increasing walking speed and is lowered by switching to running. A mechanical variable satisfied the effort–load hypothesis of the WRT if it increased with walking speed and first became larger during walking compared with running at the WRT speed such that switching gaits resulted in a favourable reduction in that variable. This definition is consistent with the criteria of Hreljac (Hreljac, 1993b; Hreljac, 1995) whereby a mechanical ‘trigger’ for the WRT is required to increase with walking speed and undergo an abrupt decrease upon switching to running, but with the addition that the stimulus for lowering the variable first occurs at the WRT speed. To assess these criteria of the effort–load hypothesis, we first examined the relationships between gait and speed using a two-way (gait: walking and running) by five repeated-measures (speed: 80%, 90%, 100%, 110% and 120% of WRT speed) ANOVA (SPSS version 21). Significant main effects of speed and gait were assessed (P<0.05) and adjusted for the false discovery rate arising from multiple comparisons (Benjamini adjustment) (Benjamini and Hochberg, 1995). The variable was deemed to increase with walking speed when a main effect of speed was detected, or in the event that an interaction effect between speed and gait existed, when a post hoc one-way ANOVA limited to walking data exhibited a main effect of speed. Secondly, a priori tests were used to establish those speeds where walking data were significantly greater than the running data. Because our criteria asked specifically when the walking data first became greater than the running data, we performed a one-tailed paired sample t-test at each speed (P<0.05 with Benjamini multiple comparison adjustment).

General joint mechanical changes between walking and running at different speeds

We also determined those variables that were generally affected by switching between a walking and running gait, but which did not satisfy the effort–load hypothesis per se. These variables demonstrated a shift between walking and running gaits at the WRT and other speeds where walking and running were compared. In order to determine whether gait had a general effect on the biomechanical variables tested, we used the two-way repeated-measures ANOVA described above to assess gait and speed main and interaction effects (P<0.05 with Benjamini multiple comparison adjustment). Variables that exhibited a main effect of gait were deemed to undergo a generalized modification between walking and running. When an interaction effect between gait and speed was found, post hoc tests were used to identify significant differences between the walking and running data at individual speeds (two-tailed paired sample t-test; P<0.05 with Benjamini multiple comparison adjustment).

The authors would like to thank James J. Dunne, Sarah M. Stearne, Luqman Aziz and Kirsty A. McDonald for providing assistance with data collection and analysis, Tony Roby and Jared Issott for assistance with equipment development and maintenance, and all of the participants for volunteering their time to this study.

Funding

This work was supported by an Australian Postgraduate Award scholarship to N.J.P.

Arnold
E. M.
,
Hamner
S. R.
,
Seth
A.
,
Millard
M.
,
Delp
S. L.
(
2013
).
How muscle fiber lengths and velocities affect muscle force generation as humans walk and run at different speeds
.
J. Exp. Biol.
216
,
2150
-
2160
.
Bartlett
J. L.
,
Kram
R.
(
2008
).
Changing the demand on specific muscle groups affects the walk-run transition speed
.
J. Exp. Biol.
211
,
1281
-
1288
.
Benjamini
Y.
,
Hochberg
Y.
(
1995
).
Controlling the false discovery rate: a practical and powerful approach to multiple testing
.
J. R. Stat. Soc.
57
,
289
-
300
.
Besier
T. F.
,
Lloyd
D. G.
,
Ackland
T. R.
,
Cochrane
J. L.
(
2001a
).
External loading of the knee joint during running and cutting maneuvers
.
Med. Sci. Sports Exerc.
33
,
1176
-
1181
.
Besier
T. F.
,
Lloyd
D. G.
,
Ackland
T. R.
,
Cochrane
J. L.
(
2001b
).
Anticipatory effects on knee joint loading during running and cutting maneuvers
.
Med. Sci. Sports Exerc.
33
,
1176
-
1181
.
Besier
T. F.
,
Sturnieks
D. L.
,
Alderson
J. A.
,
Lloyd
D. G.
(
2003
).
Repeatability of gait data using a functional hip joint centre and a mean helical knee axis
.
J. Biomech.
36
,
1159
-
1168
.
Biewener
A. A.
,
Taylor
C. R.
(
1986
).
Bone strain: a determinant of gait and speed?
J. Exp. Biol.
123
,
383
-
400
.
Clanton
T. L.
,
Hogan
M. C.
,
Gladden
L. B.
(
2013
).
Regulation of cellular gas exchange, oxygen sensing, and metabolic control
.
Compr. Physiol.
3
,
1135
-
1190
.
Farley
C. T.
,
Taylor
C. R.
(
1991
).
A mechanical trigger for the trot-gallop transition in horses
.
Science
253
,
306
-
308
.
Farris
D. J.
,
Sawicki
G. S.
(
2012a
).
The mechanics and energetics of human walking and running: a joint level perspective
.
J. R. Soc. Interface
9
,
110
-
118
.
Farris
D. J.
,
Sawicki
G. S.
(
2012b
).
Human medial gastrocnemius force-velocity signaling shifts with locomotion speed and gait
.
Proc. Natl. Acad. Sci. USA
109
,
977
-
982
.
Hoyt
D.
,
Taylor
C.
(
1981
).
Gait and the energetics of locomotion in horses
.
Nature
292
,
239
-
240
.
Hreljac
A.
(
1993a
).
Preferred and energetically optimal gait transition speeds in human locomotion
.
Med. Sci. Sports Exerc.
25
,
1158
-
1162
.
Hreljac
A.
(
1993b
).
Determinants of the gait transition speed during human locomotion: kinetic factors
.
Gait Posture
1
,
217
-
223
.
Hreljac
A.
(
1995
).
Determinants of the gait transition speed during human locomotion: kinematic factors
.
J. Biomech.
28
,
669
-
677
.
Hreljac
A.
,
Arata
A.
,
Ferber
R.
,
Mercer
J. A.
,
Row
B.
(
2001
).
An electromyographical analysis of the role of dorsiflexors on the gait transition during human locomotion
.
J. Appl. Biomech.
17
,
287
-
296
.
Hreljac
A.
,
Imamura
R. T.
,
Escamilla
R. F.
,
Edwards
W. B.
,
MacLeod
T.
(
2008
).
The relationship between joint kinetic factors and the walk-run gait transition speed during human locomotion
.
J. Appl. Biomech.
24
,
149
-
157
.
Jordan
K.
,
Challis
J. H.
,
Cusumano
J. P.
,
Newell
K. M.
(
2009
).
Stability and the time-dependent structure of gait variability in walking and running
.
Hum. Mov. Sci.
28
,
113
-
128
.
Kristianslund
E.
,
Krosshaug
T.
,
van den Bogert
A. J.
(
2012
).
Effect of low pass filtering on joint moments from inverse dynamics: implications for injury prevention
.
J. Biomech.
45
,
666
-
671
.
MacLeod
T. D.
,
Hreljac
A.
,
Imamura
R.
(
2014
).
Changes in the preferred transition speed with added mass to the foot
.
J. Appl. Biomech.
30
,
95
-
103
.
Malcolm
P.
,
Segers
V.
,
Van Caekenberghe
I.
,
De Clercq
D.
(
2009
).
Experimental study of the influence of the m. tibialis anterior on the walk-to-run transition by means of a powered ankle-foot exoskeleton
.
Gait Posture
29
,
6
-
10
.
Margaria
R.
,
Cerretelli
P.
,
Aghemo
P.
,
Sassi
G.
(
1963
).
Energy cost of running
.
J. Appl. Physiol.
18
,
367
-
370
.
Mercier
J.
,
Le Gallais
D.
,
Durand
M.
,
Goudal
C.
,
Micallef
J. P.
,
Préfaut
C.
(
1994
).
Energy expenditure and cardiorespiratory responses at the transition between walking and running
.
Eur. J. Appl. Physiol.
69
,
525
-
529
.
Minetti
A. E.
,
Ardigò
L. P.
,
Saibene
F.
(
1994
).
The transition between walking and running in humans: metabolic and mechanical aspects at different gradients
.
Acta Physiol. Scand.
150
,
315
-
323
.
Neptune
R. R.
,
Sasaki
K.
(
2005
).
Ankle plantar flexor force production is an important determinant of the preferred walk-to-run transition speed
.
J. Exp. Biol.
208
,
799
-
808
.
Prilutsky
B. I.
,
Gregor
R. J.
(
2001
).
Swing- and support-related muscle actions differentially trigger human walk-run and run-walk transitions
.
J. Exp. Biol.
204
,
2277
-
2287
.
Rubenson
J.
,
Heliams
D. B.
,
Lloyd
D. G.
,
Fournier
P. A.
(
2004
).
Gait selection in the ostrich: mechanical and metabolic characteristics of walking and running with and without an aerial phase
.
Proc. Biol. Sci.
271
,
1091
-
1099
.
Saibene
F.
,
Minetti
A. E.
(
2003
).
Biomechanical and physiological aspects of legged locomotion in humans
.
Eur. J. Appl. Physiol.
88
,
297
-
316
.
Sasaki
K.
,
Neptune
R. R.
(
2006a
).
Differences in muscle function during walking and running at the same speed
.
J. Biomech.
39
,
2005
-
2013
.
Sasaki
K.
,
Neptune
R. R.
(
2006b
).
Muscle mechanical work and elastic energy utilization during walking and running near the preferred gait transition speed
.
Gait Posture
23
,
383
-
390
.
Sasaki
K.
,
Neptune
R. R.
,
Kautz
S. A.
(
2009
).
The relationships between muscle, external, internal and joint mechanical work during normal walking
.
J. Exp. Biol.
212
,
738
-
744
.
Segers
V.
,
Lenoir
M.
,
Aerts
P.
,
De Clercq
D.
(
2007
).
Influence of M. tibialis anterior fatigue on the walk-to-run and run-to-walk transition in non-steady state locomotion
.
Gait Posture
25
,
639
-
647
.
Segers
V.
,
De Smet
K.
,
Van Caekenberghe
I.
,
Aerts
P.
,
De Clercq
D.
(
2013
).
Biomechanics of spontaneous overground walk-to-run transition
.
J. Exp. Biol.
216
,
3047
-
3054
.
Tseh
W.
,
Bennett
J.
,
Caputo
J. L.
,
Morgan
D. W.
(
2002
).
Comparison between preferred and energetically optimal transition speeds in adolescents
.
Eur. J. Appl. Physiol.
88
,
117
-
121
.
Watson
R. R.
,
Rubenson
J.
,
Coder
L.
,
Hoyt
D. F.
,
Propert
M. W. G.
,
Marsh
R. L.
(
2011
).
Gait-specific energetics contributes to economical walking and running in emus and ostriches
.
Proc. Biol. Sci.
278
,
2040
-
2046
.
Winder
W. W.
(
2001
).
Energy-sensing and signaling by AMP-activated protein kinase in skeletal muscle
.
J. Appl. Physiol.
91
,
1017
-
1028
.

Competing interests

The authors declare no competing financial interests.

Supplementary information