Load Magnitude and Locomotion Pattern Alter Locomotor System Function in Healthy Young Adult Women

Abstract

Introduction:

During cyclical steady state ambulation, such as walking, variability in stride intervals can indicate the state of the system. In order to define locomotor system function, observed variability in motor patterns, stride regulation and gait complexity must be assessed in the presence of a perturbation. Common perturbations, especially for military populations, are load carriage and an imposed locomotion pattern known as forced marching (FM). We examined the interactive effects of load magnitude and locomotion pattern on motor variability, stride regulation and gait complexity during bipedal ambulation in recruit-aged females.

Methods:

Eleven healthy physically active females (18–30 years) completed 1-min trials of running and FM at three load conditions: no additional weight/bodyweight (BW), an additional 25% of BW (BW + 25%), and an additional 45% of BW (BW + 45%). A goal equivalent manifold (GEM) approach was used to assess motor variability yielding relative variability (RV; ratio of “good” to “bad” variability) and detrended fluctuation analysis (DFA) to determine gait complexity on stride length (SL) and stride time (ST) parameters. DFA was also used on GEM outcomes to calculate stride regulation.

Results:

There was a main effect of load (p = 0.01) on RV; as load increased, RV decreased. There was a main effect of locomotion (p = 0.01), with FM exhibiting greater RV than running. Strides were regulated more tightly and corrected quicker at BW + 45% compared (p < 0.05) to BW. Stride regulation was greater for FM compared to running. There was a main effect of load for gait complexity (p = 0.002); as load increased gait complexity decreased, likewise FM had less (p = 0.02) gait complexity than running.

Discussion:

This study is the first to employ a GEM approach and a complexity analysis to gait tasks under load carriage. Reduction in “good” variability as load increases potentially exposes anatomical structures to repetitive site-specific loading. Furthermore, load carriage magnitudes of BW + 45% potentially destabilize the system making individuals less adaptable to additional perturbations. This is further evidenced by the decrease in gait complexity, which all participants demonstrated values similarly observed in neurologically impaired populations during the BW + 45% load condition.

Introduction

Bipedal ambulation requires the complex integration of multisensory information (optical, somatic and vestibular) that is used to coordinate actions within specific environments in order to achieve goal-directed movement (Alexander, 1992; Warren et al., 2001; Bent et al., 2004a; Pandy and Andriacchi, 2010; Matthis et al., 2017). Perceptions of continuously obtained multisensory information yield opportunities to act (affordances) resulting in a perception-action coupling, with a specific movements success predicated on the modulation (tuning and weighting) of the afferent signals that provide (or fail to) appropriate affordances for the task (Gibson, 1966, 1979; Hollands and Marple-Horvat, 1996; Warren et al., 2001; Bent et al., 2004a,b; Rossignol et al., 2006; Peters et al., 2017). In conjunction with sensory “reafference,” feedforward mechanisms stimulate coordinative structures or muscle synergies that produce a desired movement that achieves a locomotion task goal (Kim et al., 2011; Bizzi and Cheung, 2013; Minassian et al., 2017). Collectively, the reciprocal cooperation of feedback (afferent) and feedforward (efferent) subcomponents executing a locomotion task is known as the locomotor system. The function of the locomotor system reflects the emergent properties of the organization of the degrees of freedom during locomotor tasks, with specific macroscopic pattern of organization being influenced by the confluence of cost functions (i.e., metabolic efficiency and energy dampening), task, organism (including feedback and feedforward processes) and environmental constraints (i.e., gravity, uneven terrain) (Turvey, 1990; Newell and Vaillancourt, 2001; Davids et al., 2003; Sánchez et al., 2016; Caballero et al., 2019; Shafizadeh et al., 2019). Optimal locomotor system function is represented by biomechanical output that is both stable and adaptive to perturbation (Davids et al., 2003; West and Scafetta, 2003; Cusumano and Dingwell, 2013; Seifert et al., 2013). A common perturbation to bipedal locomotion, especially in military populations, is load carriage, especially “combat load” magnitudes of 20–30 kg (Taylor et al., 2016; Krajewski et al., 2020). How the locomotor system accommodates increasing load magnitudes to successfully execute locomotion task goals still remain unclear (LaFiandra et al., 2003; Attwells et al., 2006; Walsh et al., 2018). Thus, measuring the responses of biomechanical variables to the perturbation of additional loading during locomotory tasks provides valuable insight to the global functional state of the locomotor system.

Variability in the observed movement patterns (motor variability) represents the observed variation in a movement solution, when attempting to accomplish the same goal/task, such as using different segment coupling patterns to perform a step (Bernstein, 1967; Latash et al., 2002, 2010; Latash, 2016). The multitude of joints and muscles in the lower extremity lead to a large number of degrees of freedom that lends itself to equifinality; infinite number of movement solutions to accomplish the same task (Bernstein, 1967; Gelfand and Latash, 1998; Latash et al., 2010). A goal equivalent manifold (GEM; equifinality technique) approach seeks to quantify the “good” (plotted tangential to the GEM [δ_T]) versus “bad” (plotted perpendicular to the GEM [δ_P]) motor variability to further discriminate optimal performance (known as relative variability [the ratio of “good” motor variability to “bad” motor variability]) (Dingwell et al., 2010; Cusumano and Dingwell, 2013; Sedighi and Nussbaum, 2019). Recent theories have demonstrated that motor variability not only leverages equifinality, making the system more adaptable and stable to perturbation (i.e., overcoming varying terrain or recovering from a slip/trip) (Cusumano and Dingwell, 2013; Dingwell et al., 2017), it also has other cost function benefits (Gates and Dingwell, 2008). Specifically, by capitalizing on a larger workspace (greater relative variability) of movement patterns to perform steady-state (constant locomotion velocity) behaviors, energy can be dispersed through more supportive, anatomical structures, whereas limited motor variability (lower relative variability) may lead to site-specific mechanical overloading (cumulative mechanical stress) that can result in musculoskeletal injury (MSI) (Baida et al., 2018; Nordin and Dufek, 2019). Likewise, motor variability can distribute positive mechanical workloads across a greater number of muscle fibers improving metabolic efficiency by reducing the fatigue of a specific subset of muscle fibers (Gates and Dingwell, 2008).

Regulation of cyclical movements during steady-state behavior such as corrections of stride-to-stride fluctuations further elucidates the state of the locomotor system (Cusumano and Dingwell, 2013). Stride regulation is determined by statistical persistence assessment [alpha coefficients (Dingwell et al., 2017)] of deviations tangential (good variability) [δ_T] and orthogonal (bad variability) [δ_P] to the goal manifold (Dingwell et al., 2017). A seminal investigation by Dingwell et al. (2017) demonstrated that elderly individuals classified as low risk fallers and healthy young adults had the same amount of relative variability (ratio of “good” to “bad” motor variability) and used similar stride regulation strategies indicative of a minimal intervention principle (δ_Tα > 1; δ_Pα < 0.5) (Dingwell et al., 2017). Furthermore, it was suggested that changes in stride-regulation strategy to an absolute position control (POS) model [δ_Tα and δ_Pα < 0.5] (Cusumano and Dingwell, 2013) may be the determinant of fall risk (Dingwell et al., 2017). The latter finding was determined with computational modeling (based on a minimum intervention principle) and is still theoretical at this point (Dingwell et al., 2017), but the use of a POS regulation strategy may indicate perception heavily tuned on their exact position, neglecting/overpowering other important information which will impact affordances perception. In the case of military personnel, especially infantry, load carriage is only one perturbation that must be overcome in addition to uneven terrain, enemy threats and decision making. Thus, the quantification of regulation strategy of the system used for stride to stride fluctuations acts as an indirect assessment of the perception-action loop function namely: (i) the ability to (re)calibrate information-action in a dynamic environment, (ii) (re)weighting the relative importance of information sources as they become available, and (iii) modulate based on the relative importance in relation to the successful maintenance of a functionally useful action-response (Cusumano and Dingwell, 2013; Roerdink et al., 2019).

