The role of trunk orientation during uneven running is not well understood. This study compared the running mechanics during the approach step to and the step down for a 10 cm expected drop, positioned halfway through a 15 m runway, with that of the level step in 12 participants at a speed of 3.5 m s−1 while maintaining self-selected (17.7±4.2 deg; mean±s.d.), posterior (1.8±7.4 deg) and anterior (26.6±5.6 deg) trunk leans from the vertical. Our findings reveal that the global (i.e. the spring-mass model dynamics and centre-of-mass height) and local (i.e. knee and ankle kinematics and kinetics) biomechanical adjustments during uneven running are specific to the step nature and trunk posture. Unlike the anterior-leaning posture, running with a posterior trunk lean is characterized by increases in leg angle, leg compression, knee flexion angle and moment, resulting in a stiffer knee and a more compliant spring-leg compared with the self-selected condition. In the approach step versus the level step, reductions in leg length and stiffness through the ankle stiffness yield lower leg force and centre-of-mass position. Contrariwise, significant increases in leg length, angle and force, and ankle moment, reflect in a higher centre-of-mass position during the step down. Plus, ankle stiffness significantly decreases, owing to a substantially increased leg compression. Overall, the step down appears to be dominated by centre-of-mass height changes, regardless of having a trunk lean. Observed adjustments during uneven running can be attributed to anticipation of changes to running posture and height. These findings highlight the role of trunk posture in human perturbed locomotion relevant for the design and development of exoskeleton or humanoid bipedal robots.

