Estimating the phase of synchronized oscillators

Interlimb coordination is not strictly controlled during walking

In human walking, the left and right legs move alternately, half a stride out of phase with each other. Although various parameters, such as stride frequency and length, vary with walking speed, the antiphase relationship remains unchanged. In contrast, during walking in left-right asymmetric situations, the relative phase shifts from the antiphase condition to compensate for the asymmetry. Interlimb coordination is important for adaptive walking and we expect that interlimb coordination is strictly controlled during walking. However, the control mechanism remains unclear. In the present study, we derived a quantity that models the control of interlimb coordination during walking using two coupled oscillators based on the phase reduction theory and Bayesian inference method. The results were not what we expected. Specifically, we found that the relative phase is not actively controlled until the deviation from the antiphase condition exceeds a certain threshold. In other words, the control of interlimb coordination has a dead zone like that in the case of the steering wheel of an automobile. It is conjectured that such forgoing of control enhances energy efficiency and maneuverability. Our discovery of the dead zone in the control of interlimb coordination provides useful insight for understanding gait control in humans.

Gait signature changes with walking speed are similar among able-bodied young adults despite persistent individual-specific differences

Understanding individuals' distinct movement patterns is crucial for health, rehabilitation, and sports. Recently, we developed a machine learning-based framework to show that "gait signatures" describing the neuromechanical dynamics governing able-bodied and post-stroke gait kinematics remain individual-specific across speeds. However, we only evaluated gait signatures within a limited speed range and number of participants, using only sagittal plane (i.e., 2D) joint angles. Here we characterized changes in gait signatures across a wide range of speeds, from very slow (0.3 m/s) to exceptionally fast (above the walk-to-run transition speed) in 17 able-bodied young adults. We further assessed whether 3D kinematic and/or kinetic (ground reaction forces, joint moments, and powers) data would improve the discrimination of gait signatures. Our study showed that gait signatures remained individual-specific across walking speeds: Notably, 3D kinematic signatures achieved exceptional accuracy (99.8%, confidence interval (CI) 99.1-100%) in classifying individuals, surpassing both 2D kinematics and 3D kinetics. Moreover, participants exhibited consistent, predictable linear changes in their gait signatures across the entire speed range. These changes were associated with participants' preferred walking speeds, balance ability, cadence, and step length. These findings support gait signatures as a tool to characterize individual differences in gait and predict speed-induced changes in gait dynamics.

Discovering individual-specific gait signatures from data-driven models of neuromechanical dynamics

Locomotion results from the interactions of highly nonlinear neural and biomechanical dynamics. Accordingly, understanding gait dynamics across behavioral conditions and individuals based on detailed modeling of the underlying neuromechanical system has proven difficult. Here, we develop a data-driven and generative modeling approach that recapitulates the dynamical features of gait behaviors to enable more holistic and interpretable characterizations and comparisons of gait dynamics. Specifically, gait dynamics of multiple individuals are predicted by a dynamical model that defines a common, low-dimensional, latent space to compare group and individual differences. We find that highly individualized dynamics-i.e., gait signatures-for healthy older adults and stroke survivors during treadmill walking are conserved across gait speed. Gait signatures further reveal individual differences in gait dynamics, even in individuals with similar functional deficits. Moreover, components of gait signatures can be biomechanically interpreted and manipulated to reveal their relationships to observed spatiotemporal joint coordination patterns. Lastly, the gait dynamics model can predict the time evolution of joint coordination based on an initial static posture. Our gait signatures framework thus provides a generalizable, holistic method for characterizing and predicting cyclic, dynamical motor behavior that may generalize across species, pathologies, and gait perturbations.

Walking is like slithering: A unifying, data-driven view of locomotion

Legged movement is ubiquitous in nature and of increasing interest for robotics. Most legged animals routinely encounter foot slipping, yet detailed modeling of multiple contacts with slipping exceeds current simulation capacity. Here we present a principle that unifies multilegged walking (including that involving slipping) with slithering and Stokesian (low Reynolds number) swimming. We generated data-driven principally kinematic models of locomotion for walking in low-slip animals (Argentine ant, 4.7% slip ratio of slipping to total motion) and for high-slip robotic systems (BigANT hexapod, slip ratio 12 to 22%; Multipod robots ranging from 6 to 12 legs, slip ratio 40 to 100%). We found that principally kinematic models could explain much of the variability in body velocity and turning rate using body shape and could predict walking behaviors outside the training data. Most remarkably, walking was principally kinematic irrespective of leg number, foot slipping, and turning rate. We find that grounded walking, with or without slipping, is governed by principally kinematic equations of motion, functionally similar to frictional swimming and slithering. Geometric mechanics thus leads to a unified model for swimming, slithering, and walking. Such commonality may shed light on the evolutionary origins of animal locomotion control and offer new approaches for robotic locomotion and motion planning.

Analysis and control of a running spring-mass model with a trunk based on virtual pendulum concept

The spring-loaded inverted pendulum model has been one of the most studied conceptual models in the locomotion community. Even though it can adequately explain the center of mass trajectories of numerous legged animals, it remains insufficient in template-based control of complex robot platforms, being unable to capture additional dynamic characteristics of locomotion exhibited in additional degrees of freedom such as trunk pitch oscillations. In fact, analysis of trunk behavior during locomotion has been one of the motivations behind studying the virtual pivot point (VPP) concept, with biological inspiration and basis for both natural and synthetic systems with non-negligible trunk dynamics. This study first presents a comprehensive analysis of the VPP concept for planar running behaviors, followed by a systematic study of the existence and characteristics of periodic solutions. In particular, we investigate how periodic solutions depend on model control parameters and compare them based on stability and energetic cost. We then develop a feedback controller that can stabilize system dynamics around its periodic solutions and evaluate performance as compared to a previously introduced controller from the literature. We demonstrate the effectiveness of both controllers and find that the proposed control scheme creates larger basins of attraction with minor degradation in convergence speed. In conclusion, this study shows that the VPP concept, in conjunction with the proposed controller, could be beneficial in designing and controlling legged robots capable of running with non-trivial upper body dynamics. Our systematic analysis of periodic solutions arising from the use of the VPP concept is also an important step towards a more formal basis for comparisons of the VPP concept with bio-locomotion.