Components of the locomotor system operate/evolve over different time scales and configure in a heterarchical organization when functioning optimally (Bak et al., 1987; Turvey, 1990; Bak and Paczuski, 1995; Newell and Vaillancourt, 2001; Davids et al., 2003; Van Orden et al., 2003). A heterarchical organization of a dynamical system is considered to be complex (interaction of many independent subcomponents that yield an emergent behavior) and a perturbation of one subcomponent is less likely to affect the system globally (West and Shlesinger, 1989; Bak and Paczuski, 1995; Marks-Tarlow, 1999; Torre et al., 2007; Torre and Balasubramaniam, 2009). Thus quantification of system complexity indicates the state of dynamical system health (Iyengar et al., 1996; Gisiger, 2001; Goldberger et al., 2005; Hausdorff, 2007; Van Orden et al., 2009; Nourrit-Lucas et al., 2015; Torre et al., 2019) through non-linear signal processing techniques to determine the fractal structure of a time-series, which exhibits self-similarity at different time scales (Stadnitski, 2012). These fractals display long-range correlations, or learning behavior of current iterations from previous iterations (Hausdorff et al., 1995). Time-series structures of gait dynamics (stride length [S_L] and stride time [S_T]) that yield long range correlations (pink noise) have been linked to healthy functioning adults (Hausdorff et al., 1997, 1999; Hausdorff, 2007; Delignières and Torre, 2009; Ducharme et al., 2018); however, very strong long-range correlations exhibit over-regularity (brown noise) (Gisiger, 2001). Signals that are completely stochastic (white noise) demonstrate no correlation between strides and have been associated with individuals suffering central neurological impairment (Hausdorff et al., 1997; Hausdorff, 2009; Moon et al., 2016). Moreover, white noise has also been observed when imposing a frequency on cyclical steady-state behavior (Terrier et al., 2005; Terrier and Dériaz, 2012; Hunt et al., 2014; Ducharme et al., 2018; Roerdink et al., 2019). Interestingly, warfighters are encouraged to utilize a walking pattern during a velocity that exceeds the gait transition velocity (GTV), colloquially known as forced marching (FM) that is an unnatural (imposed frequency) gait. Little is known how load magnitude, especially military relevant loads (20–60 kg) (Taylor et al., 2016), and this imposed locomotion affect gait complexity in healthy individuals.

To date a load magnitude perturbation is evidenced only in terms of increased mechanical [greater ground reaction forces (GRF) (Birrell et al., 2007; Seay et al., 2014b) and joint kinetics (Knapik et al., 2004; Seay et al., 2014a,b; Liew et al., 2016; Willy et al., 2016, 2019; Lenton et al., 2019; Loverro et al., 2019; Wills et al., 2019; Krajewski et al., 2020)] and physiological [increased heart rate and ratings of perceived exertion (Simpson et al., 2010, 2011, 2017; Huang and Kuo, 2014)] demands compared to unloaded bipedal ambulation. The majority of these studies consisted of male dominated samples, leaving females underrepresented in load carriage research (Loverro et al., 2019). In addition, females are at twice the risk of MSI (Molloy et al., 2020), with a high incidence (∼78%) of MSI observed during basic combat training (recruits), the majority (30–64%) of those MSI suffered during load carriage conditioning in basic training (Jensen et al., 2019; Lovalekar et al., 2020) suggesting that individuals with little to no experience with load carriage tasks are of greater interest. However, there is a paucity of information regarding the effects of load carriage on motor control (LaFiandra et al., 2003; Attwells et al., 2006; Walsh et al., 2018). Importantly, previous work has focused on average behavior of spatiotemporal gait parameters (LaFiandra et al., 2003; Attwells et al., 2006) and have yet to elucidate key features of a healthy locomotor system such as motor variability, stride to stride regulation and the complexity of the system. Therefore, the purpose of this investigation was to determine the interactive effects of load magnitude and locomotion pattern on motor variability, stride regulation and gait complexity during bipedal ambulation in recruit aged females. It is hypothesized, based on an affordance-based control theory (Davids et al., 2003; Mukherjee and Yentes, 2018) that as load increases and the use of an unnatural (imposed) locomotion (FM) will constrain the locomotor system decreasing the number of affordances available which will be reflected by the reduction in relative variability (the ratio of “good” motor variability to “bad” motor variability). Likewise, increases in load and utilization of FM will lead to stricter regulation strategies. Lastly, individuals gait complexity will decrease as load increases and during the execution of the FM locomotion pattern.

Materials and Methods

Ethics Statement

All participants were read and signed informed consent that has been approved by the Institutional Review Board (IRB) of The University of Pittsburgh. They were notified of potential risks and benefits associated with participation in the study.

Subjects and Protocol

Eleven healthy, recreationally active young adult females (see Table 1 for participant characteristics) participated in this study. Recreationally active was defined as engaging in moderate physical activity a minimum of two times a week for at least 30 min, similar to comparable recruits. Moreover, women novice to load carriage and forced marching were chosen to represent a female recruit population, replicating initial exposure to load carriage tasks. Subjects were screened to exclude individuals who reported spine and lower extremity musculoskeletal injury, neurological disorder or pregnant.

The procedures for this investigation have been previously described in detail (Krajewski et al., 2020). Briefly, participants ran (RN) and forced marched (FM) on a instrumented split-belt treadmill (Bertec Corporation, Columbus, OH, United States) for 1 min at three different loaded conditions: Bodyweight (BW), plus an additional 25% of BW (BW + 25%), and plus an additional 45% of BW (BW + 45%) [which represents 20–30 kg “combat” loads in average young adult females (Taylor et al., 2016)] at 10% above their GTV (BW: 2.08 ± 0.25 m/s; BW + 25%: 2.02 ± 0.22 m/s; BW + 45%: 1.93 ± 0.23 m/s). All participants wore provided combat boots (Speed 3.0 Boot, 5.11 Tactical, Irvine, CA, United States) to control for influences of footwear on kinematics (Telfer et al., 2017) and loaded conditions were executed with an anterior-posterior loaded weight vest (Short Plus Style Vest, MIR, United States) to control for effects of center of mass (COM) location (Seay et al., 2014a; Loverro et al., 2019). All trials were randomized by load condition and then by locomotion pattern. Participants were given up to 5 min between each trial to control for effects of fatigue. During RN trials, participants were instructed to move “naturally” or how they felt most comfortable to maintain treadmill velocity. For FM trials, participants were instructed to maintain a walking gait regardless of the treadmill velocity. Each trial yielded ∼130 strides (120–180) dependent on the locomotion pattern and velocity. Prior to familiarization and data collection, participants filled out an activity questionnaire and body composition was assessed with dual energy X-ray absorptiometry (DXA) [Lunar iDXA, General Electric, Boston, MA, United States].

Data Collection and Processing

Three retro-reflective markers were placed on each boot (calcaneus, 1st and 5th metatarsophalangeal [MTP] joints) [see Figure 1 for subject experimental set up]. Kinematic data was collected via 12 infrared cameras (Vicon Motion Systems Ltd., Oxford, United Kingdom) sampling at 100 Hz. Kinetic data was collected via an instrumented split-belt treadmill sampling at 1000 Hz that was synchronized with the motion analysis system. Using the Vicon Nexus^® 2.0 software (Vicon Motion Systems Ltd., Oxford, United Kingdom), a custom labeling template was created for the marker configuration used in the study. Once all static and motion trials were reconstructed, the labeling template was used to auto label the static trials captured for each load condition (BW, BW + 25% and BW + 45%) which were then used to auto label their respective motion trials (RN and FM). Gap filling methods in Nexus 2.0 were used to correct any breaks in trajectory data due to marker occlusion. Data was then exported, and post processed in Visual 3D (C-Motion Inc., Germantown, MD, United States). Further analysis [GEM decomposition and Detrended Fluctuation Analysis (Atlas Collaboration et al., 2014)] was conducted with custom Matlab^TM 2019a (Mathworks, Inc., Natick, MA, United States) scripts. Kinematic and kinetic data were filtered with a second order Butterworth low-pass filter (cut-off frequencies of 6 Hz and 40 Hz for the kinematic and kinetic signals, respectively). Heel strike was defined as the time when vertical component of the ground reaction force exceeded a 50N threshold.

FIGURE 1

The following variables were calculated: Stride length (S_L) was computed as the distance covered from heel strike to ipsilateral heel strike; Stride time (S_T) was computed as the time elapsed from heel strike to ipsilateral heel strike; Stride speed (S_S) was computed as the quotient of S_L/S_T; Velocity (v) was computed as the average S_S over all n strides of a time-series. Average values (Means), standard deviations (SD) and DFA scaling exponents (α-value) were calculated for S_L, S_T and S_S across all trials.

Goal Equivalent Manifold Decomposition

Methods utilized for GEM decomposition have been described in detail by Dingwell et al. (2010). However, to further elaborate the process: firstly, S_L and S_T time-series for each trial was normalized to unit variance [dividing by its own standard deviation]. A specific operating point was computed for S_T as Eq. (1):

Where <■> represents the average across all n strides of the time series. The specific operating point for S_L was computed as Eq. (2):

Here v represents the velocity of the treadmill for that specific trial. The new centered operating point was then computed as Eqs. (3) and (4):

and

