<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Sports Act. Living</journal-id>
<journal-title>Frontiers in Sports and Active Living</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Sports Act. Living</abbrev-journal-title>
<issn pub-type="epub">2624-9367</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fspor.2019.00025</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Sports and Active Living</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Phase-Dependency of Medial-Lateral Balance Responses to Sensory Perturbations During Walking</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Reimann</surname> <given-names>Hendrik</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/21780/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Fettrow</surname> <given-names>Tyler</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/552185/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Grenet</surname> <given-names>David</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/750326/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Thompson</surname> <given-names>Elizabeth D.</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/585569/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Jeka</surname> <given-names>John J.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/549989/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Kinesiology and Applied Physiology, University of Delaware</institution>, <addr-line>Newark, DE</addr-line>, <country>United States</country></aff>
<aff id="aff2"><sup>2</sup><institution>Department of Kinesiology, Temple University</institution>, <addr-line>Philadelphia, PA</addr-line>, <country>United States</country></aff>
<aff id="aff3"><sup>3</sup><institution>Department of Physical Therapy, University of Delaware</institution>, <addr-line>Newark, DE</addr-line>, <country>United States</country></aff>
<aff id="aff4"><sup>4</sup><institution>Department of Physical Therapy, Temple University</institution>, <addr-line>Philadelphia, PA</addr-line>, <country>United States</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Adamantios Arampatzis, Humboldt University of Berlin, Germany</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Vassilia Hatzitaki, Aristotle University of Thessaloniki, Greece; Christopher McCrum, Maastricht University, Netherlands</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Hendrik Reimann <email>hendrikreimann&#x00040;gmail.com</email></corresp>
<fn fn-type="other" id="fn001"><p>This article was submitted to Biomechanics and Control of Human Movement, a section of the journal Frontiers in Sports and Active Living</p></fn></author-notes>
<pub-date pub-type="epub">
<day>27</day>
<month>09</month>
<year>2019</year>
</pub-date>
<pub-date pub-type="collection">
<year>2019</year>
</pub-date>
<volume>1</volume>
<elocation-id>25</elocation-id>
<history>
<date date-type="received">
<day>04</day>
<month>07</month>
<year>2019</year>
</date>
<date date-type="accepted">
<day>28</day>
<month>08</month>
<year>2019</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2019 Reimann, Fettrow, Grenet, Thompson and Jeka.</copyright-statement>
<copyright-year>2019</copyright-year>
<copyright-holder>Reimann, Fettrow, Grenet, Thompson and Jeka</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract><p>The human body is mechanically unstable during walking. Maintaining upright stability requires constant regulation of muscle force by the central nervous system to push against the ground and move the body mass in the desired way. Activation of muscles in the lower body in response to sensory or mechanical perturbations during walking is usually highly phase-dependent, because the effect any specific muscle force has on the body movement depends upon the body configuration. Yet the resulting movement patterns of the upper body after the same perturbations are largely phase-independent. This is puzzling, because any change of upper-body movement must be generated by parts of the lower body pushing against the ground. How do phase-dependent muscle activation patterns along the lower body generate phase-independent movement patterns of the upper body? We hypothesize that when a sensory system detects a deviation of the body in space from a desired state that indicates the onset of a fall, the nervous system generates a functional response by pushing against the ground in any way possible with the current body configuration. This predicts that the changes in the ground reaction force patterns following a balance perturbation should be phase-independent. Here we test this hypothesis by disturbing upright balance in the frontal plane using Galvanic vestibular stimulation at three different points in the gait cycle. We measure the resulting changes in whole-body center of mass movement and the location of the center of pressure of the ground reaction force. We find that the magnitude of the initial center of pressure shift in the direction of the perceived fall is larger for perturbations late in the gait cycle, while there is no statistically significant difference in onset time. These results contradict our hypothesis by showing that even the initial CoP shift in response to a balance perturbation depends upon the phase of the gait cycle. Contrary to expectation, we also find that the whole-body balance response is not phase-independent. Both the onset time and the magnitude of the whole-body center of mass shift depend on the phase of the perturbation. We conclude that the central nervous system recruits any available mechanism to generate a functional balance response by pushing against the ground as fast as possible in response to a perturbation, but that the different mechanisms available at different phases in the gait cycle are not equally strong, leading to phase-dependent differences in the overall response.</p></abstract> <kwd-group>
<kwd>walking</kwd>
<kwd>balance</kwd>
<kwd>motor control</kwd>
<kwd>virtual reality</kwd>
<kwd>phase dependence</kwd>
<kwd>galvanic vestibular stimulation (GVS)</kwd>
</kwd-group>
<contract-sponsor id="cn001">Deutsche Forschungsgemeinschaft<named-content content-type="fundref-id">10.13039/501100001659</named-content></contract-sponsor>
<counts>
<fig-count count="5"/>
<table-count count="4"/>
<equation-count count="0"/>
<ref-count count="47"/>
<page-count count="11"/>
<word-count count="8413"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>1. Introduction</title>
<p>The control of balance during walking for humans is an important problem because the failure to maintain balance leads to falls and often injury. In the United States alone, costs from fall-related injuries amount to 20&#x02013;30 billion USD annually (Burns et al., <xref ref-type="bibr" rid="B5">2016</xref>). The upright human body is mechanically unstable, with a high center of mass and a relatively small base of support. Balancing this body is already a challenge during standing, requiring continuous regulation of muscle activity by the central nervous system (Morasso and Schieppati, <xref ref-type="bibr" rid="B26">1999</xref>). As a result of intensive research over the last decades, we understand this neural control of balance during standing reasonably well (e.g., Winter, <xref ref-type="bibr" rid="B45">1995</xref>; Peterka, <xref ref-type="bibr" rid="B31">2002</xref>; Maurer et al., <xref ref-type="bibr" rid="B25">2006</xref>; Kiemel et al., <xref ref-type="bibr" rid="B19">2011</xref>).</p>
<p>Balance control during walking has been studied much less in comparison. One reason for this is that walking is highly nonlinear, with the configuration of the body changing substantially as contact with the ground is established and lost at different points during the gait cycle. This implies that at different points in the gait cycle, the central nervous system has different options available for how to generate force along the body and against the ground to affect balance. One well-studied mechanism is to shift the location of the foot placement when taking a step, which changes the pull of gravity on the body during the subsequent swing (Hof, <xref ref-type="bibr" rid="B13">2008</xref>; Wang and Srinivasan, <xref ref-type="bibr" rid="B43">2014</xref>; Bruijn and van Die&#x000EB;n, <xref ref-type="bibr" rid="B4">2018</xref>; Reimann et al., <xref ref-type="bibr" rid="B33">2018a</xref>). Another mechanism is to use ankle musculature to pull on the body during single stance (Hof et al., <xref ref-type="bibr" rid="B16">2010</xref>; Hof and Duysens, <xref ref-type="bibr" rid="B14">2018</xref>; Reimann et al., <xref ref-type="bibr" rid="B35">2018b</xref>). Compared to foot placement shift, the ankle mechanism is limited in effect by the relatively small area of contact under the stance foot. Moments across the ankle joint both pull the body sideways and the foot up, so excessively large moments will result in the foot rolling over. On the other hand, the ankle mechanism has the advantage of being able to act much faster throughout single-stance and even double stance, whereas the foot placement mechanism can only be used when taking a step. For optimal benefit, there is preliminary evidence that humans flexibly coordinate these two balance mechanisms (Fettrow et al., <xref ref-type="bibr" rid="B11">2018</xref>), and possibly others, to adaptively respond to different challenges of balance during walking. A recent study by Vlutters et al. (<xref ref-type="bibr" rid="B42">2019</xref>) shows that the (sagittal) ankle mechanisms is activated even when the ankle joint is blocked by an orthosis, indicating that it might be controlled by a dedicated reflex.</p>
<p>Responses to balance-related perturbations generally depend upon the phase of the gait cycle in which they are applied.Reflexes in the lower leg, probed by cutaneous electric stimulation, are highly modulated during walking, and even reverse direction (Yang and Stein, <xref ref-type="bibr" rid="B46">1990</xref>; Zehr and Stein, <xref ref-type="bibr" rid="B47">1999</xref>). Changes of foot placement in response to electric stimulation of the vestibular system are also highly phase-dependent (Bent et al., <xref ref-type="bibr" rid="B3">2004</xref>). In contrast, responses in the upper body to these stimulations do not, or only weakly, depend on the phase of the perturbation. Logan et al. (<xref ref-type="bibr" rid="B23">2014</xref>) found similar results in response to visual perturbation, where the lower body response was strongly phase-dependent and the upper body response only weakly. But during normal walking, the only way to affect the lateral translational motion of the upper body is to generate a force against the ground, using muscles along the lower body. In combination, these results pose the question of how the central nervous system achieves a phase-independent response of the upper body using a phase-dependent response of the lower body.</p>
<p>We hypothesize that the central nervous system flexibly coordinates different balance mechanisms in the lower body in a phase-dependent way to stabilize the movement of the whole body in a phase-independent way. Here we test this hypothesis experimentally by perturbing humans walking on a treadmill with Galvanic vestibular stimulation at three different points in the gait cycle, and analyzing how the response in the whole body center of mass (CoM) and the functional force response depend upon the phase of the perturbation. We define the functional response as the displacement between the center of pressure (CoP) and the CoM. This variable is proportional to the acceleration of the CoM when assuming a single-link inverted pendulum model of the body biomechanics (Hof et al., <xref ref-type="bibr" rid="B15">2005</xref>), which has proven to be reasonable approximation in a variety of cases (Kuo, <xref ref-type="bibr" rid="B21">2007</xref>) and has been shown to be related to active balance control (Vlutters et al., <xref ref-type="bibr" rid="B41">2016</xref>). We expect that the functional balance response to these stimuli does not depend on the phase of the perturbation. We focus our analysis to the frontal plane, which is more challenging for balance control compared to the sagittal plane (O&#x00027;Connor and Kuo, <xref ref-type="bibr" rid="B29">2009</xref>).</p></sec>
<sec sec-type="methods" id="s2">
<title>2. Methods</title>
<p>Twenty young, healthy subjects (10 female) volunteered for this study. Subjects were between 19 and 38 years old (22.2 &#x000B1; 4.57), 170.7 &#x000B1; 8.4 cm tall and weighed 70.5 &#x000B1; 13.5 kg. Subjects provided informed verbal and written consent to participate. Subjects with self-reported history of neurological disorder or surgical procedures involving the legs, spine or head were excluded. The experimental design was approved by the Temple University Institutional Review Board (&#x00023;22499). We used a statistical power analysis based on pilot data from another experiment (Reimann et al., <xref ref-type="bibr" rid="B35">2018b</xref>) to determine the number of subjects required to reliably detect functionally relevant differences in the use of ankle roll and foot placement, set to 1 mms integrated CoP-CoM displacement and 5 mm foot placement change (1&#x02212;&#x003B2; &#x0003D; 0.9).</p>
<sec>
<title>2.1. Experimental Design</title>
<p>The study was conducted in the virtual reality setup of the Coordination of Balance and Locomotion laboratory at Temple University. Subjects walked on a treadmill in a virtual environment projected onto a curved screen (<xref ref-type="fig" rid="F1">Figure 1</xref>, Bertec, Inc.). This virtual reality environment is highly immersive, using a number of key components to mimic natural overground walking, listed below. The treadmill is self-paced, rather than operating at a fixed speed. This means that the treadmill belt speed dynamically adapts to the walking speed of the subject, and the subject can walk at their own pace, speeding up and slowing down as desired, without having to pay attention to their anterior-posterior position on the treadmill. The self-pacing is implemented using a nonlinear PD-controller in Labview (National instruments Inc., Austin, TX, USA) to keep the markers on the posterior superior iliac spine on the anterior-posterior mid-line of the treadmill. Walking patterns on a self-paced treadmill are generally very similar to walking on a fixed-speed treadmill (Sloot et al., <xref ref-type="bibr" rid="B36">2014</xref>). The virtual world is speed-linked to the treadmill. Subjects progress through the virtual world at a rate determined by the current speed of the self-paced treadmill. This re-establishes the proper connection between walking speed and optic flow experienced in natural overground walking that is broken when walking on a normal treadmill in a fixed environment. The view in the virtual world is perspective-linked to the subject&#x00027;s head. We map the real-time positions of the two anterior head markers from the motion capture system (see below) to the view point in the virtual world. This generates a motion parallax effect, allowing subjects to change perspective in the virtual world and look around objects by moving the head. The projection screen is domed and covers almost the complete field of vision. When looking generally ahead, the subject is only just able to make out the borders of the dome. This removes any visual anchor to the fixed laboratory environment, immersing the subject in the speed- and perspective-linked virtual environment. In combination, these component create an immersive environment with very few external visual clues or constraints. Subjects can in principle stop paying attention to the fact that they are on a treadmill.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Experimental setup. Subjects walked on a self-paced treadmill immersed in a virtual environment projected onto a curved screen that covered almost their complete field of vision. The perspective in the virtual world was linked to the subject&#x00027;s head position.</p></caption>
<graphic xlink:href="fspor-01-00025-g0001.tif"/>
</fig>
<p>The virtual environment consisted of a tiled marble floor with floating cubes randomly distributed in a volume 0&#x02013;10 m above the floor, 2&#x02013;17 m to each side from the midline, and infinitely into the distance, forming a 4 m wide corridor for the subjects to walk through (see <xref ref-type="fig" rid="F1">Figure 1</xref>), implemented in Unity3d (Unity Technologies, San Francisco, CA, USA). It was designed to have visual depth at mid-range to provide visual information about the movement of the body in space, but not at close range.</p>
<p>We used Galvanic vestibular stimulation (GVS) to induce the sensation of a fall to the side. Each stimulus was a square wave current with 500 &#x003BC;A amplitude and 600 ms duration between two round electrodes (3.2 cm diameter, Axelgaard Manufacturing 103 Co., Ltd, Fallbrook, CA, USA) fixed to the mastoid processes behind the ears. The amplitude was chosen to be large enough to induce a balance perturbation, but not so large as to substantially disrupt the gait pattern. Stimuli were triggered on heel-strikes of the right foot. Between the triggering heel-strike and the stimulus onset, we added a randomized delay of 0, 150 or 450 ms, which we will refer to as EARLY (0 ms), MID (150 ms) or LATE (450 ms) stimulus. The polarity of the stimulus was randomized to induce the sensation of falling either to the RIGHT or LEFT, where the direction of the fall sensation is toward the cathode of the electric current. A randomized wash-out period of 6&#x02013;8 strides followed each stimulus. Heel-strikes were defined as downward threshold crossings of the vertical heel-marker position, where the threshold was set to 3 mm above the vertical heel-marker position of each foot during quiet standing.</p>
<p>Forty-five reflective markers were placed on the subject, using the Plug-in Gait marker set (Davis et al., <xref ref-type="bibr" rid="B7">1991</xref>) with six additional markers on the anterior thigh, anterior tibia, and 5th metatarsal head of each foot. Marker positions were recorded at 250 Hz using an infrared motion capture system with 9 cameras (Vicon, Inc.). We collected surface electromyographical data bilaterally from five muscles along the legs and hips, but the data was corrupted due to problems with the measurement device (Delsys, Inc.), precluding analysis. Ground reaction forces and moments were collected at 1,000 Hz from both sides of the split-belt treadmill and transformed into a common coordinate frame to calculate whole-body CoP (Winter, <xref ref-type="bibr" rid="B44">1990</xref>). The same speed command was sent to each belt of the treadmill, effectively treating it as single-belt.</p>
<p>After explaining the experiment, obtaining consent and placing markers and EMG sensors, subjects first walked for 15 min on the self-paced treadmill in the virtual environment to adapt to this experimental setup. We then stopped the treadmill and exposed subjects to the GVS while standing to familiarize themselves with the sensation, then asked them to respond to this balance perturbation &#x0201C;normally&#x0201D; while walking on the treadmill. Data collection blocks consisted of two alternating phases for <italic>metronome</italic> and <italic>stimulus</italic>. During metronome phases, lasting 30 s, subjects were provided an auditory metronome at 90 beats per minute and asked to use this as an &#x0201C;approximate guideline&#x0201D; for their pace, both during metronome and stimulus phases. During stimulus phases, lasting 120 s, the metronome was turned off, and subjects received GVS as described above. Data were collected only during stimulus phases. Each subject performed four blocks of walking, each block consisting of five metronome and five stimulus phases, always starting with metronome phases, for a total of 12.5 min per block. After each block, the treadmill was stopped and subjects were offered a break. This protocol was implemented in a custom Labview program that sent the head position and treadmill speed to the Unity computer via UDP and the stimulation currency to the stimulator.</p>
<sec>
<title>Data Processing</title>
<p>We filled small gaps of up to 100 ms length in the kinematic data using cubic splines, then low pass filtered with a 4th order Butterworth filter at a cut-off frequency of 10 Hz. From the marker data, we calculated joint angle trajectories based on a geometric model with 15 segments (pelvis, torso, head, thighs, lower legs, feet, upper arms, forearms, hands) and 38&#x000B0; of freedom.We estimated the hip joint centers based upon pelvis landmarks (Tylkowski et al., <xref ref-type="bibr" rid="B39">1982</xref>; Bell et al., <xref ref-type="bibr" rid="B2">1990</xref>) and the knee joint centers and knee flexion rotational axes from reference movements using the symmetrical axis of rotation approach (Ehrig et al., <xref ref-type="bibr" rid="B9">2007</xref>). We performed inverse kinematics by minimizing the distance between the measured and the model-determined marker positions (Lu and O&#x00027;Connor, <xref ref-type="bibr" rid="B24">1999</xref>). This optimization was performed first for the six pelvis degrees of freedom, which formed the root of the kinematic tree, then for the 6&#x000B0; of freedom at the lumbar and cervical joints, and last for each of the arms and legs separately. We estimated the body center of mass (CoM) position based on estimated segment CoM locations (Dumas et al., <xref ref-type="bibr" rid="B8">2007</xref>) and the inverse kinematics and calculated CoM velocities and accelerations using numerical derivation by time. Force plate data was low pass filtered with a 4th order Butterworth filter at a cut-off frequency of 50 Hz.</p>
<p>We identified heelstrike events for each foot by finding minima in the vertical positions of the heel markers with inter-peak distances &#x0003E;250 ms and peak prominence &#x0003E;2 cm, and pushoff events as the first peak in the vertical velocity of the 2nd metatarsal marker with a prominence &#x0003E;0.35 ms<sup>&#x02212;1</sup> after each heelstrike. We visually inspected the result of this automatic identification and applied manual corrections in the rare cases where events were misidentified. We then partitioned the data into strides that started and ended with a right heel-strike. For each stimulus, we extracted data from the stride before and after the triggering right heel-strike, a total of two strides per trigger. The stride after the triggering heel-strike was analyzed as the stimulus data, and the stride before the trigger was used as the unperturbed reference. For each stride, we time-normalized the data for each step between consecutive heel-strikes. Strides containing missing kinematic data were excluded from further analysis. After removing these strides, an average of 43.2 &#x000B1; 6.8 stimulus strides remained for each subject.</p>
<p>For the <italic>whole-body response</italic>, we analyze the velocity and acceleration of the CoM. Although these are not independent of each other, we chose to include both, because the acceleration is more directly related to the forces applied by the muscles, but the velocity is less susceptible to processing artifacts from numerical derivation. We represent the <italic>functional response</italic> by the displacement between CoP and CoM, which is proportional to the CoM acceleration when approximating the body by a single-link linear inverted pendulum (Hof et al., <xref ref-type="bibr" rid="B15">2005</xref>). We characterized the <italic>lower-body response</italic> by changes of the swing foot placement relative to the stance foot at each step, and the <italic>upper-body response</italic> by the head and trunk roll angle. The roll angle was defined as the angle between the segment vector and the vertical in the frontal plane, where the head segment was the vector from the seventh cervical vertebra marker to the mid-point of the two markers on the back of the head, and the trunk segment was the vector from the mid-point of the two markers on the right and left posterior superior iliac spine to the seventh cervical vertebra marker. For each subject, we subtracted the mean of the unperturbed reference data from the stimulus data to estimate the response induced by the sensory perturbation. For visualization purposes, we also estimated the 95% confidence intervals for each trajectory across all repetitions, assuming no correlation from repeated measures within subjects. This estimate was only used to generate the shaded areas in the figures.</p></sec>
<sec>
<title>Estimation of Response Onset Time and Magnitude</title>
<p>Reliably estimating the onset time of the response is difficult due to the high amount of natural variability in walking. Our approach was to extract a short interval of data starting at the stimulus onset for each trigger, where we can assume that there is no change initially, but the response starts at some point during this interval. We approximated the time-dependency of each variable over this interval by fitting a piece-wise linear model in R (R Core Team, <xref ref-type="bibr" rid="B32">2013</xref>) with two segments and a variable break-point (Muggeo, <xref ref-type="bibr" rid="B28">2008</xref>), separately for each of the twenty subjects and six combinations of phase delay and stimulus direction. To estimate the response onset time, we used the location of the break-point between the two linear segments. To estimate the response magnitude, we used the absolute slope of the second linear segment. Note that this variable is the slope of a velocity response and has units of ms<sup>&#x02212;2</sup>, but since it is the slope of a linear model fit with specific constraints, we chose to not refer to it as acceleration.</p>
<p>The model fitting process was based on the following choices and assumptions. We assumed that there is no response initially, followed by a change of the outcome variable in a pre-determined direction after the onset time. This direction of change was expected to be the same as the fall stimulus direction for the CoP-CoM displacement, and the opposite for the CoM velocity and acceleration. To reflect these assumptions, we constrained the slope of the first linear segment to 0, and treated cases where the confidence interval for the break-point or slope included 0, or where the slope had the incorrect sign as missing data. The length of the analysis interval was 450 ms for the CoP-CoM displacement and 750 ms for the CoM velocity and acceleration. This value was chosen to be long enough to encompass the initial response to each stimulus, but short enough to avoid the more complex later modulations of each variable (see <xref ref-type="fig" rid="F2">Figures 2</xref>, <xref ref-type="fig" rid="F3">3</xref>).</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Responses in the movement pattern of the whole-body medial-lateral center of mass velocity <bold>(A,C,E)</bold> and acceleration <bold>(B,D,F)</bold>. Fall stimuli to the LEFT/RIGHT, corresponding to the side of the cathode, are shown in orange/blue tones. The left column <bold>(A,B)</bold> shows changes in the CoM kinematics following a perturbation. Arrows mark the approximate perturbation onset. Time is normalized, showing the two steps containing and following the perturbation, with double-stance periods shaded gray and single-stance periods white. The thick lines are average responses, the color-shaded areas are 95%-confidence intervals (see section 2). The center and right column show box plots of the estimated onset time <bold>(C,D)</bold> and magnitude <bold>(E,F)</bold> of the responses. Horizontal lines are medians, the boxes cover the first to third quartiles, whiskers the upper and lower adjacents, and the dots are single data points, where each point corresponds to one subject.</p></caption>
<graphic xlink:href="fspor-01-00025-g0002.tif"/>
</fig>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Responses in the medial-lateral center of pressure representing the hypothetical functional response to the balance perturbation. Fall stimuli to the LEFT/RIGHT, corresponding to the side of the cathode, are shown in orange/blue tones. <bold>(A)</bold> Shows changes in the CoP location following a perturbation. Arrows mark the approximate perturbation onset. Time is normalized, showing the two steps containing, and following the perturbation, with double-stance periods shaded gray and single-stance periods white. The thick lines are average responses, the color-shaded areas are 95%-confidence intervals (see section 2). <bold>(B,C)</bold> Show box plots of the estimated onset time <bold>(B)</bold> and magnitude <bold>(C)</bold> of the responses across subjects. Horizontal lines are medians, the boxes cover the first to third quartiles, whiskers the upper, and lower adjacents, and the dots are single data points, where each point corresponds to one subject.</p></caption>
<graphic xlink:href="fspor-01-00025-g0003.tif"/>
</fig>
<p>For the upper body variables head roll and trunk roll angle, the systematic response to the stimulus was so small relative to the natural variability from walking that the segmented model approach described above did not result in reliable estimates. For these two variables, we defined the response magnitude as the maximal excursion in the stimulus direction relative to the control average over four steps following each trigger. The lower body variable foot placement change is inherently discrete, so the notion of &#x0201C;onset time&#x0201D; is not meaningful. We defined the magnitude of the foot placement response as the change of the foot placement in the direction of the fall stimulus, i.e., the cathode.</p></sec>
<sec>
<title>Statistical Analysis</title>
<p>For each outcome variable, we performed a linear mixed effects analysis using R (R Core Team, <xref ref-type="bibr" rid="B32">2013</xref>) and <italic>lme4</italic> (Bates et al., <xref ref-type="bibr" rid="B1">2014</xref>), with fixed effects <italic>phase delay</italic> (EARLY, MID, LATE) and <italic>direction</italic> (LEFT, RIGHT), and random factor <italic>subject</italic>. We tested for significance of the fixed effects using an ANOVA with Satterthwaite&#x00027;s approximation method (Fai and Cornelius, <xref ref-type="bibr" rid="B10">1996</xref>), implemented in the R-package <italic>lmerTest</italic> (Kuznetsova et al., <xref ref-type="bibr" rid="B22">2017</xref>). We addressed multiple comparisons (ANOVAs for 10 different variables) by Bonferroni correction. Where the ANOVA reported a significant effect, we performed <italic>post-hoc</italic> pairwise <italic>t</italic>-tests with Bonferroni-adjusted <italic>p</italic>-values. In all tests, we used &#x003B1; &#x0003D; 0.05 as significance threshold.</p>
<p>We initially analyzed the effect of gender as a factor in all models. Gender had no significant effect on any outcome variable (<italic>p</italic> &#x0003E; 0.05). Comparing the model with gender to a model without gender as a factor, both the Akaike information criterion and the Bayesian information criterion were consistently lower for the model without gender. Based on these results, we used the more parsimonious models without gender as a factor for further analysis.</p></sec></sec></sec>
<sec sec-type="results" id="s3">
<title>3. Results</title>
<p>Subjects were able to successfully complete the experimental task of walking under intermittent Galvanic stimulation. There was no instance of a subject falling or utilizing the safety harness. Subjects walked with an average velocity of 0.97 &#x000B1; 0.15 m s<sup>&#x02212;1</sup> across the unperturbed steps. Subjects responded to the vestibular stimulation by swaying in the direction against the fall stimulus as expected (Reimann et al., <xref ref-type="bibr" rid="B34">2017</xref>).</p>
<sec>
<title>3.1. Whole Body Response</title>
<p><xref ref-type="fig" rid="F2">Figure 2</xref> shows the average CoM velocity and acceleration trajectories for the first two steps after the triggering heelstrike. Detailed results of the ANOVAs are reported in <xref ref-type="table" rid="T1">Table 1</xref>. For both CoM velocity and acceleration, <italic>phase delay</italic> had a significant effect on both the onset time and the magnitude of the response. The effect of <italic>stimulus direction</italic> was not significant for any upper body variable, nor was there any significant interaction between <italic>phase</italic> and <italic>direction</italic>. <xref ref-type="fig" rid="F2">Figures 2C,D</xref> illustrate these results. Both onset and magnitude of the response clearly depend upon <italic>phase</italic>, whereas the pattern for the two different stimulus directions is very similar. The <italic>post-hoc</italic> pair-wise <italic>t</italic>-tests for the onset of the response show that LATE is consistently different from EARLY and MID (<italic>p</italic> &#x0003C; 0.0017 for all comparisons), while MID and EARLY are significantly different from each other for CoM velocity (<italic>p</italic> &#x0003D; 0.0072), but not acceleration (<italic>p</italic> &#x0003D; 0.3230). For the magnitude of the response, LATE is consistently different from EARLY and MID (<italic>p</italic> &#x0003C; 0.0015 for all comparisons), but there is no significant difference between EARLY and MID (<italic>p</italic> &#x0003D; 1 for velocity, <italic>p</italic> &#x0003D; 0.9283 for acceleration).</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>Anova results for CoM velocity and acceleration.</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th/>
<th valign="top" align="center"><bold>Sum Sq</bold></th>
<th valign="top" align="center"><bold>Mean Sq</bold></th>
<th valign="top" align="center"><bold>NumDF</bold></th>
<th valign="top" align="center"><bold>DenDF</bold></th>
<th valign="top" align="center"><bold><italic>F</italic>-value</bold></th>
<th valign="top" align="center"><bold>Pr(&#x0003E;F)</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left" colspan="7"><bold>Onset time of the CoM velocity response</bold></td>
</tr>
<tr>
<td valign="top" align="left">Direction</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">1.00</td>
<td valign="top" align="center">86.57</td>
<td valign="top" align="center">0.10</td>
<td valign="top" align="center">0.7548</td>
</tr>
<tr>
<td valign="top" align="left">Delay</td>
<td valign="top" align="center">0.62</td>
<td valign="top" align="center">0.31</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">88.80</td>
<td valign="top" align="center">38.04</td>
<td valign="top" align="center">0.0000</td>
</tr>
<tr>
<td valign="top" align="left">Direction:delay</td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center">0.01</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">86.57</td>
<td valign="top" align="center">1.59</td>
<td valign="top" align="center">0.2101</td>
</tr>
<tr>
<td valign="top" align="left" colspan="7"><bold>Magnitude of the CoM velocity response</bold></td>
</tr>
<tr>
<td valign="top" align="left">Direction</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">1.00</td>
<td valign="top" align="center">85.58</td>
<td valign="top" align="center">4.06</td>
<td valign="top" align="center">0.0469</td>
</tr>
<tr>
<td valign="top" align="left">Delay</td>
<td valign="top" align="center">0.05</td>
<td valign="top" align="center">0.02</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">87.14</td>
<td valign="top" align="center">23.14</td>
<td valign="top" align="center">0.0000</td>
</tr>
<tr>
<td valign="top" align="left">Direction:delay</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">85.59</td>
<td valign="top" align="center">0.35</td>
<td valign="top" align="center">0.7037</td>
</tr>
<tr>
<td valign="top" align="left" colspan="7"><bold>Onset time of the CoM acceleration response</bold></td>
</tr>
<tr>
<td valign="top" align="left">Direction</td>
<td valign="top" align="center">0.02</td>
<td valign="top" align="center">0.02</td>
<td valign="top" align="center">1.00</td>
<td valign="top" align="center">73.49</td>
<td valign="top" align="center">1.47</td>
<td valign="top" align="center">0.2296</td>
</tr>
<tr>
<td valign="top" align="left">Delay</td>
<td valign="top" align="center">0.42</td>
<td valign="top" align="center">0.21</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">77.67</td>
<td valign="top" align="center">14.06</td>
<td valign="top" align="center">0.0000</td>
</tr>
<tr>
<td valign="top" align="left">Direction:delay</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">73.48</td>
<td valign="top" align="center">0.15</td>
<td valign="top" align="center">0.8605</td>
</tr>
<tr>
<td valign="top" align="left" colspan="7"><bold>Magnitude of the CoM acceleration response</bold></td>
</tr>
<tr>
<td valign="top" align="left">Direction</td>
<td valign="top" align="center">0.12</td>
<td valign="top" align="center">0.12</td>
<td valign="top" align="center">1.00</td>
<td valign="top" align="center">72.89</td>
<td valign="top" align="center">1.69</td>
<td valign="top" align="center">0.1978</td>
</tr>
<tr>
<td valign="top" align="left">Delay</td>
<td valign="top" align="center">1.69</td>
<td valign="top" align="center">0.85</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">77.80</td>
<td valign="top" align="center">12.16</td>
<td valign="top" align="center">0.0000</td>
</tr>
<tr>
<td valign="top" align="left">Direction:delay</td>
<td valign="top" align="center">0.02</td>
<td valign="top" align="center">0.01</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">72.86</td>
<td valign="top" align="center">0.17</td>
<td valign="top" align="center">0.8405</td>
</tr>
</tbody>
</table>
</table-wrap></sec>
<sec>
<title>3.2. Functional Response</title>
<p>The CoP relative to the CoM shifted in the direction of the perceived fall, shown in <xref ref-type="fig" rid="F3">Figure 3A</xref>. The onset time of this shift did not significantly depend upon the phase delay of the stimulus (see <xref ref-type="table" rid="T2">Table 2</xref>). <xref ref-type="fig" rid="F3">Figure 3B</xref> shows the estimated onset times. Visual inspection suggests a tendency for perturbations later in the step to have shorter onset times, but this was not significant. The magnitude of the CoP-CoM shift depended significantly upon the phase delay (see <xref ref-type="table" rid="T2">Table 2</xref>). <xref ref-type="fig" rid="F3">Figure 3C</xref> shows that the magnitude of the response was larger for LATE than for EARLY and MID for both stimulus directions, and this increase is statistically significant (<italic>post-hoc</italic> pairwise <italic>t</italic>-test, <italic>p</italic> &#x0003C; 0.0024 for both comparisons). The magnitude for MID is smaller from EARLY for both stimulus directions, but this difference is not statistically significant (<italic>p</italic> &#x0003D; 0.2135).</p>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption><p>Anova results for onset time and magnitude of the function response (the CoP-CoM displacement).</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th/>
<th valign="top" align="center"><bold>Sum Sq</bold></th>
<th valign="top" align="center"><bold>Mean Sq</bold></th>
<th valign="top" align="center"><bold>NumDF</bold></th>
<th valign="top" align="center"><bold>DenDF</bold></th>
<th valign="top" align="center"><bold><italic>F</italic>-value</bold></th>
<th valign="top" align="center"><bold>Pr(&#x0003E;F)</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left" colspan="7"><bold>Onset time of the functional response</bold></td>
</tr>
<tr>
<td valign="top" align="left">Direction</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">1.00</td>
<td valign="top" align="center">70.04</td>
<td valign="top" align="center">1.04</td>
<td valign="top" align="center">0.3103</td>
</tr>
<tr>
<td valign="top" align="left">Delay</td>
<td valign="top" align="center">0.02</td>
<td valign="top" align="center">0.01</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">72.34</td>
<td valign="top" align="center">1.82</td>
<td valign="top" align="center">0.1692</td>
</tr>
<tr>
<td valign="top" align="left">Direction:delay</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">68.66</td>
<td valign="top" align="center">0.55</td>
<td valign="top" align="center">0.5789</td>
</tr>
<tr>
<td valign="top" align="left" colspan="7"><bold>Magnitude of the functional response</bold></td>
</tr>
<tr>
<td valign="top" align="left">Direction</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">1.00</td>
<td valign="top" align="center">67.79</td>
<td valign="top" align="center">0.91</td>
<td valign="top" align="center">0.3433</td>
</tr>
<tr>
<td valign="top" align="left">Delay</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">70.86</td>
<td valign="top" align="center">14.11</td>
<td valign="top" align="center">0.0000</td>
</tr>
<tr>
<td valign="top" align="left">Direction:delay</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">66.30</td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center">0.9738</td>
</tr>
</tbody>
</table>
</table-wrap></sec>
<sec>
<title>3.3. Upper Body</title>
<p>Both the head and the trunk segments leaned in the direction against the perceived fall in response to the stimulus. <xref ref-type="fig" rid="F4">Figures 4A,B</xref> show the average trajectories of the head and trunk roll angle for the first four steps after the triggering heelstrike, and <xref ref-type="fig" rid="F4">Figures 4C,D</xref> show box-plots for the magnitude of the response. For the head angle, the magnitude of the roll response does not depend upon the direction, nor the phase delay of the stimulus, in a statistically significant way (see <xref ref-type="table" rid="T3">Table 3</xref>). For the trunk angle, both phase delay and stimulus direction have a statistically significant effect (see <xref ref-type="table" rid="T3">Table 3</xref>). The <italic>post-hoc</italic> pair-wise <italic>t</italic>-test shows that trunk roll magnitude is significantly different between all combinations of delays (<italic>p</italic> &#x0003C; 0.0037 for all comparisons).</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Responses in the roll angles of the head <bold>(A,C)</bold> and trunk <bold>(B,D)</bold> segments. Fall stimuli to the LEFT/RIGHT, corresponding to the side of the cathode, are shown in orange/blue tones. <bold>(A,B)</bold> Show changes in the roll angles following a perturbation. Arrows mark the approximate perturbation onset. Time is normalized, showing the four steps containing and following the perturbation, with double-stance periods shaded gray and single-stance periods white (note that the time shown here is twice as long as in <xref ref-type="fig" rid="F2">Figures 2</xref>, <xref ref-type="fig" rid="F3">3</xref>, <xref ref-type="fig" rid="F5">5</xref>). The thick lines are average responses, the color-shaded areas are 95%-confidence intervals (see section 2). <bold>(C,D)</bold> Show box plots of the estimated magnitude of the responses. Horizontal lines are medians, the boxes cover the first to third quartiles, whiskers the upper and lower adjacents, and the dots are single data points, where each point corresponds to one step.</p></caption>
<graphic xlink:href="fspor-01-00025-g0004.tif"/>
</fig>
<table-wrap position="float" id="T3">
<label>Table 3</label>
<caption><p>Anova results for magnitude of the response in head and trunk roll angle.</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th/>
<th valign="top" align="center"><bold>Sum Sq</bold></th>
<th valign="top" align="center"><bold>Mean Sq</bold></th>
<th valign="top" align="center"><bold>NumDF</bold></th>
<th valign="top" align="center"><bold>DenDF</bold></th>
<th valign="top" align="center"><bold><italic>F</italic>-value</bold></th>
<th valign="top" align="center"><bold>Pr(&#x0003E;F)</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left" colspan="7"><bold>Magnitude of the head roll response</bold></td>
</tr>
<tr>
<td valign="top" align="left">Direction</td>
<td valign="top" align="center">1.41</td>
<td valign="top" align="center">1.41</td>
<td valign="top" align="center">1.00</td>
<td valign="top" align="center">5150.82</td>
<td valign="top" align="center">0.39</td>
<td valign="top" align="center">0.5330</td>
</tr>
<tr>
<td valign="top" align="left">Delay</td>
<td valign="top" align="center">31.92</td>
<td valign="top" align="center">15.96</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">5149.00</td>
<td valign="top" align="center">4.39</td>
<td valign="top" align="center">0.0125</td>
</tr>
<tr>
<td valign="top" align="left">Direction:delay</td>
<td valign="top" align="center">3.66</td>
<td valign="top" align="center">1.83</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">5149.06</td>
<td valign="top" align="center">0.50</td>
<td valign="top" align="center">0.6047</td>
</tr>
<tr>
<td valign="top" align="left" colspan="7"><bold>Magnitude of the trunk roll response</bold></td>
</tr>
<tr>
<td valign="top" align="left">Direction</td>
<td valign="top" align="center">224.35</td>
<td valign="top" align="center">224.35</td>
<td valign="top" align="center">1.00</td>
<td valign="top" align="center">5149.76</td>
<td valign="top" align="center">138.60</td>
<td valign="top" align="center">0.0000</td>
</tr>
<tr>
<td valign="top" align="left">Delay</td>
<td valign="top" align="center">106.06</td>
<td valign="top" align="center">53.03</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">5148.61</td>
<td valign="top" align="center">32.76</td>
<td valign="top" align="center">0.0000</td>
</tr>
<tr>
<td valign="top" align="left">Direction:delay</td>
<td valign="top" align="center">338.75</td>
<td valign="top" align="center">169.37</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">5148.64</td>
<td valign="top" align="center">104.63</td>
<td valign="top" align="center">0.0000</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec>
<title>3.4. Foot Placement</title>
<p>The foot placement at the fist post-stimulus step was shifted in the direction of the perceived fall, but only for the EARLY and MID phase delays (see <xref ref-type="fig" rid="F5">Figure 5</xref>). At the second post-stimulus step, the foot placement was strongly shifted in the opposite direction for EARLY and MID. In the LATE condition, neither the first nor the second foot placement seems to have shifted. The left side of <xref ref-type="fig" rid="F5">Figure 5</xref> show the trajectories of both heel markers for the first two steps following the trigger, and the right side shows box plots for the magnitude of the foot placement change. <xref ref-type="table" rid="T4">Table 4</xref> reports the results of the ANOVAs for magnitude of the foot placement change. Phase delay has a significant effect on both the first and the second step. Stimulus direction is not significant on the first step, but is significant on the second step. There is no interaction between delay and direction. The <italic>post-hoc</italic> pairwise <italic>t</italic>-tests show that foot placement is significantly different between all combinations of delays on both post-stimulus steps (<italic>p</italic> &#x0003C; 0.0001 for all comparisons).</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>Responses of the medial-lateral foot kinematics. Fall stimuli to the LEFT/RIGHT, corresponding to the side of the cathode, are shown in orange/blue tones. <bold>(A,B)</bold> Show changes in the average medial-lateral heel position for both legs. Arrows mark the approximate perturbation onset. Time is normalized, showing the two steps containing, and following the perturbation, with double-stance periods shaded gray, and single-stance periods white. The thick lines are average responses, the color-shaded areas are 95%-confidence intervals (see section 2). The right column shows box plots of the estimated and magnitude of the foot placement response, i.e., the medial-lateral position of the leading foot heel relative to the trailing foot heel at each heel-strike. Step one <bold>(C)</bold> corresponds to the time in the middle of the trajectories on the left. Step two <bold>(D)</bold> corresponds to the time at the end of the trajectories on the left. Horizontal lines are medians, the boxes cover the first to third quartiles, whiskers the upper and lower adjacents, and the dots are single data points, each point corresponds to one step.</p></caption>
<graphic xlink:href="fspor-01-00025-g0005.tif"/>
</fig>
<table-wrap position="float" id="T4">
<label>Table 4</label>
<caption><p>Anova results for magnitude of the foot placement response the two post-stimulus steps.</p></caption>
<table frame="hsides" rules="groups">
<thead><tr>
<th/>
<th valign="top" align="center"><bold>Sum Sq</bold></th>
<th valign="top" align="center"><bold>Mean Sq</bold></th>
<th valign="top" align="center"><bold>NumDF</bold></th>
<th valign="top" align="center"><bold>DenDF</bold></th>
<th valign="top" align="center"><bold><italic>F</italic>-value</bold></th>
<th valign="top" align="center"><bold>Pr(&#x0003E;F)</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left" colspan="7"><bold>Magnitude of the foot placement response in step ONE</bold></td>
</tr>
<tr>
<td valign="top" align="left">Direction</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">1.00</td>
<td valign="top" align="center">5177.21</td>
<td valign="top" align="center">0.20</td>
<td valign="top" align="center">0.6551</td>
</tr>
<tr>
<td valign="top" align="left">Delay</td>
<td valign="top" align="center">0.24</td>
<td valign="top" align="center">0.12</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">5168.41</td>
<td valign="top" align="center">100.97</td>
<td valign="top" align="center">0.0000</td>
</tr>
<tr>
<td valign="top" align="left">Direction:delay</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">5168.39</td>
<td valign="top" align="center">1.00</td>
<td valign="top" align="center">0.3672</td>
</tr>
<tr>
<td valign="top" align="left" colspan="7"><bold>Magnitude of the foot placement response in step TWO</bold></td>
</tr>
<tr>
<td valign="top" align="left">Direction</td>
<td valign="top" align="center">0.08</td>
<td valign="top" align="center">0.08</td>
<td valign="top" align="center">1.00</td>
<td valign="top" align="center">5162.46</td>
<td valign="top" align="center">46.39</td>
<td valign="top" align="center">0.0000</td>
</tr>
<tr>
<td valign="top" align="left">Delay</td>
<td valign="top" align="center">7.20</td>
<td valign="top" align="center">3.60</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">5161.56</td>
<td valign="top" align="center">2013.28</td>
<td valign="top" align="center">0.0000</td>
</tr>
<tr>
<td valign="top" align="left">Direction:delay</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">0.00</td>
<td valign="top" align="center">2.00</td>
<td valign="top" align="center">5161.56</td>
<td valign="top" align="center">1.33</td>
<td valign="top" align="center">0.2633</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec></sec>
<sec sec-type="discussion" id="s4">
<title>4. Discussion</title>
<p>We used Galvanic vestibular stimulation to perturb the sense of balance in walking humans at three different phases in the gait cycle. We analyzed the whole-body response to this sensory perturbation in the form of changes of the center of mass (CoM) movement, and the displacement of the center of pressure (CoP) relative to the CoM, which is one aspect of the ground reaction forces generating the whole-body movement changes. Based on previous studies, we expected that the whole-body response would be phase-independent, and hypothesized that the functional response of the CoP-CoM displacement would also be phase-independent.</p>
<p>Our results broadly support and extend the finding that vestibular information is used to maintain balance during walking and that the balance response to vestibular stimuli depends on the phase of the gait cycle (Bent et al., <xref ref-type="bibr" rid="B3">2004</xref>; Iles et al., <xref ref-type="bibr" rid="B18">2007</xref>; Dakin et al., <xref ref-type="bibr" rid="B6">2013</xref>). Bent et al. (<xref ref-type="bibr" rid="B3">2004</xref>) reported that upper-body kinematic responses to GVS were independent of the phase of the gait cycle but that the foot-placement response was phase dependent. Dakin et al. (<xref ref-type="bibr" rid="B6">2013</xref>) found that there is vestibular influence on motor output across many muscles throughout the gait cycle, but that the contribution of any individual muscle to the overall response is highly phase-dependent.</p>
<p>Subjects responded to the fall stimulus by shifting their CoP in the direction of the perceived fall, i.e., the cathode, which resulted in a body sway in the opposite direction, toward the anode, as expected. The onset time of this functional response did not depend on the phase of the perturbation (see <xref ref-type="fig" rid="F3">Figure 3</xref> and <xref ref-type="table" rid="T2">Table 2</xref>). This indicates that although different balance mechanisms are recruited at different points in the gait cycle, the neural controller tends to respond as fast as possible with any available mechanism. For the EARLY and MID perturbations, the foot placement following the perturbation was shifted in the direction of the perceived fall, but the first CoP shift occurred during single stance, indicating that the controller recruits the immediately available mechanism first and does not wait until other, potentially stronger mechanisms become available.</p>
<p>Contrary to our expectations and to results by Bent et al. (<xref ref-type="bibr" rid="B3">2004</xref>), the onset time of the whole-body response did depend on the phase of the perturbation. Responses to the LATE perturbations were consistently faster than to the EARLY and MID perturbations, with no significant difference between the latter two conditions (see <xref ref-type="fig" rid="F2">Figure 2</xref> and <xref ref-type="table" rid="T1">Table 1</xref>).</p>
<p>The magnitude of both the initial CoP shift in the direction of the perceived fall and the whole-body sway in the opposite direction depended upon the phase of the perturbation. The whole-body response was consistently smaller for the LATE perturbation compared to EARLY and MID (see <xref ref-type="fig" rid="F2">Figure 2</xref> and <xref ref-type="table" rid="T1">Table 1</xref>). The functional response, in contrast, was consistently larger for the LATE perturbations (see <xref ref-type="fig" rid="F3">Figure 3</xref>). This is not consistent with Bent et al. (<xref ref-type="bibr" rid="B3">2004</xref>), who found that the size of the upper-body response did not depend on the timing of the stimulus. This may be because we used a smaller GVS stimulus than (0.5 mA vs. 1&#x02013;1.5 mA, Bent et al., <xref ref-type="bibr" rid="B3">2004</xref>) and also different stimulus timing, with our LATE stimulus falling in between their mid-stance and toe-off stimuli.</p>
<p>On the whole, there were few significant differences between EARLY and MID perturbations, but both of these were consistently significantly different from the LATE perturbations. The CoM responses were faster and smaller, whereas the CoP shift responses were larger, but with unchanged onset time. In other words, the onset time of the functional response is phase-independent, but that of the whole-body response is not, and the magnitude of both responses is phase-dependent, but in opposite directions. This difference in response pattern is peculiar, because as noted earlier, all changes of whole-body movement have to be generated by pushing against the ground using the lower body. So how can a smaller CoP shift generate a larger CoM movement change?</p>
<p>One possible explanation is that these differences in the whole-body movement are generated by combinations of forces that do not affect the location of the CoP, so the proportional relationship between CoP and CoM movement based on the single-link inverted pendulum assumption does not hold. One example of such forces is a hip-roll mechanism, where the upper body is rotated around the hip, generating a shear force that pushes the CoM sideways without shifting the CoP (Horak and Nashner, <xref ref-type="bibr" rid="B17">1986</xref>; Reimann et al., <xref ref-type="bibr" rid="B33">2018a</xref>). While this hip mechanism is usually discussed in standing or the single-stance period during walking, a biomechanical equivalent exists during double stance. This would involve differences in the trunk roll angle response early after the stimulus. <xref ref-type="fig" rid="F4">Figure 4</xref> shows that there are indeed differences, but they are neither large nor conclusive. Explaining such details satisfactorily generally requires a detailed model of the whole sensorimotor control loop, including neural dynamics, muscle physiology, biomechanics, and interaction with the environment. While such models exist for standing (e.g., van der Kooij et al., <xref ref-type="bibr" rid="B40">1999</xref>; Peterka, <xref ref-type="bibr" rid="B30">2000</xref>), there is no comparable model to explain balance control during walking. The available models of walking are not sufficiently detailed to explain the results presented here (Taga, <xref ref-type="bibr" rid="B38">1998</xref>; Geyer and Herr, <xref ref-type="bibr" rid="B12">2010</xref>; Song and Geyer, <xref ref-type="bibr" rid="B37">2015</xref>).</p>
<p>We hypothesized that the CNS generates a functional, phase-independent motor response to balance perturbations by recruiting different balance mechanisms in a phase-dependent way. Based on this hypothesis, we predicted that the CoP shift relative to the CoM in response to a sensory fall stimulus would not depend on the phase-shift in the gait cycle in which that stimulus occurred. This prediction was partially supported by the data, in that the onset of the response was phase-independent, but also partially contradicted, in that the magnitude was phase-dependent. It is important to point out that lack of statistical significance of <italic>phase</italic> for the functional response is not sufficient evidence of phase-independence. It is possible that there is a phase-dependent difference, but this study was not sufficiently powered to detect it. Indeed, visual inspection of <xref ref-type="fig" rid="F3">Figure 3B</xref> indicates a similar onset time pattern between LEFT and RIGHT, but more data is necessary to tell if this is an actual effect.</p>
<p>Our results indicate that the CNS recruits available balance mechanisms as fast as possible to respond to a sensed threat to upright stability, but the strength of the response is initially limited during single stance and only realized after a period when other mechanisms also contribute. The reason for this limitation is likely biomechanical. Responses to the EARLY and MID perturbations begin during single stance, where only the ankle roll mechanism is available. The degree to which ankle roll can shift the CoP is limited by the surface size of the foot sole (Hof et al., <xref ref-type="bibr" rid="B16">2010</xref>), which is a possible explanation for the reduced magnitude of the CoP shift in EARLY and MID compared to LATE perturbations.</p>
<p>All stimuli were triggered on the right heel-strike, but the direction of the fall stimulus, determined by the polarity of the GVS current, was randomized. The direction of the stimulus had relatively little effect. Statistical analysis resulted in significance of <italic>direction</italic> after Bonferroni correction only for the magnitude of the trunk roll response. Visual inspection of the results supports this. For the whole body response (<xref ref-type="fig" rid="F2">Figure 2</xref>), functional response (<xref ref-type="fig" rid="F3">Figure 3</xref>) and the lower body response (<xref ref-type="fig" rid="F5">Figure 5</xref>, the patterns of both onset time and magnitude of the response are very similar between LEFT and RIGHT stimuli. Only magnitude of the trunk roll angle response (<xref ref-type="fig" rid="F4">Figure 4</xref>) in the LATE condition seems to be a distinct departure from this pattern. We observed a similar lack of dependence on direction in our previous study using visual perturbations (Reimann et al., <xref ref-type="bibr" rid="B35">2018b</xref>). This lack of effect is slightly surprising, since one might expect that a perceived fall toward the stance foot is less threatening to overall stability than away from the stance foot.</p>
<p>Such a limitation of balance responses during single stance might be the underlying reason that drives some populations with balance deficits to adapt their gait patterns in characteristic ways. For example, people with Parkinson&#x00027;s disease tend to have shorter single-stance times than age-matched controls for similar walking speed (Morris et al., <xref ref-type="bibr" rid="B27">1994</xref>). This might be because they have problems using the ankle roll mechanism for balance control and reduce the duration of the single stance period as a coping strategy.</p>
<p>The details of which balance mechanisms are recruited at which point in the gait cycle are still not well understood. One striking observation from our data is that for the LATE perturbation, the onset of the average CoP shift appears to align with the transition to double stance (<xref ref-type="fig" rid="F3">Figure 3</xref>). But since there is no change in foot placement (<xref ref-type="fig" rid="F5">Figure 5</xref>), this CoP shift cannot be generated by the well-understood foot placement mechanism. One possibility is that the transition to double stance allows recruitment of the push-off mechanism to shift weight between the two stance legs by modulating the push-off force of the trailing limb around that transition (Kim and Collins, <xref ref-type="bibr" rid="B20">2015</xref>; Reimann et al., <xref ref-type="bibr" rid="B35">2018b</xref>). Understanding this phenomenon requires a detailed kinematic analysis that we will perform in future work.</p></sec>
<sec sec-type="data-availability" id="s5">
<title>Data Availability</title>
<p>The datasets generated for this study are available on request to the corresponding author.</p></sec>
<sec id="s6">
<title>Ethics Statement</title>
<p>This studies involving human participants were reviewed and approved by University of Delaware IRB. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.</p></sec>
<sec id="s7">
<title>Author Contributions</title>
<p>HR, TF, and JJ designed the experiment. HR, TF, and ET performed the experiment. HR, TF, DG, ET, and JJ analyzed the data. HR, TF, DG, and JJ wrote the manuscript.</p>
<sec>
<title>Conflict of Interest Statement</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p></sec></sec>
</body>
<back>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bates</surname> <given-names>D.</given-names></name> <name><surname>Machler</surname> <given-names>M.</given-names></name> <name><surname>Bolker</surname> <given-names>B.</given-names></name> <name><surname>Walker</surname> <given-names>S.</given-names></name></person-group> (<year>2014</year>). <article-title>Fitting linear mixed-effects models using lme4</article-title>. <source>arXiv</source>. arXiv: 1406.5823. <pub-id pub-id-type="doi">10.18637/jss.v067.i01</pub-id></citation></ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bell</surname> <given-names>L.</given-names></name> <name><surname>Pedersen</surname> <given-names>R.</given-names></name> <name><surname>Brand</surname> <given-names>A.</given-names></name></person-group> (<year>1990</year>). <article-title>A comparison of the accuracy of several hip center location prediction methods</article-title>. <source>J. Biomechan.</source> <volume>23</volume>, <fpage>617</fpage>&#x02013;<lpage>621</lpage>. <pub-id pub-id-type="pmid">2341423</pub-id></citation></ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bent</surname> <given-names>L. R.</given-names></name> <name><surname>Inglis</surname> <given-names>J. T.</given-names></name> <name><surname>McFadyen</surname> <given-names>B. J.</given-names></name></person-group> (<year>2004</year>). <article-title>When is vestibular information important during walking?</article-title> <source>J. Neurophysiol.</source> <volume>92</volume>, <fpage>1269</fpage>&#x02013;<lpage>1275</lpage>. <pub-id pub-id-type="doi">10.1152/jn.01260.2003</pub-id><pub-id pub-id-type="pmid">15102904</pub-id></citation></ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bruijn</surname> <given-names>S. M.</given-names></name> <name><surname>van Die&#x000EB;n</surname> <given-names>J. H.</given-names></name></person-group> (<year>2018</year>). <article-title>Control of human gait stability through foot placement</article-title>. <source>J. Roy. Soc. Interface</source> <volume>15</volume>:<fpage>20170816</fpage>. <pub-id pub-id-type="doi">10.1098/rsif.2017.0816</pub-id><pub-id pub-id-type="pmid">29875279</pub-id></citation></ref>
<ref id="B5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Burns</surname> <given-names>E. R.</given-names></name> <name><surname>Stevens</surname> <given-names>J. A.</given-names></name> <name><surname>Lee</surname> <given-names>R.</given-names></name></person-group> (<year>2016</year>). <article-title>The direct costs of fatal and non-fatal falls among older adults - United States</article-title>. <source>J. Safe. Res.</source> <volume>58</volume>, <fpage>99</fpage>&#x02013;<lpage>103</lpage>. <pub-id pub-id-type="doi">10.1016/j.jsr.2016.05.001</pub-id><pub-id pub-id-type="pmid">27620939</pub-id></citation></ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dakin</surname> <given-names>C. J.</given-names></name> <name><surname>Inglis</surname> <given-names>J. T.</given-names></name> <name><surname>Chua</surname> <given-names>R.</given-names></name> <name><surname>Blouin</surname> <given-names>J.-S.</given-names></name></person-group> (<year>2013</year>). <article-title>Muscle-specific modulation of vestibular reflexes with increased locomotor velocity and cadence</article-title>. <source>J. Neurophysiol.</source> <volume>110</volume>, <fpage>86</fpage>&#x02013;<lpage>94</lpage>. <pub-id pub-id-type="doi">10.1152/jn.00843.2012</pub-id><pub-id pub-id-type="pmid">23576695</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Davis</surname> <given-names>R. B.</given-names></name> <name><surname>Ounpuu</surname> <given-names>S.</given-names></name> <name><surname>Tyburski</surname> <given-names>D.</given-names></name> <name><surname>Gage</surname> <given-names>J. R.</given-names></name></person-group> (<year>1991</year>). <article-title>A gait analysis data collection and reduction technique</article-title>. <source>Hum. Move. Sci.</source> <volume>10</volume>, <fpage>575</fpage>&#x02013;<lpage>587</lpage>. <pub-id pub-id-type="doi">10.1016/0167-9457(91)90046-Z</pub-id></citation></ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dumas</surname> <given-names>R.</given-names></name> <name><surname>Ch&#x000E8;ze</surname> <given-names>L.</given-names></name> <name><surname>Verriest</surname> <given-names>J.-P.</given-names></name></person-group> (<year>2007</year>). <article-title>Adjustments to McConville et al., and Young et al. body segment inertial parameters</article-title>. <source>J. Biomechan.</source> <volume>40</volume>, <fpage>543</fpage>&#x02013;<lpage>553</lpage>. <pub-id pub-id-type="doi">10.1016/j.jbiomech.2006.02.013</pub-id><pub-id pub-id-type="pmid">16616757</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ehrig</surname> <given-names>R. M.</given-names></name> <name><surname>Taylor</surname> <given-names>W. R.</given-names></name> <name><surname>Duda</surname> <given-names>G. N.</given-names></name> <name><surname>Heller</surname> <given-names>M. O.</given-names></name></person-group> (<year>2007</year>). <article-title>A survey of formal methods for determining functional joint axes</article-title>. <source>J. Biomechan.</source> <volume>40</volume>, <fpage>2150</fpage>&#x02013;<lpage>2157</lpage>. <pub-id pub-id-type="doi">10.1016/j.jbiomech.2006.10.026</pub-id><pub-id pub-id-type="pmid">17169365</pub-id></citation></ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fai</surname> <given-names>A. H.-T.</given-names></name> <name><surname>Cornelius</surname> <given-names>P. L.</given-names></name></person-group> (<year>1996</year>). <article-title>Approximate F-tests of multiple degree of freedom hypotheses in generalized least squares analyses of unbalanced split-plot experiments</article-title>. <source>J. Stat. Comput. Simulat.</source> <volume>54</volume>, <fpage>363</fpage>&#x02013;<lpage>378</lpage>.</citation></ref>
<ref id="B11">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Fettrow</surname> <given-names>T. D.</given-names></name> <name><surname>Reimann</surname> <given-names>H.</given-names></name> <name><surname>Thompson</surname> <given-names>E. D.</given-names></name> <name><surname>Grenet</surname> <given-names>D.</given-names></name> <name><surname>Jeka</surname> <given-names>J.</given-names></name></person-group> (<year>2018</year>). <article-title>Interdependence of balance mechanisms during walking</article-title>, in <source>Poster Presentation at Society for Neuroscience Meeting</source> (<publisher-loc>San Diego, CA</publisher-loc>).</citation></ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Geyer</surname> <given-names>H.</given-names></name> <name><surname>Herr</surname> <given-names>H.</given-names></name></person-group> (<year>2010</year>). <article-title>A Muscle-reflex model that encodes principles of legged mechanics produces human walking dynamics and muscle activities</article-title>. <source>IEEE Trans, Neural Syst. Rehabil. Eng.</source> <volume>18</volume>, <fpage>263</fpage>&#x02013;<lpage>273</lpage>. <pub-id pub-id-type="doi">10.1109/TNSRE.2010.2047592</pub-id><pub-id pub-id-type="pmid">20378480</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hof</surname> <given-names>A. L.</given-names></name></person-group> (<year>2008</year>). <article-title>The &#x00027;extrapolated center of mass&#x00027; concept suggests a simple control of balance in walking</article-title>. <source>Hum. Move. Sci.</source> <volume>27</volume>, <fpage>112</fpage>&#x02013;<lpage>125</lpage>. <pub-id pub-id-type="doi">10.1016/j.humov.2007.08.003</pub-id><pub-id pub-id-type="pmid">17935808</pub-id></citation></ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hof</surname> <given-names>A. L.</given-names></name> <name><surname>Duysens</surname> <given-names>J.</given-names></name></person-group> (<year>2018</year>). <article-title>Responses of human ankle muscles to mediolateral balance perturbations during walking</article-title>. <source>Hum. Move. Sci.</source> <volume>57</volume>, <fpage>69</fpage>&#x02013;<lpage>82</lpage>. <pub-id pub-id-type="doi">10.1016/j.humov.2017.11.009</pub-id><pub-id pub-id-type="pmid">29174418</pub-id></citation></ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hof</surname> <given-names>A. L.</given-names></name> <name><surname>Gazendam</surname> <given-names>M. G. J.</given-names></name> <name><surname>Sinke</surname> <given-names>W. E.</given-names></name></person-group> (<year>2005</year>). <article-title>The condition for dynamic stability</article-title>. <source>J. Biomechan.</source> <volume>38</volume>, <fpage>1</fpage>&#x02013;<lpage>8</lpage>. <pub-id pub-id-type="doi">10.1016/j.jbiomech.2004.03.025</pub-id><pub-id pub-id-type="pmid">15519333</pub-id></citation></ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hof</surname> <given-names>A. L.</given-names></name> <name><surname>Vermerris</surname> <given-names>S. M.</given-names></name> <name><surname>Gjaltema</surname> <given-names>W. A.</given-names></name></person-group> (<year>2010</year>). <article-title>Balance responses to lateral perturbations in human treadmill walking</article-title>. <source>J. Exp. Biol.</source> <volume>213</volume>(<issue>Pt 15</issue>), <fpage>2655</fpage>&#x02013;<lpage>2664</lpage>. <pub-id pub-id-type="doi">10.1242/jeb.042572</pub-id><pub-id pub-id-type="pmid">20639427</pub-id></citation></ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Horak</surname> <given-names>F. B.</given-names></name> <name><surname>Nashner</surname> <given-names>L. M.</given-names></name></person-group> (<year>1986</year>). <article-title>Central programming of postural movements: adaptation to altered support-surface configurations</article-title>. <source>J. Neurophysiol.</source> <volume>55</volume>, <fpage>1369</fpage>&#x02013;<lpage>1381</lpage>. <pub-id pub-id-type="doi">10.1152/jn.1986.55.6.1369</pub-id><pub-id pub-id-type="pmid">3734861</pub-id></citation></ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Iles</surname> <given-names>J. F.</given-names></name> <name><surname>Baderin</surname> <given-names>R.</given-names></name> <name><surname>Tanner</surname> <given-names>R.</given-names></name> <name><surname>Simon</surname> <given-names>A.</given-names></name></person-group> (<year>2007</year>). <article-title>Human standing and walking: comparison of the effects of stimulation of the vestibular system</article-title>. <source>Exp. Brain Res.</source> <volume>178</volume>, <fpage>151</fpage>&#x02013;<lpage>166</lpage>. <pub-id pub-id-type="doi">10.1007/s00221-006-0721-2</pub-id><pub-id pub-id-type="pmid">17031681</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kiemel</surname> <given-names>T.</given-names></name> <name><surname>Zhang</surname> <given-names>Y.</given-names></name> <name><surname>Jeka</surname> <given-names>J. J.</given-names></name></person-group> (<year>2011</year>). <article-title>Identification of neural feedback for upright stance in humans: stabilization rather than sway minimization</article-title>. <source>J. Neurosci.</source> <volume>31</volume>, <fpage>15144</fpage>&#x02013;<lpage>15153</lpage>. <pub-id pub-id-type="doi">10.1523/JNEUROSCI.1013-11.2011</pub-id><pub-id pub-id-type="pmid">22016548</pub-id></citation></ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname> <given-names>M.</given-names></name> <name><surname>Collins</surname> <given-names>S. H.</given-names></name></person-group> (<year>2015</year>). <article-title>Once-per-step control of ankle-foot prosthesis push-off work reduces effort associated with balance during walking</article-title>. <source>J. Neuroeng. Rehabil.</source> <volume>12</volume>:<fpage>43</fpage>. <pub-id pub-id-type="doi">10.1186/s12984-015-0027-3</pub-id><pub-id pub-id-type="pmid">25928176</pub-id></citation></ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kuo</surname> <given-names>A. D.</given-names></name></person-group> (<year>2007</year>). <article-title>The six determinants of gait and the inverted pendulum analogy: a dynamic walking perspective</article-title>. <source>Hum. Mov. Sci.</source> <volume>26</volume>, <fpage>617</fpage>&#x02013;<lpage>656</lpage>. <pub-id pub-id-type="doi">10.1016/j.humov.2007.04.003</pub-id><pub-id pub-id-type="pmid">17617481</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kuznetsova</surname> <given-names>A.</given-names></name> <name><surname>Brockhoff</surname> <given-names>P. B.</given-names></name> <name><surname>Christensen</surname> <given-names>R. H. B.</given-names></name></person-group> (<year>2017</year>). <article-title>lmerTest package: tests in linear mixed effects models</article-title>. <source>J. Stat. Softw.</source> <volume>82</volume>, <fpage>1</fpage>&#x02013;<lpage>26</lpage>. <pub-id pub-id-type="doi">10.18637/jss.v082.i13</pub-id></citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Logan</surname> <given-names>D.</given-names></name> <name><surname>Ivanenko</surname> <given-names>Y. P.</given-names></name> <name><surname>Kiemel</surname> <given-names>T.</given-names></name> <name><surname>Cappellini</surname> <given-names>G.</given-names></name> <name><surname>Sylos-Labini</surname> <given-names>F.</given-names></name> <name><surname>Lacquaniti</surname> <given-names>F.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Function dictates the phase dependence of vision during human locomotion</article-title>. <source>J. Neurophysiol.</source> <volume>112</volume>, <fpage>165</fpage>&#x02013;<lpage>180</lpage>. <pub-id pub-id-type="doi">10.1152/jn.01062.2012</pub-id><pub-id pub-id-type="pmid">24717345</pub-id></citation></ref>
<ref id="B24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lu</surname> <given-names>T. W.</given-names></name> <name><surname>O&#x00027;Connor</surname> <given-names>J. J.</given-names></name></person-group> (<year>1999</year>). <article-title>Bone position estimation from skin marker co-ordinates using global optimisation with joint constraints</article-title>. <source>J. Biomechan.</source> <volume>32</volume>, <fpage>129</fpage>&#x02013;<lpage>134</lpage>. <pub-id pub-id-type="pmid">10052917</pub-id></citation></ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Maurer</surname> <given-names>C.</given-names></name> <name><surname>Mergner</surname> <given-names>T.</given-names></name> <name><surname>Peterka</surname> <given-names>R. J.</given-names></name></person-group> (<year>2006</year>). <article-title>Multisensory control of human upright stance</article-title>. <source>Expe. Brain Res.</source> <volume>171</volume>, <fpage>231</fpage>&#x02013;<lpage>250</lpage>. <pub-id pub-id-type="doi">10.1007/s00221-005-0256-y</pub-id><pub-id pub-id-type="pmid">16307252</pub-id></citation></ref>
<ref id="B26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Morasso</surname> <given-names>P.</given-names></name> <name><surname>Schieppati</surname> <given-names>M.</given-names></name></person-group> (<year>1999</year>). <article-title>Can muscle stiffness alone stabilize upright standing?</article-title> <source>J. Neurophysiol.</source> <volume>82</volume>, <fpage>1622</fpage>&#x02013;<lpage>1626</lpage>. <pub-id pub-id-type="pmid">10482776</pub-id></citation></ref>
<ref id="B27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Morris</surname> <given-names>M. E.</given-names></name> <name><surname>Iansek</surname> <given-names>R.</given-names></name> <name><surname>Matyas</surname> <given-names>T. A.</given-names></name> <name><surname>Summers</surname> <given-names>J. J.</given-names></name></person-group> (<year>1994</year>). <article-title>Ability to modulate walking cadence remains intact in Parkinson&#x00027;s disease</article-title>. <source>J. Neurol. Neurosurg. Psychiat.</source> <volume>57</volume>, <fpage>1532</fpage>&#x02013;<lpage>1534</lpage>. <pub-id pub-id-type="pmid">7798986</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Muggeo</surname> <given-names>V. M. R.</given-names></name></person-group> (<year>2008</year>). <article-title>Segmented: an R package to fit regression models with broken-line relationships</article-title>. <source>R News</source> <volume>8</volume>, <fpage>20</fpage>&#x02013;<lpage>25</lpage>.</citation></ref>
<ref id="B29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>O&#x00027;Connor</surname> <given-names>S. M.</given-names></name> <name><surname>Kuo</surname> <given-names>A. D.</given-names></name></person-group> (<year>2009</year>). <article-title>Direction-dependent control of balance during walking and standing</article-title>. <source>J. Neurophysiol.</source> <volume>102</volume>, <fpage>1411</fpage>&#x02013;<lpage>1419</lpage>. <pub-id pub-id-type="doi">10.1152/jn.00131.2009</pub-id><pub-id pub-id-type="pmid">19553493</pub-id></citation></ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Peterka</surname> <given-names>R. J.</given-names></name></person-group> (<year>2000</year>). <article-title>Postural control model interpretation of stabilogram diffusion analysis</article-title>. <source>Biol. Cybernet.</source> <volume>82</volume>, <fpage>335</fpage>&#x02013;<lpage>343</lpage>. <pub-id pub-id-type="doi">10.1007/s004220050587</pub-id><pub-id pub-id-type="pmid">10804065</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Peterka</surname> <given-names>R. J.</given-names></name></person-group> (<year>2002</year>). <article-title>Sensorimotor integration in human postural control</article-title>. <source>J. Neurophysiol.</source> <volume>88</volume>, <fpage>1097</fpage>&#x02013;<lpage>118</lpage>. <pub-id pub-id-type="doi">10.1152/jn.2002.88.3.1097</pub-id><pub-id pub-id-type="pmid">12205132</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="book"><person-group person-group-type="author"><collab>R Core Team</collab></person-group> (<year>2013</year>). <source>R: A Language and Environment for Statistical Computing.</source> <publisher-loc>Vienna</publisher-loc>: <publisher-name>R Core Team</publisher-name>. <pub-id pub-id-type="pmid">30628467</pub-id></citation></ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Reimann</surname> <given-names>H.</given-names></name> <name><surname>Fettrow</surname> <given-names>T.</given-names></name> <name><surname>Jeka</surname> <given-names>J. J.</given-names></name></person-group> (<year>2018a</year>). <article-title>Strategies for the control of balance during locomotion</article-title>. <source>Kinesiol. Rev.</source> <volume>7</volume>, <fpage>18</fpage>&#x02013;<lpage>25</lpage>. <pub-id pub-id-type="doi">10.1123/kr.2017-0053</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Reimann</surname> <given-names>H.</given-names></name> <name><surname>Fettrow</surname> <given-names>T.</given-names></name> <name><surname>Thompson</surname> <given-names>E. D.</given-names></name> <name><surname>Agada</surname> <given-names>P.</given-names></name> <name><surname>McFadyen</surname> <given-names>B. J.</given-names></name> <name><surname>Jeka</surname> <given-names>J. J.</given-names></name></person-group> (<year>2017</year>). <article-title>Complementary mechanisms for upright balance during walking</article-title>. <source>PLoS ONE</source> <volume>12</volume>:<fpage>e0172215</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0172215</pub-id><pub-id pub-id-type="pmid">28234936</pub-id></citation></ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Reimann</surname> <given-names>H.</given-names></name> <name><surname>Fettrow</surname> <given-names>T.</given-names></name> <name><surname>Thompson</surname> <given-names>E. D.</given-names></name> <name><surname>Jeka</surname> <given-names>J. J.</given-names></name></person-group> (<year>2018b</year>). <article-title>Neural control of balance during walking</article-title>. <source>Front. Physiol.</source> <volume>9</volume>:<fpage>1271</fpage>. <pub-id pub-id-type="doi">10.3389/fphys.2018.01271</pub-id><pub-id pub-id-type="pmid">30271354</pub-id></citation></ref>
<ref id="B36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sloot</surname> <given-names>L.</given-names></name> <name><surname>van der Krogt</surname> <given-names>M.</given-names></name> <name><surname>Harlaar</surname> <given-names>J.</given-names></name></person-group> (<year>2014</year>). <article-title>Self-paced versus fixed speed treadmill walking</article-title>. <source>Gait Post.</source> <volume>39</volume>, <fpage>478</fpage>&#x02013;<lpage>484</lpage>. <pub-id pub-id-type="doi">10.1016/j.gaitpost.2013.08.022</pub-id><pub-id pub-id-type="pmid">24055003</pub-id></citation></ref>
<ref id="B37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Song</surname> <given-names>S.</given-names></name> <name><surname>Geyer</surname> <given-names>H.</given-names></name></person-group> (<year>2015</year>). <article-title>A neural circuitry that emphasizes spinal feedback generates diverse behaviours of human locomotion</article-title>. <source>J. Physiol.</source> <volume>593</volume>, <fpage>3493</fpage>&#x02013;<lpage>3511</lpage>. <pub-id pub-id-type="doi">10.1113/JP270228</pub-id><pub-id pub-id-type="pmid">25920414</pub-id></citation></ref>
<ref id="B38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Taga</surname> <given-names>G.</given-names></name></person-group> (<year>1998</year>). <article-title>A model of the neuro-musculo-skeletal system for anticipatory adjustment of human locomotion during obstacle avoidance</article-title>. <source>Biol. Cybernet.</source> <volume>78</volume>, <fpage>9</fpage>&#x02013;<lpage>17</lpage>. <pub-id pub-id-type="pmid">9485584</pub-id></citation></ref>
<ref id="B39">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Tylkowski</surname> <given-names>C. M.</given-names></name> <name><surname>Simon</surname> <given-names>S. R.</given-names></name> <name><surname>Mansour</surname> <given-names>J. M.</given-names></name></person-group> (<year>1982</year>). <article-title>Internal rotation gait in spastic cerebral palsy in the hip</article-title>, in <source>Proceedings of the 10th Open Scientific Meeting of the Hip Society</source>, ed <person-group person-group-type="editor"><name><surname>Nelson</surname> <given-names>J. P.</given-names></name></person-group> (St. <publisher-loc>Louis, MO</publisher-loc>: <publisher-name>Mosby</publisher-name>), <fpage>89</fpage>&#x02013;<lpage>125</lpage>.</citation></ref>
<ref id="B40">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>van der Kooij</surname> <given-names>H.</given-names></name> <name><surname>Jacobs</surname> <given-names>R.</given-names></name> <name><surname>Koopman</surname> <given-names>B.</given-names></name> <name><surname>Grootenboer</surname> <given-names>H.</given-names></name></person-group> (<year>1999</year>). <article-title>A multisensory integration model of human stance control</article-title>. <source>Biol. Cybernet.</source> <volume>80</volume>, <fpage>299</fpage>&#x02013;<lpage>308</lpage>. <pub-id pub-id-type="pmid">10365423</pub-id></citation></ref>
<ref id="B41">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Vlutters</surname> <given-names>M.</given-names></name> <name><surname>van Asseldonk</surname> <given-names>E. H. F.</given-names></name> <name><surname>Van Der Kooij</surname> <given-names>H.</given-names></name></person-group> (<year>2016</year>). <article-title>Center of mass velocity-based predictions in balance recovery following pelvis perturbations during human walking</article-title>. <source>J. Exp. Biol.</source> <volume>219</volume>. <fpage>1514</fpage>&#x02013;<lpage>1523</lpage>. <pub-id pub-id-type="doi">10.1242/jeb.129338</pub-id><pub-id pub-id-type="pmid">26994171</pub-id></citation></ref>
<ref id="B42">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Vlutters</surname> <given-names>M.</given-names></name> <name><surname>van Asseldonk</surname> <given-names>E. H. F.</given-names></name> <name><surname>van der Kooij</surname> <given-names>H.</given-names></name></person-group> (<year>2019</year>). <article-title>Ankle muscle responses during perturbed walking with blocked ankle joints</article-title>. <source>J. Neurophysiol.</source> <volume>121</volume>, <fpage>1711</fpage>&#x02013;<lpage>1717</lpage>. <pub-id pub-id-type="doi">10.1152/jn.00752.2018</pub-id><pub-id pub-id-type="pmid">30864874</pub-id></citation></ref>
<ref id="B43">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>Y.</given-names></name> <name><surname>Srinivasan</surname> <given-names>M.</given-names></name></person-group> (<year>2014</year>). <article-title>Stepping in the direction of the fall: the next foot placement can be predicted from current upper body state in steady-state walking</article-title>. <source>Biol. Lett.</source> <volume>10</volume>:<fpage>20140405</fpage>. <pub-id pub-id-type="doi">10.1098/rsbl.2014.0405</pub-id><pub-id pub-id-type="pmid">25252834</pub-id></citation></ref>
<ref id="B44">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Winter</surname> <given-names>D. A.</given-names></name></person-group> (<year>1990</year>). <source>Biomechanics and Motor Control of Human Movement</source>. <publisher-loc>Hoboken, NJ</publisher-loc>: <publisher-name>Wiley</publisher-name>.</citation></ref>
<ref id="B45">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Winter</surname> <given-names>D. A.</given-names></name></person-group> (<year>1995</year>). <article-title>Human balance and posture control during standing and walking</article-title>. <source>Gait Post.</source> <volume>3</volume>, <fpage>193</fpage>&#x02013;<lpage>214</lpage>.</citation></ref>
<ref id="B46">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname> <given-names>J.-F.</given-names></name> <name><surname>Stein</surname> <given-names>R. B.</given-names></name></person-group> (<year>1990</year>). <article-title>Phase-dependent reflex reversal in human leg muscles during walking</article-title>. <source>J. Neurophysiol.</source> <volume>63</volume>, <fpage>1109</fpage>&#x02013;<lpage>1117</lpage>. <pub-id pub-id-type="pmid">2358865</pub-id></citation></ref>
<ref id="B47">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zehr</surname> <given-names>E.</given-names></name> <name><surname>Stein</surname> <given-names>R. B.</given-names></name></person-group> (<year>1999</year>). <article-title>What functions do reflexes serve during human locomotion?</article-title> <source>Prog. Neurobiol.</source> <volume>58</volume>, <fpage>185</fpage>&#x02013;<lpage>205</lpage>. <pub-id pub-id-type="pmid">10338359</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn fn-type="financial-disclosure"><p><bold>Funding.</bold> HR was funded from a Research Fellowship from the German Research Foundation (RE 3780/1-1). HR, TF, and JJ were funded by the National Science Foundation (NSF CRCNS 1822568).</p>
</fn>
</fn-group>
</back>
</article>