Gait Entrainment to Torque Pulses from a Hip Exoskeleton Robot

Robot-aided locomotor rehabilitation has proven challenging. To facilitate progress, it is important to first understand the neuro-mechanical dynamics and control of unimpaired human locomotion. Our previous studies found that human gait entrained to periodic torque pulses at the ankle when the pulse period was close to preferred stride duration. Moreover, synchronized gait exhibited constant phase relation with the pulses so that the robot provided mechanical assistance. To test the generality of mechanical gait entrainment, this study characterized unimpaired human subjects' responses to periodic torque pulses during overground walking. The intervention was applied by a hip exoskeleton robot, Samsung GEMS-H. Gait entrainment was assessed based on the time-course of the phase at which torque pulses occurred within each stride. Experiments were conducted for two consecutive days to evaluate whether the second day elicited more entrainment. Whether entrainment was affected by the difference between pulse period and preferred stride duration was also assessed. Results indicated that the intervention evoked gait entrainment that occurred more often when the period of perturbation was closer to subjects' preferred stride duration, but the difference between consecutive days was insignificant. Entrainment was accompanied by convergence of pulse phase to a similar value across all conditions, where the robot maximized mechanical assistance. Clear evidence of motor adaptation indicated the potential of the intervention for rehabilitation. This study quantified important aspects of the nonlinear neuro-mechanical dynamics underlying unimpaired human walking, which will inform the development of effective approaches to robot-aided locomotor rehabilitation, exploiting natural dynamics in a minimally-encumbering way.

Measuring the repertoire of age-related behavioral changes in Drosophila melanogaster

Aging affects almost all aspects of an organism—its morphology, its physiology, its behavior. Isolating which biological mechanisms are regulating these changes, however, has proven difficult, potentially due to our inability to characterize the full repertoire of an animal’s behavior across the lifespan. Using data from fruit flies (D. melanogaster) we measure the full repertoire of behaviors as a function of age. We observe a sexually dimorphic pattern of changes in the behavioral repertoire during aging. Although the stereotypy of the behaviors and the complexity of the repertoire overall remains relatively unchanged, we find evidence that the observed alterations in behavior can be explained by changing the fly’s overall energy budget, suggesting potential connections between metabolism, aging, and behavior.

Aging is a ubiquitous biological phenomenon that affects many aspects of an animal’s appearance, physiology, and behavior. Our understanding of how changes in physiology lead to behavioral changes, however, has been partially limited by our ability to robustly quantify how behavior alters over timescales of days and weeks. In this study, we measure a large repertoire of behaviors of fruit flies at various ages, finding how the actions the animals perform shift with age. We observe a difference between the aging dynamics of male and female flies, and we show that many of these changes can be explained with a model of energy consumption, leading us to make predictions as to the role of metabolism in changes in aging behavior.

Data-driven geometric system identification for shape-underactuated dissipative systems

Modeling system dynamics becomes challenging when the properties of individual system components cannot be directly measured, and often requires identification of properties from observed motion. In this paper, we show that systems whose movement is highly dissipative have features which provide an opportunity to more easily identify models and more quickly optimize motions than would be possible with general techniques. Geometric mechanics provides means for reduction of the dynamics by environmental homogeneity, while the dissipative nature minimizes the role of second order (inertial) features in the dynamics. Here we extend the tools of geometric system identification to 'shape-underactuated dissipative systems (SUDS)'-systems whose motions are more dissipative than inertial, but whose actuation is restricted to a subset of the body shape coordinates. Many animal motions are SUDS, including micro-swimmers such as nematodes and flagellated bacteria, and granular locomotors such as snakes and lizards. Many soft robots are also SUDS, particularly robots that incorporate highly damped series elastic actuators to reduce the rigidity of their interactions with their environments during locomotion and manipulation. We motivate the use of SUDS models, and validate their ability to predict motion of a variety of simulated viscous swimming platforms. For a large class of SUDS, we show how the shape velocity actuation inputs can be directly converted into torque inputs, suggesting that systems with soft pneumatic or dielectric elastomer actuators can be modeled with the tools presented. Based on fundamental assumptions in the physics, we show how our model complexity scales linearly with the number of passive shape coordinates. This scaling offers a large reduction on the number of trials needed to identify the system model from experimental data, and may reduce overfitting. The sample efficiency of our method suggests its use in modeling, control, and optimization in robotics, and as a tool for the study of organismal motion in friction dominated regimes.

Rapid frequency modulation in a resonant system: aerial perturbation recovery in hawkmoths

Centimetre-scale fliers must contend with the high power requirements of flapping flight. Insects have elastic elements in their thoraxes which may reduce the inertial costs of their flapping wings. Matching wingbeat frequency to a mechanical resonance can be energetically favourable, but also poses control challenges. Many insects use frequency modulation on long timescales, but wingstroke-to-wingstroke modulation of wingbeat frequencies in a resonant spring-wing system is potentially costly because muscles must work against the elastic flight system. Nonetheless, rapid frequency and amplitude modulation may be a useful control modality. The hawkmoth Manduca sexta has an elastic thorax capable of storing and returning significant energy. However, its nervous system also has the potential to modulate the driving frequency of flapping because its flight muscles are synchronous. We tested whether hovering hawkmoths rapidly alter frequency during perturbations with vortex rings. We observed both frequency modulation (32% around mean) and amplitude modulation (37%) occurring over several wingstrokes. Instantaneous phase analysis of wing kinematics revealed that more than 85% of perturbation responses required active changes in neurogenic driving frequency. Unlike their robotic counterparts that abdicate frequency modulation for energy efficiency, synchronous insects use wingstroke-to-wingstroke frequency modulation despite the power demands required for deviating from resonance.

Dog galloping on rough terrain exhibits similar limb co-ordination patterns and gait variability to that on flat terrain