Lastly, deviations tangential to the goal manifold were represented as δ_T and deviations perpendicular to the goal manifold were represented as δ_P. These deviations were calculated with a linear coordinate transformation as Eq. (5):

Where, the σ of δ_T and δ_P were determined for each load and locomotion condition. Relative variability was calculated as the ratio between σδ_T/σδ_P. Therefore, a relative variability magnitude of 1 represents equal amounts of “good” versus “bad” variability; <1 represents more “bad” variability; and >1 represents more “good” variability. Additionally, DFA scaling exponents (α) were computed for δ_T and δ_P. Scaling exponents for δ_T and δ_P are interpreted as follows: α < 0.5 represents anti-persistence (alteration in one direction more likely followed by an alteration in opposite direction); α > 0.5 represents statistical persistence (alteration in one direction more likely followed by an alteration in same direction); and α = 0.5 represents uncorrelated (alteration in one direction has same likelihood of being followed by alteration in either direction) (Dingwell et al., 2010, 2017).

Complexity Analysis

Complexity analysis was executed utilizing fractal methods, specifically DFA (Peng et al., 1993; Hausdorff et al., 1995; Delignières et al., 2006; Stadnitski, 2012) on S_L and S_T gait variables (∼130 consecutive strides). Refer to the aforementioned references for greater detail but briefly: DFA creates a one-dimensional signal x(i), i = 1,…, L, where x is the initial signal of size L, and an integrated signal (x) is calculated according to the Eq. (6):

where B_avg is the mean value of the signal (B = signal value at specific time point). The unified time series Y is then divided into segments (boxes that don’t overlap) of length l, and the linear approximation Y_l is then obtained through a least-squares fit of each segment separately (trend of each section).

The mean fluctuation (root mean square) of the incorporated and detrended time-series is computed using Eq. (7):