Running is a low cost and easy accessibility form of physical activity that enhances health and increases longevity (Lee et al., 2017). The popularity of running continues to grow worldwide. In 2019 alone, millions of recreational participants around the world covered 1.3 billion total miles with an average distance of 4.1 miles per run (https://2019.strava.com/community/en-us/). Running outdoors in a natural environment often entails a frequent negotiation of the terrain irregularities such as variations in ground compliance, slipperiness or substrate height. In running over changing surfaces, human runners appear to use spring-mass dynamics to help passively stabilize their locomotion (Blickhan, 1989; McMahon and Cheng, 1990). The spring-mass model (Blickhan, 1989) is composed of a mass-less spring and the body (represented by a point mass), and is simply described by: leg stiffness (kleg), leg orientation (θTD) and leg length (lTD) at touchdown (TD) (Blickhan, 1989; McMahon and Cheng, 1990). Albeit the model does not account for the trunk movements, this is a simplistic approach to describe the basic mechanics of human locomotion.

The mechanics of the body's interaction with the ground is influenced by kleg, which is indicative of the average stiffness of the overall musculoskeletal system during the ground contact phase. For running on uneven ground, θTD and the kleg appear to be adjusted relative to the nature of changes in the substrate height. This involves a flatter θTD and a more compliant leg when stepping up, but a steeper θTD and a stiffer leg when stepping down. For the latter situation, runners reduce kleg in the preceding step, possibly to lower the centre of mass (CoM) in preparation for the safe negotiation of expected changes in substrate height, e.g. a drop (Ernst et al., 2014). When traversing an obstacle, kleg decreases in the penultimate step and increases in the final step prior to the barrier in order to accelerate the CoM upward and take-off velocity (Mauroy et al., 2014). Therefore, the fine adjustments of the leg properties not only in the perturbed contact but also in the preceding contacts (Mauroy et al., 2014; Müller and Blickhan, 2010; Müller et al., 2012a, 2010) potentially enhance stability in response to ground perturbations.

The global kleg is determined by the combination of the local stiffness of the joints and the leg geometry at touchdown (Farley et al., 1998; Farley and Morgenroth, 1999; Günther and Blickhan, 2002). Adjustments in both joint stiffness and kleg during the stance phase seem to be governed by initial conditions such as lTD, θTD and landing velocity (Blickhan et al., 2007; Farley et al., 1998; Farley and Morgenroth, 1999; Grimmer et al., 2008; Mauroy et al., 2014), which are previously geared by a swing-leg retraction throughout the second half of the swing phase prior to landing (Seyfarth et al., 2003). Therefore, muscle pre-activation plays a major role in adjusting the leg posture in preparation for traversing, for example, an obstacle. While adjustments in leg mechanics have been widely proven as key to the accommodation of uneven ground, in this context, little is known about the role of trunk posture during uneven running. Given the ∼50% contribution of the trunk segment to total body mass (De Leva, 1996), mechanical demands of the lower limbs can be potentially influenced by even slight alterations in trunk kinematics. Moreover, analysis of an altered mechanical behavior of the lower limbs evoked by trunk posture modifications against an expected substrate height change (i.e. induced deviations from the cyclic motion) could help the understanding of human locomotion.

There are a limited number of studies that have addressed the influence of an anterior trunk lean on knee joint energetics (Arendse et al., 2004; Teng and Powers, 2015), knee muscle activity (Teng and Powers, 2016), patellofemoral joint stress (dos Santos et al., 2019; Teng and Powers, 2014) and impact loading (Huang et al., 2019) during running over uniform ground surfaces. For running on an uneven terrain, one study showed that runners adopt a more crouched posture and larger kleg (Voloshina and Ferris, 2015), possibly to facilitate adaptations to an unfamiliar environment. As recently shown by a study on recreational runners (Shih et al., 2019), modifications in the trunk orientation are achievable with simple postural instructions following 1 month of training. However, it is unknown how the modification of sagittal plane trunk posture can potentially impact lower limb operation during uneven running. In walking with increasing trunk flexion angles against changes in substrate height, we previously demonstrated that the trunk plays a functional role through extending its angle during stepping down, which is necessary to control the angular momentum of the whole body (AminiAghdam and Blickhan, 2018; Aminiaghdam et al., 2017a). This contributed to an observed elevated CoM trajectory during the step down, presumably to ease the drop negotiation.

An overall objective of our study was to provide insight into the contribution of sagittal plane trunk lean to running mechanics against an expected change in substrate height. Anticipation of changes to running pattern can facilitate stability and manoeuvrability (Dhawale et al., 2019; Müller et al., 2015; Qiao and Jindrich, 2012). We expected runners to exhibit step-to-step modulations in local and global running mechanics to account for expected changes in the substrate height. We hypothesized that runners would exhibit strategies in the global and local running mechanics that are task specific, owing to the anticipatory nature of changes to posture and substrate height. In particular, we expected running with a posterior versus anterior trunk lean to be associated with greater compensatory kinematic adjustments, as a result of a further shift of the ground reaction force (GRF) line of action from the knee joint axis of rotation, resulting in augmented changes in the magnitude of lower limb local and global kinetic parameters.

List of symbols and abbreviations

     
  • ATL

    anterior trunk lean

  •  
  • BW

    body weight

  •  
  • CoMTD, CoMTO

    normalized centre of mass at touchdown and at toe-off

  •  
  • Fz

    peak normalized vertical ground reaction force

  •  
  • GRF

    ground reaction force

  •  
  • kankle

    normalized ankle stiffness

  •  
  • kknee

    normalized knee stiffness

  •  
  • kleg

    dimensionless leg stiffness

  •  
  • lTD, lTO

    normalized leg length at touchdown and at toe-off

  •  
  • Mankle

    peak normalized ankle extension moment

  •  
  • Mknee

    peak normalized knee extension moment

  •  
  • PTL

    posterior trunk lean

  •  
  • STL

    self-selected trunk lean

  •  
  • control step

  •  
  • approach step

  •  
  • step↓

    step down

  •  
  • t

    time

  •  
  • ΔLmax

    dimensionless maximum leg compression

  •  
  • θleg,TD, θleg,TO

    leg orientation at touchdown and at toe-off

Participants

A convenience sample of 12 (six females, six males) volunteer recreational runners (mean±s.d. age 28.5±5.7 years, body mass 65.5±8.6 kg, height 168.9±6.4 cm) gave written, informed consent to participate in the study. The experimental protocol was approved by the local Ethics Committee of Friedrich-Schiller-University Jena and conducted according to the Declaration of Helsinki.

Experimental design and protocol

Data were collected on a 15 m long instrumented track, with two consecutive force plates embedded halfway along (1000 Hz; 9281B, 9287BA, Kistler, Winterthur, Switzerland) and 12 infrared system cameras (250 Hz; MCU1000, Qualisys, Gothenburg, Sweden). Force plates were set at a distance of one step from each other, allowing step lengths ranging from 1.40 to 2.30 m. Kinematics and GRF data were synchronized using a Kistler external trigger and BioWare data acquisition software (Kistler). A 12 body segment model was defined using 19 reflective markers. Joint coordinate standards of the International Society of Biomechanics (Wu et al., 2002, 2005) were applied. Trunk angle was defined by the angle sustained by the line connecting the midpoint between the L5–S1 junction (L5) and the seventh cervical spinous process (C7) with respect to the vertical. Following running with self-selected trunk lean (STL; mean±s.d. 17.7±4.2 deg), participants were instructed to run with anterior trunk lean (ATL; 26.6±5.6 deg) and posterior trunk lean (PTL; 1.8±7.4 deg) within a range in which they felt comfortable when running across even/uneven runways (Fig. 1). Mean trunk flexion angle was calculated as the average sagittal plane trunk posture during the stance phase of running over the level step across all trials and all participants. The order of the ATL and PTL conditions was randomized for each participant. After running on the even, uniform runway, the variable-height force plate at the site of the second contact was visibly lowered by 10 cm and participants ran along the uneven runway (Fig. 1). Practice trials were permitted to allow participants to become familiar with the running velocity and with the desired trunk postures. The participants accomplished 10 valid trials per condition in which each force plate was fully struck with a single foot, whilst the first force plate was set to be hit by the left foot.

Fig. 1.

Schematic illustration of drop negotiation during running with altered trunk postures. The figure depicts the definition of the trunk angle (θtrunk) and the spring–mass model properties as used in this study. A visible drop was created after the approach step () by lowering the second force plate (step down; step↓) by 10 cm. STL, self-selected trunk lean; ATL, anterior trunk lean; PTL, posterior trunk lean; CoM, centre of mass; kleg, leg stiffness; lTD, leg length at touchdown; θleg,TD, leg angle at touchdown.

Fig. 1.

Schematic illustration of drop negotiation during running with altered trunk postures. The figure depicts the definition of the trunk angle (θtrunk) and the spring–mass model properties as used in this study. A visible drop was created after the approach step () by lowering the second force plate (step down; step↓) by 10 cm. STL, self-selected trunk lean; ATL, anterior trunk lean; PTL, posterior trunk lean; CoM, centre of mass; kleg, leg stiffness; lTD, leg length at touchdown; θleg,TD, leg angle at touchdown.

Parameters of interest

The ensemble average of the following parameters in the sagittal plane during the stance phase across the level, unperturbed step (), approach step () and step down (step↓) were determined: (1) vertical position of the CoM relative to the ground, determined by the body segmental analysis method relative to the laboratory coordinate system (De Leva, 1996; Gard et al., 2004); (2) leg length, defined as the length between the hip and the ball of the foot marker (Fig. 1); (3) leg angle, the angle between the leg and the ground (Fig. 1), calculated with respect to the negative x-axis; (4) knee and ankle joint angles; (5) peak vertical GRF (Fz); (6) maximum leg compression (Δlmax) by subtracting the minimum leg length between TD and toe-off (TO) from lTD; (7) leg stiffness (kleg) as the ratio between the Fz and Δlmax; (8) net knee and ankle joint moments calculated using a rigid linked segment model, anthropomorphic data and inverse dynamics analysis (Zatsiorsky, 1983); and (9) knee (kknee) and ankle (kankle) joint stiffness calculated from the ratio of the change in net muscle moment to joint angular displacement between TD and the instant when the joints reached maximal flexion. A vertical GRF threshold of 0.03× body weight was used to determine the instants of TD and TO at each contact (Aminiaghdam et al., 2017b). Leg length, Δlmax and vertical position of the CoM were normalized to the distance between the greater trochanter marker and the lateral malleoli marker at the instant of TD (l0). Fz was normalized to the participant's body weight (N/BW). Each joint moment (N m kg−1) and joint stiffness (N m deg−1 kg−1) were normalized to the subject's body mass.

Data processing

For data analysis, we chose all trials that were distributed in a narrow range of each participant's running steady-state velocity (3.5 m s−1). From the calculated mean horizontal velocity of the L5 marker for each of two force plates, we discarded those trials that contained two values differing by more than 5%. Kinetic and kinematic data of all successful trials were analysed using custom written Matlab (Mathworks Inc., Natick, MA, USA) code. The raw coordinate data were filtered using a fourth-order low-pass, zero-lag Butterworth filter with 12 Hz cutoff frequency (Aminiaghdam et al., 2017b). Following confirmation that the data were normally distributed, a two-way repeated measurements ANOVA was used to examine interactions between posture (STL, ATL and PTL) and step (, and step↓) on the parameters of interest (outlined earlier) in SPSS (v. 21.0, IBM® Co.). In the case of a significant interaction, simple main effects were used to compare between-posture changes across each step as well as between-step changes for each individual posture using one-way ANOVA and post hoc comparisons with Bonferroni adjustments for multiple comparisons. In the case of a non-significant interaction, the main effects of posture and step were evaluated on each dependent variable of interest. Where Mauchly's test indicated a violation of sphericity, P-values and degrees of freedom were corrected using the Greenhouse–Geisser correction factor. Results are expressed as means±s.d. over all participants and parameters. The statistical significance level of all tests was set to P=0.05.

The data analysed included 720 trials with a total of 2160 step cycles. Participants were successful in maintaining their stability (no falls) on every trial while running across the level and uneven track. Table 1 summarizes the kinematic and kinetics parameters of uneven running with altered sagittal plane trunk orientations along with the main effects of step and posture and their interaction effects. Fig. 2 illustrates ensemble-averaged global and local sagittal plane kinematic waveforms, and Fig. 3 illustrates ensemble-averaged joint moment and vertical GRF waveform of running with three trunk orientations across , and step↓.

Fig. 2.

Global and local sagittal plane kinematic waveforms. Shown are ensemble-averaged trunk angular displacement (A–C), normalized body CoM (D–F), normalized leg length (G–I), knee angular displacement (J–L) and ankle angular displacement (M–O) for STL, ATL and PTL running conditions across the stance phase of (control step), and step↓ (N=12). The contact time was normalized to 100%. The grey shaded area represents the corresponding s.d. for the STL condition.

Fig. 2.

Global and local sagittal plane kinematic waveforms. Shown are ensemble-averaged trunk angular displacement (A–C), normalized body CoM (D–F), normalized leg length (G–I), knee angular displacement (J–L) and ankle angular displacement (M–O) for STL, ATL and PTL running conditions across the stance phase of (control step), and step↓ (N=12). The contact time was normalized to 100%. The grey shaded area represents the corresponding s.d. for the STL condition.

Fig. 3.

Joint moment and vertical ground reaction force (GRF) waveforms. Shown are ensemble-averaged normalized knee extensor moments (A–C), ankle extensor moments (D–F) and vertical GRF (Fz) normalized to body weight (BW) (G–I) for STL, ATL and PTL running conditions across the stance phase of , and step↓ (N=12). The contact time was normalized to 100%. The grey shaded area represents the corresponding s.d. for the STL condition.

Fig. 3.

Joint moment and vertical ground reaction force (GRF) waveforms. Shown are ensemble-averaged normalized knee extensor moments (A–C), ankle extensor moments (D–F) and vertical GRF (Fz) normalized to body weight (BW) (G–I) for STL, ATL and PTL running conditions across the stance phase of , and step↓ (N=12). The contact time was normalized to 100%. The grey shaded area represents the corresponding s.d. for the STL condition.

Table 1.

Kinematics and kinetics of uneven running with altered sagittal plane trunk orientations

Kinematics and kinetics of uneven running with altered sagittal plane trunk orientations
Kinematics and kinetics of uneven running with altered sagittal plane trunk orientations

Two-way repeated measures ANOVA indicated step-specific effects of posture on θankle,TD, ΔLmax, θankle,TO, θleg,TO and tcontact (Table 1). Post hoc comparisons (simple main effects) revealed significant changes in some kinematic parameters during the step down (step↓) as compared with those during the level step () and the approach step (), with no between-posture changes across each step (Table 1). In step↓, the ankle angle demonstrated substantial increases in plantarflexion at TD for all running conditions. Likewise, ΔLmax in step↓ significantly increased by ∼35%, ∼70% and 75% during ATL, STL and PTL postures, respectively, compared with those of and (Table 1). The ankle angle at TO of the step↓ increased by ∼45% and ∼60% during STL and ATL conditions, respectively, compared with that of (Table 1). Furthermore, the leg became more vertical across all running conditions, as demonstrated by significant decreases in leg angle (∼4.5 deg during the STL and ATL conditions, and ∼7 deg during the PTL condition), presumably in preparation for accelerating the body CoM for surmounting the drop. The contact time significantly decreased in step↓ compared with , irrespective to the trunk orientation (Table 1).

Analyses considering the main effect of posture revealed significant alterations in running mechanics due to changes in the sagittal trunk lean angle. Looking at the global leg parameters, kleg was significantly reduced during the PTL condition by 15% compared with the STL condition, and by ∼12% compared with the ATL condition (Fig. 4A), whilst Fz remained unchanged. However, kleg did not significantly differ during ATL running compared with the STL condition (Fig. 4A). In addition, PTL running was characterized by a slightly more vertical leg at TD as compared with that in the ATL condition (Fig. 4D). As for the local leg parameters, the knee flexion angle at TD was slightly greater in both PTL and ATL conditions when compared with the STL condition (Fig. 4E). The peak knee extension moment (Mknee) was significantly lower by ∼7% during the ATL condition, and was significantly higher by ∼15% during the PTL condition when compared with the STL condition (Fig. 4C). As compared with the STL condition, kknee was significantly greater by ∼13% during the PTL condition, whereas it remained unchanged during the ATL condition (Fig. 4B). At TO, the knee flexion angle was greater by ∼11% during the PTL condition when compared with the STL condition (Fig. 4F).

Fig. 4.

Main effects of posture. The main effects of posture (mean±s.d.) on variables for which two-way repeated-measurement ANOVA revealed no step×posture interaction (N=12). Leg (A) and knee (B) stiffness (k), peak knee extension moment (Mknee; C), leg (D) and knee (E) angle (θ) at touchdown, and knee angle at toe-off (F). Significant differences from STL and ATL are indicated by * and **, respectively (P<0.05; Bonferroni post hoc test). Error bars indicate s.d.

Fig. 4.

Main effects of posture. The main effects of posture (mean±s.d.) on variables for which two-way repeated-measurement ANOVA revealed no step×posture interaction (N=12). Leg (A) and knee (B) stiffness (k), peak knee extension moment (Mknee; C), leg (D) and knee (E) angle (θ) at touchdown, and knee angle at toe-off (F). Significant differences from STL and ATL are indicated by * and **, respectively (P<0.05; Bonferroni post hoc test). Error bars indicate s.d.

Furthermore, analyses considering the main effect of step revealed significant changes in kinetic and kinematic parameters of running during and step↓ when compared with . In , runners reduced their kleg by ∼24% (Fig. 5G), resulting in a ∼13% decrease in the magnitude of Fz (Fig. 5H). At the ankle level, kankle (Fig. 5I) and Mankle (Fig. 5J) were reduced by ∼24% and 12%, respectively, in . At the end of , runners lowered their CoM by ∼6% (Fig. 5B) using an increased knee flexion (Fig. 5F). This was associated with a ∼2% significant decrease in leg length (Fig. 5D). At TD of step↓, landing involved an increased leg length (Fig. 5C) and a more extended knee joint angle (Fig. 5E), leading to a ∼5% higher vertical position of the CoM (Fig. 5A). As compared with values for , Mankle increased by ∼40% (Fig. 5J), while kankle demonstrated a ∼36% decrease (Fig. 5I) in step↓. kleg was reduced (Fig. 5G), despite a ∼27% increase in Fz (Fig. 5H); however, the difference did not reach significance. This was mainly due to a much greater increase in ΔLmax (Table 1) across all running conditions. At the end of the stance phase, the vertical position of the CoM was enhanced by ∼4% (Fig. 5B), possibly in preparation to surmount the drop.

Fig. 5.

Main effects of step. The main effects of step (mean±s.d.) on variables for which two-way repeated-measurement ANOVA revealed no step×posture interaction (N=12). CoM at touchdown (A) and toe-off (B) normalized to leg length (l0), normalized leg length at touchdown (C) and toe-off (D), knee angle at touchdown (E) and toe-off (F), leg stiffness normalized to body weight (G), peak vertical GRF normalized to body weight (H), ankle stiffness (I), peak normalized ankle extension moment (J) and leg angle at touchdown (K). Significant differences from and are indicated by * and **, respectively (P<0.05; Bonferroni post hoc test). Error bars indicate s.d.

Fig. 5.

Main effects of step. The main effects of step (mean±s.d.) on variables for which two-way repeated-measurement ANOVA revealed no step×posture interaction (N=12). CoM at touchdown (A) and toe-off (B) normalized to leg length (l0), normalized leg length at touchdown (C) and toe-off (D), knee angle at touchdown (E) and toe-off (F), leg stiffness normalized to body weight (G), peak vertical GRF normalized to body weight (H), ankle stiffness (I), peak normalized ankle extension moment (J) and leg angle at touchdown (K). Significant differences from and are indicated by * and **, respectively (P<0.05; Bonferroni post hoc test). Error bars indicate s.d.

An overall objective of our study was providing insight into how changing the sagittal plane trunk orientation influences running mechanics against an expected substrate change. It was postulated that runners would exhibit strategies in global and local running mechanics that are task specific, facilitated by the anticipatory nature of changes to posture and substrate height. This was confirmed by the findings of the present study, as running with a posterior versus anterior trunk lean was associated with greater compensatory kinematic adjustments, resulting in augmented changes in the magnitude of lower limb local and global kinetic parameters. However, runners exhibited step-to-step modulations in local and global running mechanics to account for expected change in the substrate height. Interestingly, the control of body motion for the accommodation of expected substrate height changes appears not to be affected by trunk orientation in running humans.

Posture-related adjustments in running mechanics

When considering the main effect of posture (averaging across the steps), we observed posture-specific adjustments in locomotor behaviour. These results are consistent with those of our previous studies analysing trunk-flexed walking (with the same experimental setup), which revealed that leaning the trunk forward up to 30 deg does not lead to between-step changes in the kinematics (Aminiaghdam et al., 2017a) and kinetics of the braking phase of stance (Aminiaghdam and Rode, 2017) in able-bodied gait. Running with ATL is characterized with a slightly more crouched configuration than that of the STL condition, induced possibly by a forward shift of the CoM. Further, along with an increased knee flexion at TD (Fig. 4E), the knee moment reduced significantly by ∼7% (Fig. 4C) when compared with the STL condition. This finding is in agree­ment with a previous study by Teng and Powers (2014) showing a ∼7% decrease in knee extensor moment at the time of peak patellofemoral joint stress during level running when increasing sagittal trunk lean by ∼7 deg from a self-selected trunk posture (Teng and Powers, 2014). By contrast, the study by dos Santos et al. (2019) revealed no significant decrease in the peak knee extensor moment (1.55±0.32 versus 1.62±0.36 N m kg−1) when runners increase the trunk inclination angle (by ∼9 deg) from their self-selected trunk posture. Moreover, an anterior lean of the trunk does not induce significant changes in joint stiffness across the knee and ankle joints and thus the leg stiffness tends to remain virtually the same as for STL running (Fig. 4A, Table 1). Likewise, the ankle moment appears not to be influenced by alterations in trunk orientation. These results match those observed in an earlier study (Teng and Powers, 2014) that an increase of the sagittal plane trunk inclination redistributes the energetics across the hip and knee without imposing biomechanical demands on the ankle plantar-flexor muscles during level running.

In contrast, a posterior lean of the trunk yields increases in knee flexion at TD (3.4 deg) and TO (2.2 deg) as well as in the knee moment (∼15%) when compared with STL running (Fig. 4). An increased extension moment implies a higher muscle activation of the knee extensors, leading to a higher kknee (Fig. 4B). Teng and Powers (2014) reported a ∼5% increase in the peak knee extensor moment when running with a more extended sagittal trunk lean (Teng and Powers, 2014). A lower extent of the decrease in the extension moment compared with that of our study can be explained by the fact that the degree of posterior lean in our study was profoundly greater. This, in turn, induces more compensatory kinematic adaptations in leg posture. In another similar study, Teng and Powers (2016) demonstrated that weak hip-extensor muscles prompt the adoption of a more upright trunk posture during running (Teng and Powers, 2016). Such postural adaptation appears to reduce the demand on the hip extensors, but at the cost of imposing higher demands on the knee extensors. Moving the trunk backwards shifts the line of action of the GRF away from the axis of rotation at the knee joint, and therefore the larger moment arms result in greater net external moments. This, in turn, increases the angular displacement at the knee joint, and subsequently the compression of the spring-leg. Runners demonstrate adaptations in local stiffness by decreasing kankle by ∼10% and increasing kknee by ∼13% (Fig. 4B). The local regulation of the joint stiffness coupled with the combination of increased leg compression and decreased Fz yields a reduction in kleg for the PTL condition (Fig. 4B). Consequently, the sagittal plane trunk posture has an influence on biomechanical demands (global and local) of the lower limb during running, but this is independent of changes in substrate height.

Step-related adjustments in running mechanics

It appears that the anticipation of substrate height change governs the control scheme of running. Runners exhibit substantial adjustments in body mechanics, namely feed-forward strategies to accommodate substrate height changes while adopting various trunk orientations. This includes changes in the spring-mass properties (i.e. lTD, θleg,TD and kleg) and local mechanics across the knee and ankle joints. The properties of the spring-mass model are modulated with an attenuation bias in the preceding contact when coping with drops and downward steps, and with an augmentation bias in the lowered level (Müller et al., 2012a). These adjustments are utilized to smooth the CoM kinematics (Blickhan et al., 2007). In our study, analysis of the main effect of step revealed the global and local adjustments in the locomotor behaviour specific to the step nature during uneven running. In the approach step, runners demonstrated a ∼24% reduction in leg stiffness in preparation to step down when compared with the control step (Fig. 5G). Likewise, ankle stiffness was reduced by ∼24% (Fig. 5I). As a result, global leg stiffness appears to be mainly adjusted locally at the level of the ankle joint. These modulations in local and global stiffness gave rise to changes in the magnitude of the leg force, which decreased by ∼13% (Fig. 5H). Given the joint stiffness also depends on the activation level of the muscles acting on the joint (Müller et al., 2010, 2012b), an attenuation of global leg stiffness may be due to decreased muscle activation, and a more flexed leg configuration. Here, runners demonstrated a 12% reduction in the ankle moment (Fig. 5J). However, as they approached the end of the stance, knee flexion (Figs 2K and 5F) increased so that the leg length (Figs 2H and 5D) decreased by ∼2% and subsequently the CoM was lowered by ∼6% at TO (Figs 2E and 5B). As for obstacle negotiation, kleg decreased in the penultimate step and then increased in the final step prior to a 0.65 m-high barrier in an attempt to enhance the movement of the CoM while leaping the obstacle (Mauroy et al., 2013). It seems that the runners resort to spring-leg stiffness attenuation at the obstacle or upward step, but also in the preparatory step(s) prior to a drop or downward step (Müller et al., 2012a). The findings observed in our study mirror those of previous studies that have reported decreases in the stiffness of the ankle joint proportional to the elevation of the next step during running (Grimmer et al., 2008; Mauroy et al., 2013, 2014; Müller et al., 2012a, 2010). Müller et al. (2012a, 2010) suggested that the modulation of the leg stiffness is actively achieved.

In the step down (drop), the perturbed leg lands in a more vertical orientation with more extended knee and ankle joint positions, resulting in an elongated leg (Figs 2I and 5C). Similarly, Müller et al. (2012a) demonstrated increases in leg angle and leg length when accommodating expected/unexpected drops and expected steps during running. However, they noted that the negotiation of the unexpected versus expected perturbations is associated with greater leg adjustments at TD, due possibly to the lack of feed-forward control of motion prior to the perturbation. It seems that the modulation of the TD leg angle plays a crucial role in the dynamics of stance following the perturbation. Such mechanical behaviour appears to also be exploited by small birds as Daley and Biewener (2006) noted that the extended limb posture at TD of an unexpected drop accounts for 80% of their stance limb loading. Furthermore, given a relatively straight architecture of the human leg, the longer leg at TD is achieved mainly by the ankle plantar flexion (Müller and Blickhan, 2010; Müller et al., 2012a) and slightly by the knee extension. The coupling of these adjustments in the perturbed leg with a posterior motion of the trunk across all running conditions (Fig. 2C) slightly enhances the vertical position of the CoM at TD (Fig. 2F). Seethapathi and Srinivasan (2019) showed that runners adjust leg force over the whole stance and foot placement in a manner to recover the CoM trajectory back to steady state as they run. They further demonstrated the application of these control strategies for locomotion in simulation, showing the robustness of a simple biped against larger discrete perturbations or constant noise-like perturbations. Here, the peak leg force increased by 27% in the perturbed step but decreased by 13% in the preparatory step when compared with the level step. In agreement with Seethapathi and Srinivasan (2019), running humans exhibit adjustments in leg loading and CoM motion in the approach step and the step down (Figs 2D–F and 5A,B) to facilitate the traverse of expected substrate height changes.

In locomotion, the foot segment enables forward propulsion by serving as the link between the ground and the kinetic chain of the lower limb and trunk. Compared with other lower limb moments and forces, the primary contributor to the vertical GRF is the ankle moment (Chen, 2006). The ankle extensor moment exhibits a ∼40% increase (Figs 3F and 5J), due potentially to a higher activation of the plantar-flexor muscles, which stems particularly from an increased pre-activation of the m. gastrocnemius medialis (Müller et al., 2015). Our findings accord with those of previous studies in terms of the modulation of leg stiffness (Müller and Blickhan, 2010; Müller et al., 2012a), which have shown that stepping down off an elevation (e.g. single drop of different elevations) is not associated with changes in leg stiffness. This can be due to the extensive leg compression displayed across all running conditions. However, the duration of the ground contact during the lowered (perturbed) step is significantly shorter than during the level steps, irrespective to the running postures (Table 1). Towards the end of the drop, the CoM rises (by ∼4%) through greater plantar flexion and a more vertical leg orientation (regardless of trunk posture), possibly in preparation to surmount the drop (Figs 2F and 5B). The step down therefore appears to be dominated by CoM height changes, regardless of trunk lean.

When interpreting the findings of this study, several limitations need to be acknowledged. First, the runway utilized in the present study may not fully resemble running outdoors, which represents a wide range of expected/unexpected and/or varied magnitudes of surface perturbations, and thus possibly elicits different locomotor outputs. Second, the sample size of the study might be relatively small to reach a convincing conclusion. Third, our results also do not exclude the possibility of being influenced by using an overly conservative Bonferroni correction for multiple comparisons which preserves type I error of the global null hypothesis. Overall, the step-to-step adjustments in various running mechanical parameters exhibit no dependency on sagittal trunk orientation. It is most likely that the observed adjustments are influenced by feed-forward control as both postural and environmental changes were imposed experimentally in an anticipatory fashion. The reliance of the stability of the sagittal-plane dynamics of running on anticipatory strategies when traversing rough terrains (i.e. slope and height variations) has also been validated by a simulation study (Dhawale et al., 2019). In agreement with a previous study (Qiao and Jindrich, 2012), our findings demonstrate that human locomotion employs task-level strategies to account for either or both external and internal perturbations. Characterizing the mechanical behaviour governed by feed-forward control strategies could help better understanding of sensory-motor mechanisms underlying the stabilization of human perturbed locomotion. By illustrating adjustments in global leg-spring stiffness specific to alterations in posture and ground configurations, our findings further support the notion of dependency of the stability of spring-mass running through the global leg stiffness on local joint elasticity, and on leg geometry at TD (Arampatzis et al., 1999; Farley et al., 1998; Farley and Morgenroth, 1999; Günther and Blickhan, 2002; Mauroy et al., 2014). These findings may improve understanding of the role of posture in human perturbed locomotion relevant for the design and development of exoskeleton or humanoid bipedal robots.

Author contributions

Conceptualization: S.A., R.B.; Methodology: S.A., R.B.; Software: S.A., R.B.; Validation: S.A., R.B., K.K.; Formal analysis: S.A.; Investigation: S.A.; Resources: S.A., R.B., K.K.; Data curation: S.A.; Writing - original draft: S.A.; Writing - review & editing: S.A., R.B., K.K.; Visualization: S.A.; Supervision: S.A., R.B., K.K.; Project administration: S.A., R.B., K.K.

Funding

This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors.

Aminiaghdam
,
S.
and
Rode
,
C.
(
2017
).
Effects of altered sagittal trunk orientation on kinetic pattern in able-bodied walking on uneven ground
.
Biol. Open
6
,
1000
-
1007
.
AminiAghdam
,
S.
and
Blickhan
,
R.
(
2018
).
The effects of an expected twofold perturbation on able-bodied gait: Trunk flexion and uneven ground surface
.
Gait Posture
61
,
431
-
438
.
Aminiaghdam
,
S.
,
Blickhan
,
R.
,
Muller
,
R.
and
Rode
,
C.
(
2017a
).
Posture alteration as a measure to accommodate uneven ground in able-bodied gait
.
PLoS ONE
12
,
e0190135
.
Aminiaghdam
,
S.
,
Rode
,
C.
,
Müller
,
R.
and
Blickhan
,
R.
(
2017b
).
Increasing trunk flexion transforms human leg function into that of birds despite different leg morphology
.
J. Exp. Biol.
220
,
478
-
486
.
Arampatzis
,
A.
,
Brüggemann
,
G.-P.
and
Metzler
,
V.
(
1999
).
The effect of speed on leg stiffness and joint kinetics in human running
.
J. Biomech.
32
,
1349
-
1353
.
Arendse
,
R. E.
,
Noakes
,
T. D.
,
Azevedo
,
L. B.
,
Romanov
,
N.
,
Schwellnus
,
M. P.
and
Fletcher
,
G.
(
2004
).
Reduced eccentric loading of the knee with the pose running method
.
Med. Sci. Sports Exerc.
36
,
272
-
277
.
Blickhan
,
R.
(
1989
).
The spring-mass model for running and hopping
.
J. Biomech.
22
,
1217
-
1227
.
Blickhan
,
R.
,
Seyfarth
,
A.
,
Geyer
,
H.
,
Grimmer
,
S.
,
Wagner
,
H.
and
Günther
,
M.
(
2007
).
Intelligence by mechanics
.
Phil. Trans. R. Soc. A Math. Phys. Eng. Sci.
365
,
199
-
220
.
Chen
,
G.
(
2006
).
Induced acceleration contributions to locomotion dynamics are not physically well defined
.
Gait Posture
23
,
37
-
44
.
Daley
,
M. A.
and
Biewener
,
A. A.
(
2006
).
Running over rough terrain reveals limb control for intrinsic stability
.
Proc. Natl Acad. Sci. USA
103
,
15681
-
15686
.
De Leva
,
P.
(
1996
).
Adjustments to Zatsiorsky-Seluyanov's segment inertia parameters
.
J. Biomech.
29
,
1223
-
1230
.
Dhawale
,
N.
,
Mandre
,
S.
and
Venkadesan
,
M.
(
2019
).
Dynamics and stability of running on rough terrains
.
R. Soc. Open Sci.
6
,
181729
.
dos Santos
,
A. F.
,
Nakagawa
,
T. H.
,
Serrão
,
F. V.
and
Ferber
,
R.
(
2019
).
Patellofemoral joint stress measured across three different running techniques
.
Gait Posture
68
,
37
-
43
.
Ernst
,
M.
,
Götze
,
M.
,
Müller
,
R.
and
Blickhan
,
R.
(
2014
).
Vertical adaptation of the center of mass in human running on uneven ground
.
Hum. Mov. Sci.
38
,
293
-
304
.
Farley
,
C. T.
,
Houdijk
,
H. H.
,
Van Strien
,
C.
and
Louie
,
M.
(
1998
).
Mechanism of leg stiffness adjustment for hopping on surfaces of different stiffnesses
.
J. Appl. Physiol.
85
,
1044
-
1055
.
Farley
,
C. T.
and
Morgenroth
,
D. C.
(
1999
).
Leg stiffness primarily depends on ankle stiffness during human hopping
.
J. Biomech.
32
,
267
-
273
.
Gard
,
S. A.
,
Miff
,
S. C.
and
Kuo
,
A. D.
(
2004
).
Comparison of kinematic and kinetic methods for computing the vertical motion of the body center of mass during walking
.
Hum. Mov. Sci.
22
,
597
-
610
.
Grimmer
,
S.
,
Ernst
,
M.
,
Günther
,
M.
and
Blickhan
,
R.
(
2008
).
Running on uneven ground: leg adjustment to vertical steps and self-stability
.
J. Exp. Biol.
211
,
2989
-
3000
.
Günther
,
M.
and
Blickhan
,
R.
(
2002
).
Joint stiffness of the ankle and the knee in running
.
J. Biomech.
35
,
1459
-
1474
.
Huang
,
Y.
,
Xia
,
H.
,
Chen
,
G.
,
Cheng
,
S.
,
Cheung
,
R. T. H.
and
Shull
,
P. B.
(
2019
).
Foot strike pattern, step rate, and trunk posture combined gait modifications to reduce impact loading during running
.
J. Biomech.
86
,
102
-
109
.
Lee
,
D.-c.
,
Brellenthin
,
A. G.
,
Thompson
,
P. D.
,
Sui
,
X.
,
Lee
,
I.-M.
and
Lavie
,
C. J.
(
2017
).
Running as a key lifestyle medicine for longevity
.
Prog. Cardiovasc. Dis.
60
,
45
-
55
.
Mauroy
,
G.
,
Schepens
,
B.
and
Willems
,
P. A.
(
2013
).
The mechanics of running while approaching and jumping over an obstacle
.
Eur. J. Appl. Physiol.
113
,
1043
-
1057
.
Mauroy
,
G.
,
Schepens
,
B.
and
Willems
,
P.
(
2014
).
Leg stiffness and joint stiffness while running to and jumping over an obstacle
.
J. Biomech.
47
,
526
-
535
.
McMahon
,
T. A.
and
Cheng
,
G. C.
(
1990
).
The mechanics of running: how does stiffness couple with speed?
J. Biomech.
23
,
65
-
78
.
Müller
,
R.
and
Blickhan
,
R.
(
2010
).
Running on uneven ground: leg adjustments to altered ground level
.
Hum. Mov. Sci.
29
,
578
-
589
.
Müller
,
R.
,
Grimmer
,
S.
and
Blickhan
,
R.
(
2010
).
Running on uneven ground: leg adjustments by muscle pre-activation control
.
Hum. Mov. Sci.
29
,
299
-
310
.
Müller
,
R.
,
Ernst
,
M.
and
Blickhan
,
R.
(
2012a
).
Leg adjustments during running across visible and camouflaged incidental changes in ground level
.
J. Exp. Biol.
215
,
3072
-
3079
.
Müller
,
R.
,
Siebert
,
T.
and
Blickhan
,
R.
(
2012b
).
Muscle preactivation control: simulation of ankle joint adjustments at touchdown during running on uneven ground
.
J. Appl. Biomech.
28
,
718
-
725
.
Müller
,
R.
,
Häufle
,
D. F. B.
and
Blickhan
,
R.
(
2015
).
Preparing the leg for ground contact in running: the contribution of feed-forward and visual feedback
.
J. Exp. Biol.
218
,
451
-
457
.
Qiao
,
M.
and
Jindrich
,
D. L.
(
2012
).
Task-level strategies for human sagittal-plane running maneuvers are consistent with robotic control policies
.
PLoS ONE
7
,
e51888
.
Seethapathi
,
N.
and
Srinivasan
,
M.
(
2019
).
Step-to-step variations in human running reveal how humans run without falling
.
ELife
8
,
e38371
.
Seyfarth
,
A.
,
Geyer
,
H.
and
Herr
,
H.
(
2003
).
Swing-leg retraction: a simple control model for stable running
.
J. Exp. Biol.
206
,
2547
-
2555
.
Shih
,
H.-T.
,
Teng
,
H.-L.
,
Gray
,
C.
,
Poggemiller
,
M.
,
Tracy
,
I.
and
Lee
,
S.-P.
(
2019
).
Four weeks of training with simple postural instructions changes trunk posture and foot strike pattern in recreational runners
.
Phys. Ther. Sport
35
,
89
-
96
.
Teng
,
H.-L.
and
Powers
,
C. M.
(
2014
).
Sagittal plane trunk posture influences patellofemoral joint stress during running
.
J. Orthop. Sports Phys. Ther.
44
,
785
-
792
.
Teng
,
H.-L.
and
Powers
,
C. M.
(
2015
).
Influence of trunk posture on lower extremity energetics during running
.
Med. Sci. Sports Exerc.
47
,
625
-
630
.
Teng
,
H.-L.
and
Powers
,
C. M.
(
2016
).
Hip-extensor strength, trunk posture, and use of the knee-extensor muscles during running
.
J. Athl. Train.
51
,
519
-
524
.
Voloshina
,
A. S.
and
Ferris
,
D. P.
(
2015
).
Biomechanics and energetics of running on uneven terrain
.
J. Exp. Biol.
218
,
711
-
719
.
Wu
,
G.
,
Siegler
,
S.
,
Allard
,
P.
,
Kirtley
,
C.
,
Leardini
,
A.
,
Rosenbaum
,
D.
,
Whittle
,
M.
,
D'Lima
,
D. D.
,
Cristofolini
,
L.
and
Witte
,
H.
(
2002
).
ISB recommendation on definitions of joint coordinate system of various joints for the reporting of human joint motion—part I: ankle, hip, and spine
.
J. Biomech.
35
,
543
-
548
.
Wu
,
G.
,
van der Helm
,
F. C.
,
Veeger
,
H. D.
,
Makhsous
,
M.
,
Van Roy
,
P.
,
Anglin
,
C.
,
Nagels
,
J.
,
Karduna
,
A. R.
,
McQuade
,
K.
and
Wang
,
X.
(
2005
).
ISB recommendation on definitions of joint coordinate systems of various joints for the reporting of human joint motion—Part II: shoulder, elbow, wrist and hand
.
J. Biomech.
38
,
981
-
992
.
Zatsiorsky
,
V.
(
1983
).
The mass and inertia characteristics of the main segments of the human body
.
Biomechanics
4B
,
1152
-
1159
.

Competing interests

The authors declare no competing or financial interests.