Understanding how animals regulate their gait during locomotion can give biological insight and inspire controllers for robots. Why animals use the gallop at the highest speeds remains incompletely explained. Hypothesized reasons for galloping include that it enables recruitment of spinal musculoskeletal structures, that it minimizes energy losses as predicted by collisional theory, or that it provides extended flight phases with more time for leg placement and hence enhances or provides necessary maneuverability [Alexander 1988 Am. Zool. 28 237-45; Ruina, Bertram and Srinivasan 2005 J. Theor. Biol. 237 170-92; Usherwood 2019 J. Exp. Zool. Part A 333 9-19; Hildebrand1989 Bioscience 39 766-75]. The latter-most hypothesis has implications in robotics, where controllers based on the concept of multistability have gained some traction. Here we examine this hypothesis by studying the dynamics of dog gait on flat and rough terrain. This hypothesis predicts that injection of noise into timing and location of ground contacts during the galloping gait by rough terrain will result in an isotropically more noisy gallop gait, centered around the gallop used on flat terrain. We find that dog gait in terms of leg swing timing on rough terrain is not consistently more variable about the mean gait, and constrain the upper limits of this variability to values that are unlikely to be biologically relevant. However the location of the mean gait indeed only shifts by a small amount. Therefore, we find limited support for this hypothesis. This suggests that achieving a target gallop gait with tight regulation is still the desired behavior, and that large amounts of variability in gait are not a desired feature of the gallop. For robotics, our results suggest that the emergent animal-environment dynamics on rough terrain do not exhibit uniformly wider basins of attraction. Future robotics work could test whether controllers that do or do not allow shifts in mean gait and gait variability produce more economical and/or stable gallops.

Drosophila uses a tripod gait across all walking speeds, and the geometry of the tripod is important for speed control

Changes in walking speed are characterized by changes in both the animal's gait and the mechanics of its interaction with the ground. Here we study these changes in walking Drosophila. We measured the fly's center of mass movement with high spatial resolution and the position of its footprints. Flies predominantly employ a modified tripod gait that only changes marginally with speed. The mechanics of a tripod gait can be approximated with a simple model - angular and radial spring-loaded inverted pendulum (ARSLIP) - which is characterized by two springs of an effective leg that become stiffer as the speed increases. Surprisingly, the change in the stiffness of the spring is mediated by the change in tripod shape rather than a change in stiffness of individual legs. The effect of tripod shape on mechanics can also explain the large variation in kinematics among insects, and ARSLIP can model these variations.

Predicting walking response to ankle exoskeletons using data-driven models

Despite recent innovations in exoskeleton design and control, predicting subject-specific impacts of exoskeletons on gait remains challenging. We evaluated the ability of three classes of subject-specific phase-varying (PV) models to predict kinematic and myoelectric responses to ankle exoskeletons during walking, without requiring prior knowledge of specific user characteristics. Each model-PV, linear PV (LPV) and nonlinear PV (NPV)-leveraged Floquet theory to predict deviations from a nominal gait cycle due to exoskeleton torque, though the models differed in complexity and expected prediction accuracy. For 12 unimpaired adults walking with bilateral passive ankle exoskeletons, we predicted kinematics and muscle activity in response to three exoskeleton torque conditions. The LPV model's predictions were more accurate than the PV model when predicting less than 12.5% of a stride in the future and explained 49-70% of the variance in hip, knee and ankle kinematic responses to torque. The LPV model also predicted kinematic responses with similar accuracy to the more-complex NPV model. Myoelectric responses were challenging to predict with all models, explaining at most 10% of the variance in responses. This work highlights the potential of data-driven PV models to predict complex subject-specific responses to ankle exoskeletons and inform device design and control.

Gait Phase Recognition Using Deep Convolutional Neural Network with Inertial Measurement Units

Gait phase recognition is of great importance in the development of assistance-as-needed robotic devices, such as exoskeletons. In order for a powered exoskeleton with phase-based control to determine and provide proper assistance to the wearer during gait, the user’s current gait phase must first be identified accurately. Gait phase recognition can potentially be achieved through input from wearable sensors. Deep convolutional neural networks (DCNN) is a machine learning approach that is widely used in image recognition. User kinematics, measured from inertial measurement unit (IMU) output, can be considered as an ‘image’ since it exhibits some local ‘spatial’ pattern when the sensor data is arranged in sequence. We propose a specialized DCNN to distinguish five phases in a gait cycle, based on IMU data and classified with foot switch information. The DCNN showed approximately 97% accuracy during an offline evaluation of gait phase recognition. Accuracy was highest in the swing phase and lowest in terminal stance.

Information-based centralization of locomotion in animals and robots

The centralization of locomotor control from weak and local coupling to strong and global is hard to assess outside of particular modeling frameworks. We developed an empirical, model-free measure of centralization that compares information between control signals and both global and local states. A second measure, co-information, quantifies the net redundancy in global and local control. We first validate that our measures predict centralization in simulations of phase-coupled oscillators. We then test how centralization changes with speed in freely running cockroaches. Surprisingly, across all speeds centralization is constant and muscle activity is more informative of the global kinematic state (the averages of all legs) than the local state of that muscle's leg. Finally we use a legged robot to show that mechanical coupling alone can change the centralization of legged locomotion. The results of these systems span a design space of centralization and co-information for biological and robotic systems.

The manifold structure of limb coordination in walking Drosophila

Terrestrial locomotion requires animals to coordinate their limb movements to efficiently traverse their environment. While previous studies in hexapods have reported that limb coordination patterns can vary substantially, the structure of this variability is not yet well understood. Here, we characterized the symmetric and asymmetric components of variation in walking kinematics in the genetic model organism Drosophila. We found that Drosophila use a single continuum of coordination patterns without evidence for preferred configurations. Spontaneous symmetric variability was associated with modulation of a single control parameter-stance duration-while asymmetric variability consisted of small, limb-specific modulations along multiple dimensions of the underlying symmetric pattern. Commands that modulated walking speed, originating from artificial neural activation or from the visual system, evoked modulations consistent with spontaneous behavior. Our findings suggest that Drosophila employ a low-dimensional control architecture, which provides a framework for understanding the neural circuits that regulate hexapod legged locomotion.

Locomotion: exploiting noise for state estimation

Running, walking, flying and swimming are all processes in which animals produce propulsion by executing rhythmic motions of their bodies. Dynamical stability of the locomotion is hardly automatic: millions of older people are injured by falling each year. Stability frequently requires sensory feedback. We investigate how organisms obtain the information they use in maintaining their stability. Assessing stability of a periodic orbit of a dynamical system requires information about the dynamics of the system off the orbit. For locomotion driven by a periodic orbit, perturbations that "kick" the trajectory off the orbit must occur in order to observe convergence rates toward the orbit. We propose that organisms generate excitations in order to set the gains for stabilizing feedback. We hypothesize further that these excitations are stochastic but have heavy-tailed, non-Gaussian probability distributions. Compared to Gaussian distributions, we argue that these are more effective for estimating stability characteristics of the orbit. Finally, we propose experiments to test the efficacy of these ideas.

Phase reduction and phase-based optimal control for biological systems: a tutorial

A powerful technique for the analysis of nonlinear oscillators is the rigorous reduction to phase models, with a single variable describing the phase of the oscillation with respect to some reference state. An analog to phase reduction has recently been proposed for systems with a stable fixed point, and phase reduction for periodic orbits has recently been extended to take into account transverse directions and higher-order terms. This tutorial gives a unified treatment of such phase reduction techniques and illustrates their use through mathematical and biological examples. It also covers the use of phase reduction for designing control algorithms which optimally change properties of the system, such as the phase of the oscillation. The control techniques are illustrated for example neural and cardiac systems.

A model of temporal scaling correctly predicts that motor timing improves with speed

Timing is fundamental to complex motor behaviors: from tying a knot to playing the piano. A general feature of motor timing is temporal scaling: the ability to produce motor patterns at different speeds. One theory of temporal processing proposes that the brain encodes time in dynamic patterns of neural activity (population clocks), here we first examine whether recurrent neural network (RNN) models can account for temporal scaling. Appropriately trained RNNs exhibit temporal scaling over a range similar to that of humans and capture a signature of motor timing, Weber's law, but predict that temporal precision improves at faster speeds. Human psychophysics experiments confirm this prediction: the variability of responses in absolute time are lower at faster speeds. These results establish that RNNs can account for temporal scaling and suggest a novel psychophysical principle: the Weber-Speed effect.

Limping following limb loss increases locomotor stability

Although many arthropods have the ability to voluntarily lose limbs, how these animals rapidly adapt to such an extreme perturbation remains poorly understood. It is thought that moving with certain gaits can enable efficient, stable locomotion; however, switching gaits requires complex information flow between and coordination of an animal's limbs. We show here that upon losing two legs, spiders can switch to a novel, more statically stable gait, or use temporal adjustments without a gait change. The resulting gaits have higher overall static stability than the gaits that would be imposed by limb loss. By decreasing the time spent in a low-stability configuration - effectively 'limping' over less-stable phases of the stride - spiders increased the overall stability of the less statically stable gait with no observable reduction in speed, as compared with the intact condition. Our results shed light on how voluntary limb loss could have persisted evolutionarily among many animals, and provide bioinspired solutions for robots when they break or lose limbs.

Uncovering the structure of the mouse gait controller: Mice respond to substrate perturbations with adaptations in gait on a continuum between trot and bound

Animals, including humans, have been shown to maintain a gait during locomotion that minimizes the risk of injury and energetic cost. Despite the importance of understanding the mechanisms of gait regulation, ethical and experimental challenges have prevented full exploration of these. Here we present data on the gait response of mice to rapid, precisely timed, spatially confined mechanical perturbations. Our data elucidate that after the mechanical perturbation, the mouse gait response is anisotropic, preferring deviations away from the trot towards bounding, over those towards other gaits, such as walk or pace. We quantified this shift by projecting the observed gait onto the line between trot and bound, in the space of quadrupedal gaits. We call this projection λ. For λ=0, the gait is the ideal trot; for λ=±π, it is the ideal bound. We found that the substrate perturbation caused a significant shift in λ towards bound during the stride in which the perturbation occurred and the following stride (linear mixed effects model: Δλ=0.26±0.07 and Δλ=0.21±0.07, respectively; random effect for animal, p < 0.05 for both strides, n = 8 mice). We hypothesize that this is because the bounding gait is better suited to rapid acceleration or deceleration, and an exploratory analysis of jerk showed that it was significantly correlated with λ (p < 0.05). Understanding how gait is controlled under perturbations can aid in diagnosing gait pathologies and in the design of more agile robots.

Improving the quality of a collective signal in a consumer EEG headset

This work focuses on the experimental data analysis of electroencephalography (EEG) data, in which multiple sensors are recording oscillatory voltage time series. The EEG data analyzed in this manuscript has been acquired using a low-cost commercial headset, the Emotiv EPOC+. Our goal is to compare different techniques for the optimal estimation of collective rhythms from EEG data. To this end, a traditional method such as the principal component analysis (PCA) is compared to more recent approaches to extract a collective rhythm from phase-synchronized data. Here, we extend the work by Schwabedal and Kantz (PRL 116, 104101 (2016)) evaluating the performance of the Kosambi-Hilbert torsion (KHT) method to extract a collective rhythm from multivariate oscillatory time series and compare it to results obtained from PCA. The KHT method takes advantage of the singular value decomposition algorithm and accounts for possible phase lags among different time series and allows to focus the analysis on a specific spectral band, optimally amplifying the signal-to-noise ratio of a common rhythm. We evaluate the performance of these methods for two particular sets of data: EEG data recorded with closed eyes and EEG data recorded while observing a screen flickering at 15 Hz. We found an improvement in the signal-to-noise ratio of the collective signal for the KHT over the PCA, particularly when random temporal shifts are added to the channels.

Behaviorally Selective Engagement of Short-Latency Effector Pathways by Motor Cortex

Blocking motor cortical output with lesions or pharmacological inactivation has identified movements that require motor cortex. Yet, when and how motor cortex influences muscle activity during movement execution remains unresolved. We addressed this ambiguity using measurement and perturbation of motor cortical activity together with electromyography in mice during two forelimb movements that differ in their requirement for cortical involvement. Rapid optogenetic silencing and electrical stimulation indicated that short-latency pathways linking motor cortex with spinal motor neurons are selectively activated during one behavior. Analysis of motor cortical activity revealed a dramatic change between behaviors in the coordination of firing patterns across neurons that could account for this differential influence. Thus, our results suggest that changes in motor cortical output patterns enable a behaviorally selective engagement of short-latency effector pathways. The model of motor cortical influence implied by our findings helps reconcile previous observations on the function of motor cortex.

Longitudinal quasi-static stability predicts changes in dog gait on rough terrain

Legged animals utilize gait selection to move effectively and must recover from environmental perturbations. We show that on rough terrain, domestic dogs, Canis lupus familiaris, spend more time in longitudinal quasi-statically stable patterns of movement. Here, longitudinal refers to the rostro-caudal axis. We used an existing model in the literature to quantify the longitudinal quasi-static stability of gaits neighbouring the walk, and found that trot-like gaits are more stable. We thus hypothesized that when perturbed, the rate of return to a stable gait would depend on the direction of perturbation, such that perturbations towards less quasi-statically stable patterns of movement would be more rapid than those towards more stable patterns of movement. The net result of this would be greater time spent in longitudinally quasi-statically stable patterns of movement. Limb movement patterns in which diagonal limbs were more synchronized (those more like a trot) have higher longitudinal quasi-static stability. We therefore predicted that as dogs explored possible limb configurations on rough terrain at walking speeds, the walk would shift towards trot. We gathered experimental data quantifying dog gait when perturbed by rough terrain and confirmed this prediction using GPS and inertial sensors (n=6, P<0.05). By formulating gaits as trajectories on the n-torus we are able to make tractable the analysis of gait similarity. These methods can be applied in a comparative study of gait control which will inform the ultimate role of the constraints and costs impacting locomotion, and have applications in diagnostic procedures for gait abnormalities, and in the development of agile legged robots.

Summary: Dogs co-ordinate their limbs on rough terrain in a manner consistent with optimization for quasi-static longitudinal stability.

A Robust Parameterization of Human Gait Patterns Across Phase-Shifting Perturbations

The phase of human gait is difficult to quantify accurately in the presence of disturbances. In contrast, recent bipedal robots use time-independent controllers relying on a mechanical phase variable to synchronize joint patterns through the gait cycle. This concept has inspired studies to determine if human joint patterns can also be parameterized by a mechanical variable. Although many phase variable candidates have been proposed, it remains unclear which, if any, provide a robust representation of phase for human gait analysis or control. In this paper we analytically derive an ideal phase variable (the hip phase angle) that is provably monotonic and bounded throughout the gait cycle. To examine the robustness of this phase variable, ten able-bodied human subjects walked over a platform that randomly applied phase-shifting perturbations to the stance leg. A statistical analysis found the correlations between nominal and perturbed joint trajectories to be significantly greater when parameterized by the hip phase angle (0.95+) than by time or a different phase variable. The hip phase angle also best parameterized the transient errors about the nominal periodic orbit. Finally, interlimb phasing was best explained by local (ipsilateral) hip phase angles that are synchronized during the double-support period.

A Perturbation Mechanism for Investigations of Phase Variables in Human Locomotion

The concept of a phase variable, a mechanical measurement of the body's progression through the gait cycle, has been used to parameterize the leg joint patterns of autonomous bipedal robots, producing human-like gaits with robustness to external perturbations. It was recently proposed that the kinematic response of humans to a perturbation could also be parameterized by a phase variable. In order to properly study this phase variable hypothesis with human subjects, a custom perturbation mechanism was built to cause phase shifts in the gait cycle. The main goals of this study are to introduce the design of a novel perturbation mechanism and experimentally demonstrate its ability to effect phase changes during the gait cycle.

Optimal Extraction of Collective Oscillations from Unreliable Measurements

We present a method that facilitates a phase description of collective, irregular-oscillatory dynamics from unreliable multichannel recordings. The collective oscillations may be represented in each channel with fluctuating amplitude, phase offsets, and substantial amounts of measurement noise. Our method performs well under such realistic conditions, as we exemplify with collective brain rhythms in multichannel electroencephalogram recordings.

Mapping the stereotyped behaviour of freely moving fruit flies

A frequent assumption in behavioural science is that most of an animal's activities can be described in terms of a small set of stereotyped motifs. Here, we introduce a method for mapping an animal's actions, relying only upon the underlying structure of postural movement data to organize and classify behaviours. Applying this method to the ground-based behaviour of the fruit fly, Drosophila melanogaster, we find that flies perform stereotyped actions roughly 50% of the time, discovering over 100 distinguishable, stereotyped behavioural states. These include multiple modes of locomotion and grooming. We use the resulting measurements as the basis for identifying subtle sex-specific behavioural differences and revealing the low-dimensional nature of animal motions.

Planar Covariation of Hindlimb and Forelimb Elevation Angles during Terrestrial and Aquatic Locomotion of Dogs

The rich repertoire of locomotor behaviors in quadrupedal animals requires flexible inter-limb and inter-segmental coordination. Here we studied the kinematic coordination of different gaits (walk, trot, gallop, and swim) of six dogs (Canis lupus familiaris) and, in particular, the planar covariation of limb segment elevation angles. The results showed significant variations in the relative duration of rearward limb movement, amplitude of angular motion, and inter-limb coordination, with gait patterns ranging from a lateral sequence of footfalls during walking to a diagonal sequence in swimming. Despite these differences, the planar law of inter-segmental coordination was maintained across different gaits in both forelimbs and hindlimbs. Notably, phase relationships and orientation of the covariation plane were highly limb specific, consistent with the functional differences in their neural control. Factor analysis of published muscle activity data also demonstrated differences in the characteristic timing of basic activation patterns of the forelimbs and hindlimbs. Overall, the results demonstrate that the planar covariation of inter-segmental coordination has emerged for both fore- and hindlimbs and all gaits, although in a limb-specific manner.

Dual dimensionality reduction reveals independent encoding of motor features in a muscle synergy for insect flight control

What are the features of movement encoded by changing motor commands? Do motor commands encode movement independently or can they be represented in a reduced set of signals (i.e. synergies)? Motor encoding poses a computational and practical challenge because many muscles typically drive movement, and simultaneous electrophysiology recordings of all motor commands are typically not available. Moreover, during a single locomotor period (a stride or wingstroke) the variation in movement may have high dimensionality, even if only a few discrete signals activate the muscles. Here, we apply the method of partial least squares (PLS) to extract the encoded features of movement based on the cross-covariance of motor signals and movement. PLS simultaneously decomposes both datasets and identifies only the variation in movement that relates to the specific muscles of interest. We use this approach to explore how the main downstroke flight muscles of an insect, the hawkmoth Manduca sexta, encode torque during yaw turns. We simultaneously record muscle activity and turning torque in tethered flying moths experiencing wide-field visual stimuli. We ask whether this pair of muscles acts as a muscle synergy (a single linear combination of activity) consistent with their hypothesized function of producing a left-right power differential. Alternatively, each muscle might individually encode variation in movement. We show that PLS feature analysis produces an efficient reduction of dimensionality in torque variation within a wingstroke. At first, the two muscles appear to behave as a synergy when we consider only their wingstroke-averaged torque. However, when we consider the PLS features, the muscles reveal independent encoding of torque. Using these features we can predictably reconstruct the variation in torque corresponding to changes in muscle activation. PLS-based feature analysis provides a general two-sided dimensionality reduction that reveals encoding in high dimensional sensory or motor transformations.

Understanding movement control is challenging because the brains of nearly all animals send motor command signals to many muscles, and these signals produce complex movements. In studying animal movement, one cannot always record all the motor commands an animal uses or know all the ways in which movement varies in response. A combined approach is necessary to find the relevant patterns: the changes in movement that correspond to changes in the recorded motor commands. Techniques exist to identify simple patterns in either the motor commands or the movements, but in this paper we develop an approach that identifies patterns in both simultaneously. We use this technique to understand how agile flying insects control aerial turns. The two main downstroke muscles of moths are thought to produce turns by creating a power difference between the left and right wings. The moth’s brain may only need to specify the difference in activation between the two muscles. We discover that moth’s brain actually has independent control over each muscle, and this separate control increases the moth’s ability to adjust turning within a single wingstroke. Our computational approach reveals sophisticated patterns of movement processing even in the small nervous systems of insects.

Constructing predictive models of human running

Running is an essential mode of human locomotion, during which ballistic aerial phases alternate with phases when a single foot contacts the ground. The spring-loaded inverted pendulum (SLIP) provides a starting point for modelling running, and generates ground reaction forces that resemble those of the centre of mass (CoM) of a human runner. Here, we show that while SLIP reproduces within-step kinematics of the CoM in three dimensions, it fails to reproduce stability and predict future motions. We construct SLIP control models using data-driven Floquet analysis, and show how these models may be used to obtain predictive models of human running with six additional states comprising the position and velocity of the swing-leg ankle. Our methods are general, and may be applied to any rhythmic physical system. We provide an approach for identifying an event-driven linear controller that approximates an observed stabilization strategy, and for producing a reduced-state model which closely recovers the observed dynamics.

Intersegmental coupling and recovery from perturbations in freely running cockroaches

Cockroaches are remarkably stable runners, exhibiting rapid recovery from external perturbations. To uncover the mechanisms behind this important behavioral trait, we recorded leg kinematics of freely running animals in both undisturbed and perturbed trials. Functional coupling underlying inter-leg coordination was monitored before and during localized perturbations, which were applied to single legs via magnetic impulses. The resulting transient effects on all legs and the recovery times to normal pre-perturbation kinematics were studied. We estimated coupling architecture and strength by fitting experimental data to a six-leg-unit phase oscillator model. Using maximum-likelihood techniques, we found that a network with nearest-neighbor inter-leg coupling best fitted the data and that, although coupling strengths vary among preparations, the overall inputs entering each leg are approximately balanced and consistent. Simulations of models with different coupling strengths encountering perturbations suggest that the coupling schemes estimated from our experiments allow animals relatively fast and uniform recoveries from perturbations.

A survey of phase variable candidates of human locomotion

Studies show that the human nervous system is able to parameterize gait cycle phase using sensory feedback. In the field of bipedal robots, the concept of a phase variable has been successfully used to mimic this behavior by parameterizing the gait cycle in a time-independent manner. This approach has been applied to control a powered transfemoral prosthetic leg, but the proposed phase variable was limited to the stance period of the prosthesis only. In order to achieve a more robust controller, we attempt to find a new phase variable that fully parameterizes the gait cycle of a prosthetic leg. The angle with respect to a global reference frame at the hip is able to monotonically parameterize both the stance and swing periods of the gait cycle. This survey looks at multiple phase variable candidates involving the hip angle with respect to a global reference frame across multiple tasks including level-ground walking, running, and stair negotiation. In particular, we propose a novel phase variable candidate that monotonically parameterizes the whole gait cycle across all tasks, and does so particularly well across level-ground walking. In addition to furthering the design of robust robotic prosthetic leg controllers, this survey could help neuroscientists and physicians study human locomotion across tasks from a time-independent perspective.

Lag, lock, sync, slip: the many 'phases' of coupled flagella

In a multitude of life's processes, cilia and flagella are found indispensable. Recently, the biflagellated chlorophyte alga Chlamydomonas has become a model organism for the study of ciliary motility and synchronization. Here, we use high-speed, high-resolution imaging of single pipette-held cells to quantify the rich dynamics exhibited by their flagella. Underlying this variability in behaviour are biological dissimilarities between the two flagella-termed cis and trans, with respect to a unique eyespot. With emphasis on the wild-type, we derive limit cycles and phase parametrizations for self-sustained flagellar oscillations from digitally tracked flagellar waveforms. Characterizing interflagellar phase synchrony via a simple model of coupled oscillators with noise, we find that during the canonical swimming breaststroke the cis flagellum is consistently phase-lagged relative to, while remaining robustly phase-locked with, the trans flagellum. Transient loss of synchrony, or phase slippage, may be triggered stochastically, in which the trans flagellum transitions to a second mode of beating with attenuated beat envelope and increased frequency. Further, exploiting this alga's ability for flagellar regeneration, we mechanically induced removal of one or the other flagellum of the same cell to reveal a striking disparity between the beatings of the cis and trans flagella, in isolation. These results are evaluated in the context of the dynamic coordination of Chlamydomonas flagella.

Evidence for a time-invariant phase variable in human ankle control

Human locomotion is a rhythmic task in which patterns of muscle activity are modulated by state-dependent feedback to accommodate perturbations. Two popular theories have been proposed for the underlying embodiment of phase in the human pattern generator: a time-dependent internal representation or a time-invariant feedback representation (i.e., reflex mechanisms). In either case the neuromuscular system must update or represent the phase of locomotor patterns based on the system state, which can include measurements of hundreds of variables. However, a much simpler representation of phase has emerged in recent designs for legged robots, which control joint patterns as functions of a single monotonic mechanical variable, termed a phase variable. We propose that human joint patterns may similarly depend on a physical phase variable, specifically the heel-to-toe movement of the Center of Pressure under the foot. We found that when the ankle is unexpectedly rotated to a position it would have encountered later in the step, the Center of Pressure also shifts forward to the corresponding later position, and the remaining portion of the gait pattern ensues. This phase shift suggests that the progression of the stance ankle is controlled by a biomechanical phase variable, motivating future investigations of phase variables in human locomotor control.

Haptic feedback enhances rhythmic motor control by reducing variability, not improving convergence rate

Stability and performance during rhythmic motor behaviors such as locomotion are critical for survival across taxa: falling down would bode well for neither cheetah nor gazelle. Little is known about how haptic feedback, particularly during discrete events such as the heel-strike event during walking, enhances rhythmic behavior. To determine the effect of haptic cues on rhythmic motor performance, we investigated a virtual paddle juggling behavior, analogous to bouncing a table tennis ball on a paddle. Here, we show that a force impulse to the hand at the moment of ball-paddle collision categorically improves performance over visual feedback alone, not by regulating the rate of convergence to steady state (e.g., via higher gain feedback or modifying the steady-state hand motion), but rather by reducing cycle-to-cycle variability. This suggests that the timing and state cues afforded by haptic feedback decrease the nervous system's uncertainty of the state of the ball to enable more accurate control but that the feedback gain itself is unaltered. This decrease in variability leads to a substantial increase in the mean first passage time, a measure of the long-term metastability of a stochastic dynamical system. Rhythmic tasks such as locomotion and juggling involve intermittent contact with the environment (i.e., hybrid transitions), and the timing of such transitions is generally easy to sense via haptic feedback. This timing information may improve metastability, equating to less frequent falls or other failures depending on the task.

Instantaneous kinematic phase reflects neuromechanical response to lateral perturbations of running cockroaches

Instantaneous kinematic phase calculation allows the development of reduced-order oscillator models useful in generating hypotheses of neuromechanical control. When perturbed, changes in instantaneous kinematic phase and frequency of rhythmic movements can provide details of movement and evidence for neural feedback to a system-level neural oscillator with a time resolution not possible with traditional approaches. We elicited an escape response in cockroaches (Blaberus discoidalis) that ran onto a movable cart accelerated laterally with respect to the animals' motion causing a perturbation. The specific impulse imposed on animals (0.50 [Formula: see text] 0.04 m s[Formula: see text]; mean, SD) was nearly twice their forward speed (0.25 [Formula: see text] 0.06 m s[Formula: see text]. Instantaneous residual phase computed from kinematic phase remained constant for 110 ms after the onset of perturbation, but then decreased representing a decrease in stride frequency. Results from direct muscle action potential recordings supported kinematic phase results in showing that recovery begins with self-stabilizing mechanical feedback followed by neural feedback to an abstracted neural oscillator or central pattern generator. Trials fell into two classes of forward velocity changes, while exhibiting statistically indistinguishable frequency changes. Animals pulled away from the side with front and hind legs of the tripod in stance recovered heading within 300 ms, whereas animals that only had a middle leg of the tripod resisting the pull did not recover within this period. Animals with eight or more legs might be more robust to lateral perturbations than hexapods.

Using computational and mechanical models to study animal locomotion

Recent advances in computational methods have made realistic large-scale simulations of animal locomotion possible. This has resulted in numerous mathematical and computational studies of animal movement through fluids and over substrates with the purpose of better understanding organisms' performance and improving the design of vehicles moving through air and water and on land. This work has also motivated the development of improved numerical methods and modeling techniques for animal locomotion that is characterized by the interactions of fluids, substrates, and structures. Despite the large body of recent work in this area, the application of mathematical and numerical methods to improve our understanding of organisms in the context of their environment and physiology has remained relatively unexplored. Nature has evolved a wide variety of fascinating mechanisms of locomotion that exploit the properties of complex materials and fluids, but only recently are the mathematical, computational, and robotic tools available to rigorously compare the relative advantages and disadvantages of different methods of locomotion in variable environments. Similarly, advances in computational physiology have only recently allowed investigators to explore how changes at the molecular, cellular, and tissue levels might lead to changes in performance at the organismal level. In this article, we highlight recent examples of how computational, mathematical, and experimental tools can be combined to ultimately answer the questions posed in one of the grand challenges in organismal biology: "Integrating living and physical systems."

Finding the dimension of slow dynamics in a rhythmic system

Dynamical systems with asymptotically stable periodic orbits are generic models for rhythmic processes in dissipative physical systems. This paper presents a method for reconstructing the dynamics near a periodic orbit from multivariate time-series data. It is used to test theories about the control of legged locomotion, a context in which time series are short when compared with previous work in nonlinear time-series analysis. The method presented here identifies appropriate dimensions of reduced order models for the deterministic portion of the dynamics. The paper also addresses challenges inherent in identifying dynamical models with data from different individuals.

A bioinspired dynamical vertical climbing robot

This paper describes the inspiration for, design, analysis, and implementation of, and experimentation with the first dynamical vertical climbing robot. Biologists have proposed a pendulous climbing model that abstracts remarkable similarities in dynamical wall scaling behavior exhibited by radically different animal species. We numerically study a version of that pendulous climbing template dynamically scaled for applicability to utilitarian payloads with conventional electronics and actuation. This simulation study reveals that the incorporation of passive compliance can compensate for the scaled model’s poorer power density and scale disadvantages relative to biology. However, the introduction of additional dynamical elements raises new concerns about stability regarding both the power stroke and limb coordination schemes that we allay via mathematical analysis of further simplified models. Combining these numerical and analytical insights into a series of design prototypes, we document the correspondence of the various models to the scaled platforms and report that our final prototype climbs dynamically at vertical speeds up to 0.67 m/s (1.5 body-lengths per second, in rough agreement with our models’ predictions).

A single muscle's multifunctional control potential of body dynamics for postural control and running

A neuromechanical approach to control requires understanding how mechanics alters the potential of neural feedback to control body dynamics. Here, we rewrite activation of individual motor units of a behaving animal to mimic the effects of neural feedback without concomitant changes in other muscles. We target a putative control muscle in the cockroach, Blaberus discoidalis (L.), and simultaneously capture limb and body dynamics through high-speed videography and a micro-accelerometer backpack. We test four neuromechanical control hypotheses. We supported the hypothesis that mechanics linearly translates neural feedback into accelerations and rotations during static postural control. However, during running, the same neural feedback produced a nonlinear acceleration control potential restricted to the vertical plane. Using this, we reject the hypothesis from previous work that this muscle acts primarily to absorb energy from the body. The conversion of the control potential is paralleled by nonlinear changes in limb kinematics, supporting the hypothesis that significant mechanical feedback filters the graded neural feedback for running control. Finally, we insert the same neural feedback signal but at different phases in the dynamics. In this context, mechanical feedback enables turning by changing the timing and direction of the accelerations produced by the graded neural feedback.

Stimulus predictability mediates a switch in locomotor smooth pursuit performance for Eigenmannia virescens

The weakly electric glass knifefish, Eigenmannia virescens, will swim forward and backward, using propulsion from an anal ribbon fin, in response to motion of a computer-controlled moving refuge. Fish were recorded performing a refuge-tracking behavior for sinusoidal (predictable) and sum-of-sines (pseudo-random) refuge trajectories. For all trials, we observed high coherence between refuge and fish trajectories, suggesting linearity of the tracking dynamics. But superposition failed: we observed categorical differences in tracking between the predictable single-sine stimuli and the unpredictable sum-of-sines stimuli. This nonlinearity suggests a stimulus-mediated adaptation. At all frequencies tested, fish demonstrated reduced tracking error when tracking single-sine trajectories and this was typically accompanied by a reduction in overall movement. Most notably, fish demonstrated reduced phase lag when tracking single-sine trajectories. These data support the hypothesis that fish generate an internal dynamical model of the stimulus motion, hence improving tracking of predictable trajectories (relative to unpredictable ones) despite similar or reduced motor cost. Similar predictive mechanisms based on the dynamics of stimulus movement have been proposed recently, but almost exclusively for nonlocomotor tasks by humans, such as oculomotor target tracking and posture control. These data suggest that such mechanisms might be common across taxa and behaviors.

Insects running on elastic surfaces

In nature, cockroaches run rapidly over complex terrain such as leaf litter. These substrates are rarely rigid, and are frequently very compliant. Whether and how compliant surfaces change the dynamics of rapid insect locomotion has not been investigated to date largely due to experimental limitations. We tested the hypothesis that a running insect can maintain average forward speed over an extremely soft elastic surface (10 N m(-1)) equal to 2/3 of its virtual leg stiffness (15 N m(-1)). Cockroaches Blaberus discoidalis were able to maintain forward speed (mean +/- s.e.m., 37.2+/-0.6 cm s(-1) rigid surface versus 38.0+/-0.7 cm s(-1) elastic surface; repeated-measures ANOVA, P=0.45). Step frequency was unchanged (24.5+/-0.6 steps s(-1) rigid surface versus 24.7+/-0.4 steps s(-1) elastic surface; P=0.54). To uncover the mechanism, we measured the animal's centre of mass (COM) dynamics using a novel accelerometer backpack, attached very near the COM. Vertical acceleration of the COM on the elastic surface had a smaller peak-to-peak amplitude (11.50+/-0.33 m s(-2), rigid versus 7.7+/-0.14 m s(-2), elastic; P=0.04). The observed change in COM acceleration over an elastic surface required no change in effective stiffness when duty factor and ground stiffness were taken into account. Lowering of the COM towards the elastic surface caused the swing legs to land earlier, increasing the period of double support. A feedforward control model was consistent with the experimental results and provided one plausible, simple explanation of the mechanism.

Fruit flies modulate passive wing pitching to generate in-flight turns

Flying insects execute aerial maneuvers through subtle manipulations of their wing motions. Here, we measure the free-flight kinematics of fruit flies and determine how they modulate their wing pitching to induce sharp turns. By analyzing the torques these insects exert to pitch their wings, we infer that the wing hinge acts as a torsional spring that passively resists the wing's tendency to flip in response to aerodynamic and inertial forces. To turn, the insects asymmetrically change the spring rest angles to generate asymmetric rowing motions of their wings. Thus, insects can generate these maneuvers using only a slight active actuation that biases their wing motion.