The aforementioned calculations are repeated for a range of l. The goal of this analysis is to identify the relation between F(l) and the size of segment l because this relationship serves as an indicator of a scaling phenomenon. In general, F(l) increases with increases in the range of segment l. A double plot logarithmic graph (log (F(l) vs log_l) is then formed, and this graph is used to acquire the scaling exponent (α). A linear dependency implies the existence of self-fluctuations, and F(l) which is the slope of line outlines the scaling α exponent, which increases with l based on a power law, as detailed in Eq. (6):

DFA ultimately yields a scaling exponent (α) which represent the correlational structure of the signal. White noise (uncorrelated or completely stochastic) is represented as α = 0.5; Pink noise (positive long-range correlations) is represented as α = 1.0; Brown noise (persistent long-range correlations or too much regularity) is represented as α = 1.5 (Peng et al., 1995). Classifications based upon a range of α were employed to provide greater clarity as values are rarely the exact values listed above. “Suboptimal self-organization” was represented by α < 0.75; “Optimal self-organization” was represented by α = 0.75 – 1.30; “Impaired self-organization” was represented by α > 1.30. These values were based upon previously established ranges that classified populations (healthy, elderly, and impaired) as either white, pink or brown noise (Hausdorff et al., 1996, 1997, 1999; Ravi et al., 2020). The values attained during the RN with BW condition was considered the baseline because it is the natural locomotion pattern used 10% above GTV and unperturbed by an external load (no added load carriage). Change classifications were then determined for each individual (change from the baseline condition) as either “positive change” in complexity (“Suboptimal self-organization” or “Impaired self-organization” to “Optimal self-organization”), “negative change” (“Optimal self-organization” to “Suboptimal self-organization” or “Impaired self-organization”), “no change positive” (“Optimal self-organization” to “Optimal self-organization”) and “no change negative” (“Suboptimal self-organization” or “Impaired self-organization” to “Suboptimal self-organization” or “Impaired self-organization”).

Statistical Analysis

Descriptive statistics (mean and SD) were reported for all the variables. In order to determine interactive effects of load and locomotion on relative variability and gait complexity a two-way repeated measure analysis of variance (RMANOVA) for Load × Locomotion (3 × 2) was conducted separately. Additionally, to further elucidate findings regarding relative variability, tangential and perpendicular variability were assessed within locomotion pattern with a 3 × 2 (Load × Direction) RMANOVA. If interactions were significant, simple main effects were performed (paired t-tests for locomotion/direction stratified by load and RMANOVA for load stratified by locomotion/direction). If no significant interaction was observed, only main effects were analyzed. Post hoc analysis using Bonferroni-corrected pairwise comparisons were conducted when necessary.

To determine interactive effects of load and locomotion on stride regulation a three-way Load × Locomotion × Direction (3 × 2 × 2) RMANOVA was conducted on DFA scaling exponents of δ_T and δ_P. If a significant three-way interaction was observed then two-way RMANOVAs were conducted for load by locomotion (3 × 2), load by direction (3 × 2) and locomotion by direction (2 × 2). If a two-way interaction was observed, then simple main effects were analyzed (RMANOVA for load and paired t-tests for locomotion and direction). If no significant two-way interaction was observed, only main effects were analyzed. Post hoc analysis using Bonferroni-corrected pairwise comparisons were conducted when necessary. Lastly, if no significant three-way interaction was observed, only main effects were analyzed. Post hoc analysis using Bonferroni-corrected pairwise comparisons were conducted when necessary.

Partial eta squared (η²_P) was calculated as a measure of effect size given the within-subject design (Bakeman, 2005; Richardson, 2011), with magnitudes of effect interpreted as: 0.01–0.085 (small effect); 0.09–0.24 (moderate effect); and > 0.25 (large effect) (Cohen et al., 2003). Additionally, frequencies of complexity classifications and change classifications are reported to qualitatively examine individual responses. The alpha level was set at 0.05 (p ≤ 0.05).

Results

Relative Variability

See Table 2 for mean and SD of all GEM related outcomes. There was no significant interaction between load and locomotion for relative variability of motor control (F₂,₂₀ = 0.167, p = 0.85, η²_P = 0.02). Load had a significant influence on relative variability reducing the number of successful or “good” movement solutions, exemplified as relative variability magnitude decreasing as load magnitude increased confirmed by the main effect of load (F₂,₂₀ = 5.50, p = 0.01, η²_P = 0.36); with post hoc analysis revealing BW + 45% (1.28 ± 0.05) being significantly (p = 0.02) less than BW (1.55 ± 0.07). Additionally, FM demonstrated a more relative variability compared to running indicated by the main effect of locomotion (F₁,₁₀ = 8.90, p = 0.01, η²_P = 0.47) with estimated marginal means for FM (1.53 ± 0.08) being greater than RN (1.27 ± 0.04).

The interaction between load and direction during RN was not statistically significant (F₂,₂₀ = 2.33, p = 0.12, η²_p = 0.19). There was no significant main effect of load (F₂,₂₀ = 3.05, p = 0.07, η²_p = 0.23). While not significant, as load increased mean tangential (“good”) variability decreased (BW = 1.13 ± 0.10, BW + 25% = 1.10 ± 0.06, BW + 45% = 1.06 ± 0.07) and mean perpendicular (“bad”) variability increased (BW = 0.84 ± 0.13, BW + 25% = 0.88 ± 0.08, BW + 45% = 0.94 ± 0.08). However, regardless of load, tangential variability was always greater evidenced by the main effect of direction (F₁,₁₀ = 60.91, p < 0.001, η²_p = 0.86), with estimated marginal means revealing variability along the tangential (1.10 ± 0.01) was greater than along the perpendicular (0.88 ± 0.02).

There was no significant interaction between load and direction for FM (F₂,₂₀ = 1.79, p = 0.19, η²_p = 0.15). There was no main effect of load (F₂,₂₀ = 2.75, p = 0.09, η²_p = 0.22). While not significant, as load increased mean tangential (“good”) variability decreased slightly (BW = 1.20 ± 0.08, BW + 25% = 1.15 ± 0.08, BW + 45% = 1.15 ± 0.06) and mean perpendicular (“bad”) variability increased (BW = 0.74 ± 0.013, BW + 25% = 0.81 ± 0.12, BW + 45% = 0.81 ± 08). Lastly, regardless of load, tangential variability was always greater than perpendicular variability indicated by the main effect of direction (F₁,₁₀ = 90.10, p < 0.001, η²_p = 0.90), with estimated marginal means revealing tangential variability (1.17 ± 0.02) was greater than perpendicular variability (0.79 ± 0.02). See Figure 2 for exemplar plots of S_L and S_T combinations along the goal manifold.

FIGURE 2

Scaling Exponents (α): S_T & S_p

See Table 2 for all means and SD of scaling exponents. There was no significant 3-way interaction between load, locomotion and direction (F₂,₂₀ = 1.96, p = 0.17, η²_p = 0.16). As load increased control of stride-to-stride fluctuations became more strict (and therefore corrected more quickly) evidenced by the main effect of load (F₂,₂₀ = 8.87, p = 0.002, η²_p = 0.47), with post hoc analysis revealing that BW (0.6 ± 0.08) was significantly (p = 0.02) greater than BW + 45% (0.04 ± 0.11). Additionally, running exhibited less control of stride-to-stride fluctuations indicated by the main effect of locomotion (F₁,₁₀ = 8.57, p = 0.02, η²_p = 0.46), with estimated marginal means revealing that running (0.46 ± 0.07) was greater than FM (0.19 ± 0.08). Lastly, “bad” variations (perpendicular to the goal manifold) were controlled and corrected more quickly than “good” variations (tangential to goal manifold) as evidenced by the main effect of direction (F₁,₁₀ = 67.12, p < 0.001, η²_p = 0.87), with estimated marginal means revealing that persistence along the tangential (0.47 ± 0.06) was greater than the perpendicular (0.18 ± 0.07).

Complexity Analysis (α): S_L, S_T

See Table 3 for mean and SD of scaling exponents. There was no significant interaction between load and locomotion for S_L (F₂,₂₀ = 0.03, p = 0.97, η²_p = 0.003). As load magnitude increased, gait complexity decreased independent of locomotion pattern as evidence by the main effect of load (F₂,₂₀ = 8.74, p = 0.002, η²_p = 0.4 = 7), with post hoc pairwise comparisons revealing BW (0.69 ± 0.06) was significantly (p = 0.02) greater than BW + 45% (0.1 ± 0.16). Additionally, FM reduced gait complexity compared to running independent of load magnitude indicated by the main effect of locomotion (F₁,₁₀ = 7.59, p = 0.02, η²_p = 0.43), with estimated marginal means revealed running (0.59 ± 0.06) was greater than FM (0.23 ± 0.13).

There was no significant interaction between load and locomotion for S_T (F₂,₂₀ = 0.43, p = 0.36, eta = 0.10). The increase in load magnitude decreased gait complexity independent of locomotion pattern as evidenced by the main effect of load (F₂,₂₀ = 6.52, p = 0.007, η²_p = 0.40), with post hoc analysis revealing that BW (0.67 ± 0.10) was greater (p = 0.03) than BW + 45% (−0.10 ± 0.19). Lastly, FM reduced gait complexity compared to running independent of load magnitude as indicated by the main effect of locomotion (F₁,₁₀ = 9.66, p < 0.001, η²_p = 0.75), with estimated marginal means revealing running (0.76 ± 0.11) was greater than FM (−0.003 ± 0.15).

See Table 3 for frequency of observed complexity classifications by condition and Figure 3 for the frequency of each classification change of all conditions combined. For S_L and S_T 43.64% and 49.09% of the observed changes from baseline were negative changes (“Optimal self-organization” to “Suboptimal’ or ‘Impaired.” Only two participants demonstrated a positive change with 63% of those occurrences accounted for by one participant. Positive change classifications only occurred at the BW + 25% load condition and not the BW + 45%. Two of the participants with positive changes from baseline reported having prior experience in a multitude of exercise modalities that would improve lower limb muscular endurance and anaerobic conditioning (see Table 1). For S_L only 10.91% of changes were classified as “no change positive” with 58% of the occurrences accounted for by two individuals. Both individuals performed a greater variety of exercise modalities including resistance training and trained more frequently and for longer periods of time. Lastly, only a single subject demonstrated an “Impaired” classification at baseline (RN at BW). This participant trained for the shortest durations (30 min maximum) only 2–3 times a week (see Table 1 for more detailed characteristics).

FIGURE 3

Discussion

The objectives of this investigation were to determine the interactive effects of load magnitude and locomotion pattern on motor variability, stride-to-stride control/regulation and complexity. Load magnitude significantly altered relative variability independent of locomotion pattern evidenced by the significant main effect of load. As load increased relative variability decreased with BW + 45% (1.28σ) having 21% less relative variability compared to BW (1.55σ), suggesting individuals are better able to leverage either the system degeneracy, or indeed, redundancy of coordinative patterns employed to execute a stride during unloaded bipedal ambulatory tasks. The higher relative variability ratios are achieved by a greater variance tangential the goal manifold (“good” variability) and a reduction in variance perpendicular to the goal manifold (“bad” variability) [see Figure 2]. Importantly, these findings alone do not necessarily indicate more stable performance (task execution), but that during unloaded conditions there is a greater workspace of solutions that can be utilized to accomplish the goal task (maintain velocity). Coupled with the observed complexity scaling exponents ∼1 for S_L/S_T at BW and BW + 25% load conditions suggest that individuals display an optimally organized locomotor system fully leveraging equifinality. Furthermore, the lack of significant differences between BW and BW + 25% load conditions indicates that the latter load condition is insufficient to impose a constraint on the locomotor system that results in altered motor variability. At BW + 45%, the locomotor system has a diminished ability to leverage “good” motor variability. The lack of observed “good” variability, may limit energy dispersion across multiple structures up the kinetic chain when energy is highest [due to increased forces from load (Birrell et al., 2007; Lloyd, 2010; Liew et al., 2016; Willy et al., 2016, 2019; Krajewski et al., 2020)] increasing cumulative mechanic stress MSI risk (Nordin and Dufek, 2019; Baggaley et al., 2020).

Contrary to the proposed hypothesis, the FM (RV = 1.52) resulted in 19.6% greater relative variability than running (RV = 1.27). However, this finding alone should be interpreted with caution; the greater variability simply denotes a larger set of motor solutions observed for FM compared to running. Similar observations were made by Black et al. (2007) demonstrating greater motor variability in adolescents with down syndrome compared to aged-matched neurotypical adolescents, indicating that during unperturbed steady-state gait, those with down-syndrome needed to utilize their full movement solution potential in order to maintain velocity (Black et al., 2007). In addition, Black et al. (2007) demonstrated that the greater motor variability exhibited during unperturbed gait suggested a locomotor system operating closer to the action boundaries that would be more likely to fail (fall) with additional pertubrations (Black et al., 2007). Thus, greater relative variability alone may not indicate more task performance stability or motor-system health. In females during FM there is a greater contribution of frontal plane moments (adduction/abduction) to knee total joint moment, indicative of potentially deleterious movements (Krajewski et al., 2020). Although a greater range of movement solutions were utilized during FM compared to running, they may have included more movement solutions that while successful in executing the goal task (maintaining velocity), are maladaptive/deleterious with respect to joint/tissue health (Zabala et al., 2013; Asay et al., 2018; Baggaley et al., 2020; Krajewski et al., 2020). Future work should focus on elucidating the movement pattern structures (segment couplings) in more detail, which coincide with relative variability, its magnitude and structure.

Additionally, the unfamiliarity of the task for this sample (they had little to no experience with FM) may have demonstrated a greater exploration of state-space as a consequence of learning to perform this locomotion pattern (Newell and Vaillancourt, 2001; Pacheco and Newell, 2017). The greater relative variability exhibited by FM compared to running may be indicative of a locomotor system trying to adapt to improve mechanical efficiency or reduce pain sensations associated with each stride. By obstructing the natural bifurcation of locomotion (participants executing FM at a velocity 10% above their GTV), the system is in search of an order parameter (segment coupling pattern) that adheres to cost function (decrease metabolic cost/increase mechanical efficiency) of the task control parameter (horizontal velocity) (McGinnis and Newell, 1982; Newell and Vaillancourt, 2001; Davids et al., 2003; Caballero et al., 2019; Shafizadeh et al., 2019). Consequently, the observed complexity for FM was α = 0.11 (“suboptimal”), supporting the notion that the observed greater relative variability in FM compared to running, represented a state-space exploratory behavior. In conjunction with the low system complexity, the locomotor system is less stable to additional perturbations. Whilst imposed locomotion (FM) and load carriage are themselves perturbations dynamic military environments present even more perturbations (uneven terrain, enemy threats, orders from local commanders). If motor learning occupies a large percentage of the system workspace, it potentially results in the failure to identify external perturbations [due to competition over feedback/feedforward resources (Woollacott and Shumway-Cook, 2002; Al-Yahya et al., 2011; De Sanctis et al., 2014)] making the locomotor system less stable and more susceptible to failure (slips/trips/falls, identification of threats). Future research should compare novice and experienced individuals with load carriage to determine if there is a difference in relative variability and complexity while completing ambulatory tasks with load.

Stride-to-stride regulation demonstrated less control for unloaded conditions as evidenced by the main effect of load (p = 0.002). The BW + 45% load condition exhibited (δ_Tα = 0.22, δ_pα = −0.11) significantly more regulation than BW (δ_Tα = 0.74, δ_pα = 0.45). Alpha coefficient (α) values less than.5 indicate statistical anti-persistence, representing much stricter control because the subsequent stride variation is more likely to be the opposite composition (combination of S_L and S_T) than the previous. Thus, movements were corrected much quicker and more often with the addition of substantial (BW + 45%) load carriage independent of locomotion pattern. Moreover, the main effect of direction (p < 0.001) exhibited more strict regulation perpendicular to the manifold (“bad” variability) compared to tangential (“good” variability) [see Table 1 for δ_T (α) and δ_p (α) means]. The combination of less control of “good” variability and more control of “bad” variability [resembling an ideal minimum intervention principle (MIP) model [δ_Tα > 1; δ_pα < 0.5] (Cusumano and Dingwell, 2013)] is indicative of a stride-regulation strategy of a healthy system recognizing movement variations that impede the execution of the task goal (maintain velocity). Further, the “system controller” minimally intervenes, minimizing control effort theoretically freeing up system capacity for other components of the locomotor system (Cusumano and Dingwell, 2013). In fact, the BW results (δ_Tα = 0.74, δ_pα = 0.45) of this investigation were similar to stride regulation findings of healthy adults (δ_Tα = ∼0.90, δ_pα = ∼0.42) by Dingwell et al. (2017). If the individual were unhealthy, a failure to properly regulate “bad” motor variability potentially results in them unable to successfully execute the task. However, as load increased both δ_Tα and δ_pα decreased indicating a change in stride-regulation strategy that reflected an absolute position control (POS) model [δ_Tα and δ_pα < 0.5] (Cusumano and Dingwell, 2013), postulating that maximal control effort was used thus reducing the capacity of other locomotor system components (Cusumano and Dingwell, 2013).

Locomotion pattern also affected stride-to-stride regulation independent of load, with FM (δ_Tα = 0.36, δ_pα = 0.04) demonstrating stricter control compared to running (δ_Tα = 0.60, δ_pα = 0.35). Even “good” motor variability (δ_Tα) was tightly regulated for FM evidenced by the δ_Tα < 0.5. Once again, the stride-regulation strategy utilized for FM at all load conditions more closely mimicked that of a POS model (Cusumano and Dingwell, 2013). A stricter regulation of stride-to-stride intervals coupled with greater relative variability of FM at BW + 45% compared to RN BW + 45% may further evidence the participant learning by exploring state-space and freezing/unlocking different degrees of freedom quickly (Bernstein, 1967; Sporns and Edelman, 1993; Newell and Vaillancourt, 2001). Thus in response to the load and imposed locomotion perturbations a structural reorganization of the locomotor system occurs, with potentially greater reliance on supraspinal input that disrupt feed forward mechanisms of gait (further supported by the observed α < 0.5 for FM of S_L and S_T) due to the competition over feedback/feedforward resources (Woollacott and Shumway-Cook, 2002; Al-Yahya et al., 2011; De Sanctis et al., 2014). It is likely that the perturbations associated with FM and load elicit more system capacity dedicated to control in a relatively unperturbed state (walking on a treadmill) potentially overpowering other system components important to navigating dynamic environments. This may have important consequences to military populations, as loaded ambulatory tasks are often undertaken in dynamic circumstances requiring the integration of multiple information sources to understand the context within which the action takes place, i.e., quickly navigating across an open area of the battle space while aiming and firing a weapon.

When assessing complexity of gait dynamics (S_L and S_T) there was a main effect of load (p = 0.002 and p = 0.007), with complexity decreasing as load increased. The BW + 45% condition (S_Lα = 0.1; S_Tα = −0.1) exhibited significant less complexity than BW (S_Lα = 0.69; S_Tα = 0.67). Whilst, group mean comparison demonstrated a decrease in complexity, the group mean for BW still represents a “suboptimal” organization (α < 0.75) (Hausdorff et al., 1997; Hausdorff, 2007). However, when observing the individual results, eight of the participants had “optimal” complexity or.75 < α < 1.30 during BW load conditions (see Tables 1, 3). Moreover, only five participants had “optimal” complexity for BW + 25% and every participant exhibited “suboptimal” complexity for the BW + 45% load condition. “Optimal” complexity or scaling exponents (α) ∼1 represent long-range correlations (pink noise) that is indicative of skilled performance and utilization of prior stride information (i.e., proprioceptive) to influence future strides (Hausdorff et al., 1996, 1997; Delignières and Torre, 2009; Nourrit-Lucas et al., 2015). Likewise, the fractal structure of “optimal” represents an independence of fluctuations at different time scales meaning a perturbation of one system component will not likely affect the global system (locomotor system as a whole) (West and Shlesinger, 1989; Bak and Paczuski, 1995; Marks-Tarlow, 1999; Torre et al., 2007). Therefore, the data observed within the present study indicates that a load carriage magnitude of at least BW + 45% reduces the system stability and adaptability, predisposing the system to failure (falling) in the presence of additional perturbations (i.e., increased fall risk with uneven terrain) in females with limited/no load carriage experience. Importantly, BW + 45% for this female sample was 26.6 ± 4.7kg, which represents a typical combat load (20 kg) (Taylor et al., 2016), a load used during operations that will bombard individuals with perturbations (terrain obstacles[debris], enemy threats, officer commands/questions) that if not actioned correctly could result in serious MSI or death. The latter suggests the need for further research to determine if experience/training with load carriage improves system complexity that can better handle additional adaptations.

In addition to a main effect of load, locomotion patterns also affected the complexity of S_L and S_T (p = 0.02 and p < 0.001, respectively). As hypothesized, the natural locomotion pattern running (S_Lα = 0.59; S_Tα = 0.76), exhibited greater complexity than FM (S_Lα = 0.23; S_Tα < −0.01). Considering that running is the extant locomotion observed at a velocity exceeding GTV, it is not surprising more participants (seven) had “optimal” complexity during RN conditions (BW and BW + 25% only). Five of the participants did have “optimal” complexity for FM but at the BW and BW + 25% conditions only. These five participants (S1, S7, S8, S10, and S11) engaged in moderate amounts of exercise time per week and engaged in a multitude of exercise modalities that included a form of anaerobic conditioning and/or strength training (see Table 1). Furthermore, one participant (S3) had an “impaired” complexity at baseline (running at BW). Interestingly, this participant performed the least amount of exercise time per week and engaged in the least varied modes of exercise [see Table 1]. Thus, the movement poverty in terms of time and coordination diversity, may impact system adaptability as reflected in the “impaired” and “suboptimal” complexity during minimally perturbed steady-state behaviors. Moreover, the complexity outcomes demonstrated varying individual responses highlighting its potential utility as a mechanical “biomarker” of the current state of the locomotor system and training adaptation response. Future research should compare complexity of different fitness level groups (i.e., highly aerobic versus highly anaerobic fit individuals) during loaded ambulatory tasks to further elucidate characteristics of “optimal” performers.

The primary limitations of this study are its small, sex specific (female) sample and the number of consecutive data points (∼130). While this investigation cannot make inferences about sex specific responses to load carriage in the absence of a male cohort, females are an under represented population in load carriage literature (Loverro et al., 2019). Therefore, findings regarding motor variability, stride regulation and gait complexity of this study can only be generalized to healthy recruit-aged females (18–33 years old) whom are novices to load carriage and forced marching. Future research should compare males versus females and recruit versus experienced individuals (deployable soldiers) to determine if the unfamiliarity to the tasks and equipment were confounding the observed responses to load magnitude and locomotion pattern. Likewise, it is advised to perform fractal analysis with a minimum of 512 consecutive data points (Delignières et al., 2006), however, 128 consecutive strides has been determined to be within 6% of the actual scaling value (Hausdorff et al., 1997). In addition, performing steady-state behavior with load carriage while minimizing the effects of fatigue is difficult beyond several minute trials. Moreover, our findings on the effect of locomotion pattern on gait complexity are similar to investigations regarding imposed frequencies and complexity (Delignieres et al., 2004; Ducharme et al., 2018). Nevertheless, while not ideal for complexity, the ∼130 consecutive data points is a robust time-series for GEM decomposition (Dingwell et al., 2010) and provides the first quantitative data on load magnitude’s influence on gait complexity. Lastly, upon the completion of the study we learned that one participant (S5) completed the experimental protocols with a distal avulsion of the semitendinosus. Interestingly, this participant exhibited “optimal” complexity for running at BW and BW + 25%. Although alone this data is inconclusive and may represent an isolated incident, but it does suggest that gait complexity of S_L and S_T may be an inappropriate factor to assess musculoskeletal health. The complexity of S_L and S_T specifically may only represent the global function of the locomotor system at achieving the task goal.

In conclusion, there are no interactive effects of load magnitude and locomotion pattern on motor variability, stride regulation and gait complexity. But load and locomotion do independently alter the function of the locomotor system. As load increases there is a reduction in relative variability (good:bad motor variability), gait complexity (α < 0.5) and strides become more tightly controlled. FM further reduces gait complexity and mimics stride regulation strategies of BW + 45% load conditions. Moreover, despite more relative variability for FM compared to running, this appears to be a consequence of state-space exploration, as supported by the increased stride control/error correction and “suboptimal” complexity. While complexity of S_L and S_T may be indicative of locomotor system function in terms of achieving a task goal (maintaining horizontal velocity in the case of this investigation), the variability/complexity of these factors do not appear to represent the health of the musculoskeletal system (i.e., state of joint/tissue health and whether an injurious movement pattern is being used). Therefore, additional order parameters (different gait variables) should be investigated to identify a marker of global MSI risk. Likewise, research should be conducted with longer trials to confirm the complexity findings of this investigation and elucidate the role of fatigue. Soldiers ultimately operate in dynamic environments with lots of perturbations under substantial load carriage where poor movement/slips, trips and fall can produce MSI or even death. The findings of this study indicate that locomotor system function is altered by FM and BW + 45% load, resulting in reduced motor variability and a system with less stability/adaptability that is potentially more susceptible to failure with additional perturbations.

References

• 1

  AlexanderR. M. (1992). A model of bipedal locomotion on compliant legs.Philos. Trans. R. Soc. Lond. B Biol. Sci.338189–198. 10.1098/rstb.1992.0138

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  Al-YahyaE.DawesH.SmithL.DennisA.HowellsK.CockburnJ. (2011). Cognitive motor interference while walking: a systematic review and meta-analysis.Neurosci. Biobehav. Rev.35715–728. 10.1016/j.neubiorev.2010.08.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  AsayJ. L.Erhart-HledikJ. C.AndriacchiT. P. (2018). Changes in the total knee joint moment in patients with medial compartment knee osteoarthritis over 5 years.J. Orthop Res.362373–2379. 10.1002/jor.23908

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  Atlas Collaboration, AadG.AbbottB.AbdallahJ.Abdel KhalekS.AbdinovO.et al (2014). Measurement of the muon reconstruction performance of the ATLAS detector using 2011 and 2012 LHC proton-proton collision data.Eur. Phys. J. C Part Fields74:3130.

  • Google Scholar
  • View reference in article

• 5

  AttwellsR. L.BirrellS. A.HooperR. H.MansfieldN. J. (2006). Influence of carrying heavy loads on soldiers’ posture, movements and gait.Ergonomics491527–1537. 10.1080/00140130600757237

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  BaggaleyM.EspositoM.XuC.UnnikrishnanG.ReifmanJ.EdwardsW. B. (2020). Effects of load carriage on biomechanical variables associated with tibial stress fractures in running.Gait Posture77190–194. 10.1016/j.gaitpost.2020.01.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  BaidaS. R.GoreS. J.Franklyn-MillerA. D.MoranK. A. (2018). Does the amount of lower extremity movement variability differ between injured and uninjured populations? A systematic review.Scand. J. Med. Sci. Sports281320–1338. 10.1111/sms.13036

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  BakP.PaczuskiM. (1995). Complexity, contingency, and criticality.Proc. Natl. Acad. Sci. U.S.A.926689–6696. 10.1073/pnas.92.15.6689

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  BakP.TangC.WiesenfeldK. (1987). Self-organized criticality: an explanation of the 1/f noise.Phys. Rev. Lett.59381–384. 10.1103/physrevlett.59.381

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  BakemanR. (2005). Recommended effect size statistics for repeated measures designs.Behav. Res. Methods37379–384. 10.3758/bf03192707

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  BentL. R.InglisJ. T.McFadyenB. J. (2004a). When is vestibular information important during walking?J. Neurophysiol.921269–1275. 10.1152/jn.01260.2003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  BentL. R.McFadyenB. J.InglisJ. T. (2004b). Is the use of vestibular information weighted differently across the initiation of walking?Exp. Brain Res.157407–416.

  • Google Scholar
  • View reference in article

• 13

  BernsteinN. A. (1967). The Coordination and Regulation of Movements.London: Pergamon Press.

  • Google Scholar
  • View reference in article

• 14

  BirrellS. A.HooperR. H.HaslamR. A. (2007). The effect of military load carriage on ground reaction forces.Gait Posture26611–614. 10.1016/j.gaitpost.2006.12.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  BizziE.CheungV. C. (2013). The neural origin of muscle synergies.Front. Comput. Neurosci.7:51. 10.3389/fncom.2013.00051

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  BlackD.SmithB. A.WuJ.UlrichB. D. (2007). Uncontrolled manifold analysis of segmental angle variability during walking: preadolescents with and without down syndrome.Exp. Brain Res.183511–521. 10.1007/s00221-007-1066-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  CaballeroC.DavidsK.HellerB.WheatJ.MorenoF. J. (2019). Movement variability emerges in gait as adaptation to task constraints in dynamic environments.Gait Posture701–5. 10.1016/j.gaitpost.2019.02.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  CohenJ.CohenP.WestS. G.AikenL. S. (2003). Applied Multiple Regression/Correlation Analysis for the Behavioral Sciences, 3rd Edn. New York, NY: Routledge.

  • Google Scholar
  • View reference in article

• 19

  CusumanoJ. P.DingwellJ. B. (2013). Movement variability near goal equivalent manifolds: fluctuations, control, and model-based analysis.Hum. Mov. Sci.32899–923. 10.1016/j.humov.2013.07.019

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  DavidsK.GlazierP.AraujoD.BartlettR. (2003). Movement systems as dynamical systems: the functional role of variability and its implications for sports medicine.Sports Med.33245–260. 10.2165/00007256-200333040-00001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  De SanctisP.ButlerJ. S.MalcolmB. R.FoxeJ. J. (2014). Recalibration of inhibitory control systems during walking-related dual-task interference: a mobile brain-body imaging (MOBI) study.Neuroimage9455–64. 10.1016/j.neuroimage.2014.03.016

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  DelignieresD.LemoineL.TorreK. (2004). Time intervals production in tapping and oscillatory motion.Hum. Mov. Sci.2387–103. 10.1016/j.humov.2004.07.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  DelignièresD.RamdaniS.LemoineL.TorreK.FortesM.NinotG. (2006). Fractal analyses for ‘short’ time series: a re-assessment of classical methods.J. Math. Psychol.50525–544. 10.1016/j.jmp.2006.07.004

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  DelignièresD.TorreK. (2009). Fractal dynamics of human gait: a reassessment of Hausdorff et al. (1996)’s data.J. Appl. Physiol.1061272–1279. 10.1152/japplphysiol.90757.2008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  DingwellJ.JohnJ.CusumanoJ. P. (2010). Do humans optimally exploit redundancy to control step variability in walking?PLoS Comput. Biol.6:e1000856. 10.1371/journal.pcbi.1000856

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  DingwellJ.SalinasM.CusumanoJ. (2017). Increased gait variability may not imply impaired stride-to-stride control of walking in healthy older adults.Gait Posture55131–137. 10.1016/j.gaitpost.2017.03.018

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  DucharmeS.LiddyJ. J.HaddadJ. M.BusaM. A.ClaxtonL. J.Van EmmerikR. E. (2018). Association between stride time fractality and gait adaptability during unperturbed and asymmetric walking.Hum. Mov. Sci.58248–259. 10.1016/j.humov.2018.02.011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  GatesD. H.DingwellJ. B. (2008). The effects of neuromuscular fatigue on task performance during repetitive goal-directed movements.Exp. Brain Res.187573–585. 10.1007/s00221-008-1326-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  GelfandI. M.LatashM. L. (1998). On the problem of adequate language in motor control.Motor Control2306–313. 10.1123/mcj.2.4.306

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  GibsonJ. J. (1966). The Senses Considered as Perceptual Systems.Boston, MA: Houghton Mifflin.

  • Google Scholar
  • View reference in article

• 31

  GibsonJ. J. (1979). The Ecological Approach to Visual Perception.Boston, MA: Houghton Mifflin.

  • Google Scholar
  • View reference in article

• 32

  GisigerT. (2001). Scale invariance in biology: Coincidence or footprint of a universal mechanism?Biol. Rev. Camb. Philos. Soc.76161–209. 10.1017/s1464793101005607

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  GoldbergerA.PengC. K.LipstizL. A. (2005). What is physiologic complexity and how does it change with aging and disease?Neurobiol. Aging2323–26. 10.1016/s0197-4580(01)00266-4

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  HausdorffJ.PengC. K.LadinZ. V. I.WeiJ. Y.GoldbergerA. L. (1995). Is walking a random walk? Evidence for long-range correlations in stride interval of human gait.J. Appl. Physiol.349–358. 10.1152/jappl.1995.78.1.349

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  HausdorffJ. M. (2007). Gait dynamics, fractals and falls: finding meaning in the stride-to-stride fluctuations of human walking.Hum. Mov. Sci.26555–589. 10.1016/j.humov.2007.05.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  HausdorffJ. M. (2009). Gait dynamics in Parkinson’s disease: common and distinct behavior among stride length, gait variability, and fractal-like scaling.Chaos19:026113. 10.1063/1.3147408

  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  HausdorffJ. M.MitchellS. L.FirtionR.PengC. K.CudkowiczM. E.WeiJ. Y.et al (1997). Altered fractal dynamics of gait: reduced stride-interval correlations with aging and Huntington’s disease.J. Appl. Physiol.82262–269. 10.1152/jappl.1997.82.1.262

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  HausdorffJ. M.PurdonP. L.PengC. K.LadinZ.WeiJ. Y.GoldbergerA. L. (1996). Fractal dynamics of human gait: stability of long-range correlations in stride interval fluctuations.J. Appl. Physiol.801448–1457. 10.1152/jappl.1996.80.5.1448

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  HausdorffJ. M.ZemanyL.PengC.GoldbergerA. L. (1999). Maturation of gait dynamics: stride-to-stride variability and its temporal organization in children.J. Appl. Physiol.861040–1047.

  • Google Scholar
  • View reference in article

• 40

  HollandsM. A.Marple-HorvatD. E. (1996). Visually guided stepping under conditions of step cycle-related denial of visual information.Exp. Brain Res.109343–356.

  • Google Scholar

• 41

  HuangT.KuoA. D. (2014). Mechanics and energetics of load carriage during human walking.J. Exp. Biol.217605–613. 10.1242/jeb.091587

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  HuntN.McGrathD.StergiouN. (2014). The influence of auditory-motor coupling on fractal dynamics in human gait.Sci. Rep.4:5879.

  • Google Scholar

• 43

  IyengarN.PengC. K.MorinR.GoldbergerA. L.LipsitzL. A. (1996). Age-related alterations in the fractal scaling of cardiac interbeat interval dynamics.Am. J. Physiol. Regul. Integr. Comp. Physiol.2711078–1084.

  • Google Scholar

• 44

  JensenA.LairdM.JamesonJ. T.KellyK. R. (2019). Prevalence of musculoskeletal injuries sustained during marine corps recruit training.Mil. Med.184(Suppl. 1), 511–520. 10.1093/milmed/usy387

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  KimY.TagawaY.ObinataG.HaseK. (2011). Robust control of CPG-based 3D neuromusculoskeletal walking model.Biol. Cybern.105269–282. 10.1007/s00422-011-0464-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  KnapikJ. J.ReynoldsK. L.HarmanE. (2004). Soldier load carriage: historical, physiological, biomechanical, and medical aspects.Mil. Med.16945–56.

  • Google Scholar

• 47

  KrajewskiK. T.DeverD. E.JohnsonC. C.RawcliffeA. J.AhamedN. U.FlanaganS. D.et al (2020). Load carriage magnitude and locomotion strategy alter knee total joint moment during bipedal ambulatory tasks in recruit-aged women.J. Biomech.105:109772. 10.1016/j.jbiomech.2020.109772

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  LaFiandraM.WagenaarR. C.HolfK. G.ObusekJ. P. (2003). How do load carriage and walking speed influence trunk coordination and stride parameters?J. Biomech.3687–95. 10.1016/s0021-9290(02)00243-9

  • CrossRef
  • Google Scholar

• 49

  LatashM. (2016). Biomechanics as a window into the neural control of movement.J. Hum. Kinet.527–20. 10.1515/hukin-2015-0190

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  LatashM.LevinM. F.ScholzJ. P.SchönerG. (2010). Motor control theories and their applications.Medicina46382–392. 10.3390/medicina46060054

  • CrossRef
  • Google Scholar

• 51

  LatashM.ScholzJ. P.SchönerG. (2002). Motor control strategies revealed in the structure of motor variability.Exerc. Sports Sci. Rev.3026–31. 10.1097/00003677-200201000-00006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  LentonG.DoyleT. L. A.LloydD. G.HiggsJ.BillingD.SaxbyD. J. (2019). Lower-limb joint work and power are modulated during load carriage based on load configuration and walking speed.J. Biomech.83174–180. 10.1016/j.jbiomech.2018.11.036

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  LiewB.MorrisS.NettoK. (2016). The effects of load carriage on joint work at different running velocities.J. Biomech.493275–3280. 10.1016/j.jbiomech.2016.08.012

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  LloydR. (2010). Kinetic changes associated with load carriage using two rucksack designs.Ergonomics431331–1341. 10.1080/001401300421770

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  LovalekarM.KeenanK. A.BealsK.NindlB. C.PihokerA. A.ColemanL. C.et al (2020). Incidence and pattern of musculoskeletal injuries among women and men during marine corps training in sex-integrated units.J. Sci. Med. Sport19S1440–S2440.

  • Google Scholar

• 56

  LoverroK.HasselquistL.LewisC. L. (2019). Females and males use different hip and knee mechanics in response to symmertric military-relevant loads.J. Biomech.95:109280. 10.1016/j.jbiomech.2019.07.024

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  Marks-TarlowT. (1999). The self as a dynamical system nonlinear dynamics.Psychol. Life Sci.3311–345.

  • Google Scholar

• 58

  MatthisJ. S.BartonS. L.FajenB. R. (2017). The critical phase for visual control of human walking over complex terrain.Proc. Natl. Acad. Sci. U.S.A.114E6720–E6729.

  • Google Scholar

• 59

  McGinnisP.NewellK. M. (1982). Topological dynamics: a framework for describing movement and its constraints.Hum. Mov. Sci.1289–305. 10.1016/0167-9457(82)90017-3

  • CrossRef
  • Google Scholar

• 60

  MinassianK.HofstoetterU. S.DzeladiniF.GuertinP. A.IjspeertA. (2017). The human central pattern generator for locomotion: Does it exist and contribute to walking?Neuroscientist23649–663. 10.1177/1073858417699790

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  MolloyJ. M.PendergrassT. L.LeeI. E.ChervakM. C.HauretK. G.RhonD. I. (2020). Musculoskeletal injuries and united states army readiness part I: overview of injuries and their strategic impact.Mil. Med.10.1093/milmed/usaa027

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  MoonY.SungJ.AnR.HernandezM. E.SosnoffJ. J. (2016). Gait variability in people with neurological disorders: a systematic review and meta-analysis.Hum. Mov. Sci.47197–208. 10.1016/j.humov.2016.03.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63

  MukherjeeM.YentesJ. M. (2018). Movement variability: a perspective on success in sports, health, and life.Scand. J. Med. Sci. Sports28758–759. 10.1111/sms.13038

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  NewellK. M.VaillancourtD. E. (2001). Dimensional change in motor learning.Hum. Mov. Sci.20695–715. 10.1016/s0167-9457(01)00073-2

  • CrossRef
  • Google Scholar

• 65

  NordinA.DufekJ. S. (2019). Reviewing the variability-overuse injury hypothesis: Does movement variability relate to landing injuries?Res. Q. Exerc. Sport90190–205. 10.1080/02701367.2019.1576837

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66

  Nourrit-LucasD.TossaA. O.ZelicG.DelignieresD. (2015). Learning, motor skill, and long-range correlations.J. Mot. Behav.47182–189. 10.1080/00222895.2014.967655

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 67

  PachecoM.NewellK. M. (2017). Search strategies in practice: Influence of information and task constraints.Acta Psychol.1829–20. 10.1016/j.actpsy.2017.11.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68

  PandyM. G.AndriacchiT. P. (2010). Muscle and joint function in human locomotion.Annu. Rev. Biomed. Eng.12401–433. 10.1146/annurev-bioeng-070909-105259

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69

  PengC.MietusJ.HausdorffJ. M.HavlinS.StanleyH. E.GoldbergerA. L. (1993). Long-range anti-correlations and non-Gaussian behavior of the heartbeat.Phys. Rev. Lett.701343–1346. 10.1103/physrevlett.70.1343

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70

  PengC. K.HavlinS.StanleyH. E.GoldbergerA. L. (1995). Quantification of scaling exponents and crossove phenomena in nonstationary heartbeat time series.Chaos582–87. 10.1063/1.166141

  • CrossRef
  • Google Scholar

• 71

  PetersR. M.DaltonB. H.BlouinJ. S.InglisJ. T. (2017). Precise coding of ankle angle and velocity by human calf muscle spindles.Neuroscience34998–105. 10.1016/j.neuroscience.2017.02.034

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72

  RaviD. K.MarmelatV.TaylorW. R.NewellK. M.StergiouN.SinghN. B. (2020). Assessing the temporal organization of walking variability: a systematic review and consensus guidelines on detrended fluctuation analysis.Front. Physiol.11:562. 10.3389/fphys.2020.00562

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  RichardsonJ. T. E. (2011). Eta squared and partial eta squared as measures of effect size in educational research.Educ. Res. Rev.6135–147. 10.1016/j.edurev.2010.12.001

  • CrossRef
  • Google Scholar

• 74

  RoerdinkM.de JongeC. P.SmidL. M.DaffertshoferA. (2019). Tightening up the control of treadmill walking: effects of maneuverability range and acoustic pacing on stride-to-stride fluctuations.Front. Physiol.10:257. 10.3389/fphys.2019.00257

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 75

  RossignolS.DubucR.GossardJ. P. (2006). Dynamic sensorimotor interactions in locomotion.Physiol. Rev.8689–154. 10.1152/physrev.00028.2005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76

  SánchezC. C.MurilloD. B.DavidsK.HernándezF. J. M. (2016). Variations in task constraints shape emergent performance outcomes and complexity levels in balancing.Exp. Brain Res.2341611–1622. 10.1007/s00221-016-4563-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 77

  SeayJ.FellinR. E.SauerS. G.FrykmanP. N.BenselC. K. (2014a). Lower extremity biomechanical changes associated with symmetrical torso loading during simulated marching.Mil. Med.17985–91. 10.7205/milmed-d-13-00090

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 78

  SeayJ.FrykmanP. N.SauerS. G.GutekunstD. J. (2014b). Lower extremity mechanics during marching at three different cadences for 60 minutes.J. Appl. Biomech.3021–30. 10.1123/jab.2012-0090

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 79

  SedighiA.NussbaumM. A. (2019). Exploration of different classes of metrics to characterize motor variability during repetitive symmetric and asymmetric lifting tasks.Sci. Rep.9:9821.

  • Google Scholar

• 80

  SeifertL.ButtonC.DavidsK. (2013). Key properties of expert movement systems in sport : an ecological dynamics perspective.Sports Med.43167–178. 10.1007/s40279-012-0011-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 81

  ShafizadehM.CrowtherR.WheatJ.DavidsK. (2019). Effects of personal and task constraints on limb coordination during walking: a systematic review and meta-analysis.Clin. Biomech.611–10. 10.1016/j.clinbiomech.2018.10.024

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 82

  SimpsonK. M.MunroB. J.SteeleJ. R. (2011). Effect of load mass on posture, heart rate and subjective responses of recreational female hikers to prolonged load carriage.Appl. Ergon.42403–410. 10.1016/j.apergo.2010.08.018

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 83

  SimpsonR.GrahamS. M.Forida-JamesG. D.ConnaboyC.ClementR.JacksonA. S. (2010). Perceived exertion and heart rate models for estimating metabolic workload in elite British soldiers performing a backpack load-carriage task.Appl. Physiol. Nutr. Metab.35650–656. 10.1139/h10-053

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 84

  SimpsonR. J.GrahamS. M.ConnaboyC.ClementR.PolloniniL.Florida-JamesG. D. (2017). Blood lactate thresholds and walking/running economy are determinants of backpack-running performance in trained soldiers.Appl. Ergon.58566–572. 10.1016/j.apergo.2016.04.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 85

  SpornsO.EdelmanG. M. (1993). Solving Bernstein’s problem: a proposal for the development of coordinated movement by selection.Child Dev.64960–981. 10.2307/1131321

  • CrossRef
  • Google Scholar

• 86

  StadnitskiT. (2012). Measuring fractality.Front. Physiol.3:127. 10.3389/fphys.2012.00127

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 87

  TaylorN. A.PeoplesG. E.PetersenS. R. (2016). Load carriage, human performance, and employment standards.Appl. Physiol. Nutr. Metab.41(6 Suppl. 2), S131–S147.

  • Google Scholar

• 88

  TelferS.LangeM. J.SudduthA. S. M. (2017). Factors influencing knee adduction moment measurement: a systematic review and meta-regression analysis.Gait Posture58333–339. 10.1016/j.gaitpost.2017.08.025

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 89

  TerrierP.DériazO. (2012). Persistent and anti-persistent pattern in stride-to-stride variability of treadmill walking: influence of rhythmic auditory cueing.Hum. Mov. Sci.311585–1597. 10.1016/j.humov.2012.05.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 90

  TerrierP.TurnerV.SchutzY. (2005). GPS analysis of human locomotion: further evidence for long-range correlations in stride-to-stride fluctuations of gait parameters.Hum. Mov. Sci.2497–115. 10.1016/j.humov.2005.03.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 91

  TorreK.BalasubramaniamR. (2009). Disentangling stability, variability and adaptability in human performance: focus on the interplay between local variance and serial correlation.J. Exp. Psychol. Hum. Percept. Perform.37539–50. 10.1037/a0020556

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 92

  TorreK.DelignièresD.LemoineL. (2007). Detection of long-range dependence and estimation of fractal exponents through ARFIMA modelling.Br. J. Math. Stat. Psychol.6085–106. 10.1348/000711005x89513

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 93

  TorreK.VergotteG.VielÉ.PerreyS.DupeyronA. (2019). Fractal properties in sensorimotor variability unveil internal adaptations of the organism before symptomatic functional decline.Sci. Rep.9:15736.

  • Google Scholar

• 94

  TurveyM. (1990). Coordination.Am. Psychol.45938–953.

  • Google Scholar

• 95

  Van OrdenG.KloosH.WallotS. (2009). Living in the pink: intentionality wellbeing and complexity.Handb. Philos. Sci.10639–689.

  • Google Scholar

• 96

  Van OrdenG. C.HoldenJ. G.TurveyM. T. (2003). Self-organization of cognitive performance.J. Exp. Psychol. Gen.132331–350.

  • Google Scholar

• 97

  WalshG. S.LowD. C.ArkesteijnM. (2018). Effect of stable and unstable load carriage on walking gait variability, dynamic stability and muscle activity of older adults.J. Biomech.7318–23. 10.1016/j.jbiomech.2018.03.018

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 98

  WarrenW. H.KayB. A.ZoshW. D.DuchonA. P.SahucS. (2001). Optic flow is used to control human walking.Nat. Neurosci.4213–216. 10.1038/84054

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 99

  WestB. J.ScafettaN. (2003). Nonlinear dynamical model of human gait.Phys. Rev. E Stat. Nonlin. Soft Matter Phys.67(5 Pt 1), 051917.

  • Google Scholar

• 100

  WestB. J.ShlesingerM. F. (1989). On the ubiquity of 1/f noise.Int. J. Mod. Phys. B.3795–819.

  • Google Scholar

• 101

  WillsJ.SaxbyD. J.LentonG. K.DoyleT. L. A. (2019). Ankle and knee moment and power adaptations are elicited through load carriage conditioning in males.J. Biomech.97:109341. 10.1016/j.jbiomech.2019.109341

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 102

  WillyR. W.DeVitaP.MeardonS. A.BaggaleyM.WombleC. C.WillsonJ. D. (2019). Effects of load carriage and step length manipulation on Achilles tendon and knee loads.Mil. Med.184e482–e489.

  • Google Scholar

• 103

  WillyR. W.WillsonJ. D.ClowersK.BaggaleyM.MurrayN. (2016). The effects of body-borne loads and cadence manipulation on patellofemoral and tibiofemoral joint kinetics during running.J. Biomech.494028–4033. 10.1016/j.jbiomech.2016.10.043

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 104

  WoollacottM.Shumway-CookA. (2002). Attention and the control of posture and gait: a review of an emerging area of research.Gait Posture161–14. 10.1016/s0966-6362(01)00156-4

  • CrossRef
  • Google Scholar

• 105

  ZabalaM. E.FavreJ.ScanlanS. F.DonahueJ.AndriacchiT. P. (2013). Three-dimensional knee moments of ACL reconstructed and control subjects during gait, stair ascent, and stair descent.J. Biomech.46515–520. 10.1016/j.jbiomech.2012.10.010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar