Parabolic movement primitives and cortical states: merging optimality with geometric invariance

Previous studies have suggested that several types of rules govern the generation of complex arm movements. One class of rules consists of optimizing an objective function (e.g., maximizing motion smoothness). Another class consists of geometric and kinematic constraints, for instance the coupling between speed and curvature during drawing movements as expressed by the two-thirds power law. It has also been suggested that complex movements are composed of simpler elements or primitives. However, the ability to unify the different rules has remained an open problem. We address this issue by identifying movement paths whose generation according to the two-thirds power law yields maximally smooth trajectories. Using equi-affine differential geometry we derive a mathematical condition which these paths must obey. Among all possible solutions only parabolic paths minimize hand jerk, obey the two-thirds power law and are invariant under equi-affine transformations (which preserve the fit to the two-thirds power law). Affine transformations can be used to generate any parabolic stroke from an arbitrary parabolic template, and a few parabolic strokes may be concatenated to compactly form a complex path. To test the possibility that parabolic elements are used to generate planar movements, we analyze monkeys’ scribbling trajectories. Practiced scribbles are well approximated by long parabolic strokes. Of the motor cortical neurons recorded during scribbling more were related to equi-affine than to Euclidean speed. Unsupervised segmentation of simulta- neously recorded multiple neuron activity yields states related to distinct parabolic elements. We thus suggest that the cortical representation of movements is state-dependent and that parabolic elements are building blocks used by the motor system to generate complex movements.     

A Compact Representation of Drawing Movements with Sequences of Parabolic Primitives

Some studies suggest that complex arm movements in humans and monkeys may optimize several objective functions, while others claim that arm movements satisfy geometric constraints and are composed of elementary components. However, the ability to unify different constraints has remained an open question. The criterion for a maximally smooth (minimizing jerk) motion is satisfied for parabolic trajectories having constant equi-affine speed, which thus comply with the geometric constraint known as the two-thirds power law. Here we empirically test the hypothesis that parabolic segments provide a compact representation of spontaneous drawing movements. Monkey scribblings performed during a period of practice were recorded. Practiced hand paths could be approximated well by relatively long parabolic segments. Following practice, the orientations and spatial locations of the fitted parabolic segments could be drawn from only 2–4 clusters, and there was less discrepancy between the fitted parabolic segments and the executed paths. This enabled us to show that well-practiced spontaneous scribbling movements can be represented as sequences (“words”) of a small number of elementary parabolic primitives (“letters”). A movement primitive can be defined as a movement entity that cannot be intentionally stopped before its completion. We found that in a well-trained monkey a movement was usually decelerated after receiving a reward, but it stopped only after the completion of a sequence composed of several parabolic segments. Piece-wise parabolic segments can be generated by applying affine geometric transformations to a single parabolic template. Thus, complex movements might be constructed by applying sequences of suitable geometric transformations to a few templates. Our findings therefore suggest that the motor system aims at achieving more parsimonious internal representations through practice, that parabolas serve as geometric primitives and that non-Euclidean variables are employed in internal movement representations (due to the special role of parabolas in equi-affine geometry).     

The Relation between Geometry and Time in Mental Actions

Mental imagery is a cognitive tool that helps humans take decisions by simulating past and future events. The hypothesis has been advanced that there is a functional equivalence between actual and mental movements. Yet, we do not know whether there are any limitations to its validity even in terms of some fundamental features of actual movements, such as the relationship between space and time. Although it is impossible to directly measure the spatiotemporal features of mental actions, an indirect investigation can be conducted by taking advantage of the constraints existing in planar drawing movements and described by the two-thirds power law (2/3PL). This kinematic law describes one of the most impressive regularities observed in biological movements: movement speed decreases when curvature increases. Here, we compared the duration of identical actual and mental arm movements by changing the constraints imposed by the 2/3PL. In the first two experiments, the length of the trajectory remained constant, while its curvature (Experiment 1) or its number of inflexions (Experiment 2) was manipulated. The results showed that curvature, but not the number of inflexions, proportionally and similarly affected actual and mental movement duration, as expected from the 2/3PL. Two other control experiments confirmed that the results of Experiment 1 were not attributable to eye movements (Experiment 3) or to the perceived length of the displayed trajectory (Experiment 4). Altogether, our findings suggest that mental movement simulation is tuned to the kinematic laws characterizing actions and that kinematics of actual and mental movements is completely specified by the representation of their geometry.     

Elegant‐crested Tinamous Eudromia elegans do not synchronize head and leg movements during head‐bobbing

Head‐bobbing is the fore–aft movement of the head relative to the body during terrestrial locomotion in birds. It is considered to be a behaviour that helps to stabilize images on the retina during locomotion, yet some studies have suggested biomechanical links between the movements of the head and legs. This study analysed terrestrial locomotion and head‐bobbing in the Elegant‐crested Tinamou Eudromia elegans at a range of speeds by synchronously recording high‐speed video and ground reaction forces in a laboratory setting. The results indicate that the timing of head and leg movements are dissociated from one another. Nonetheless, head and neck movements do affect stance duration, ground reaction forces and body pitch and, as a result, the movement of the centre of mass in head‐bobbing birds. This study does not support the hypothesis that head‐bobbing is itself constrained by terrestrial locomotion. Instead, it suggests that visual cues are the primary trigger for head‐bobbing in birds, and locomotion is, in turn, constrained by a need for image stabilization and depth perception.     

Muscle activity in rapid multi-degree-of-freedom elbow movements: solutions from a musculoskeletal model

The activity of certain muscles that cross the elbow joint complex (EJC) are affected by forearm position and forearm movement during elbow flexion/extension. To investigate whether these changes are based on the musculoskeletal geometry of the joint, a three-dimensional musculotendinoskeletal computer model of the EJC was used to estimate individual muscle activity in multi-degree-of-freedom (df) rapid (ballistic) elbow movements. It is hypothesized that this model could reproduce the major features of elbow muscle activity during multi-df elbow movements using dynamic optimal control theory, given a minimum-time performance criterion. Results from the model are presented and verified with experimental kinematic and electromyographic data from movements that involved both one-df elbow flexion/extension and two-df flexion/extension with forearm pronation/supination. The model demonstrated how the activity of particular muscles is affected by both forearm position and movement, as measured in these experiments and as previously reported by others. These changes were most evident in the flexor muscles and least evident in the extensor muscles. The model also indicated that, for specific one- and two-df movements, activating a muscle that is antagonistic or noncontributory to the movement could reduce the movement time. The major features of muscle activity in multi-df elbow movements appear to be highly dependent on the joint's musculoskeletal geometry and are not strictly based on neural influences or neuroanatomical substrates. Received: 9 May 1997 / Accepted in revised form: 8 December 1998     

Analysis of kinematic invariances of multijoint reaching movement

 There is a no unique relationship between the trajectory of the hand, represented in cartesian or extrinsic space, and its trajectory in joint angle or intrinsic space in the general condition of joint redundancy. The goal of this work is to analyze the relation between planning the trajectory of a multijoint movement in these two coordinate systems. We show that the cartesian trajectory can be planned based on the task parameters (target coordinates, etc.) prior to and independently of angular trajectories. Angular time profiles are calculated from the cartesian trajectory to serve as a basis for muscle control commands. A unified differential equation that allows planning trajectories in cartesian and angular spaces simultaneously is proposed. Due to joint redundancy, each cartesian trajectory corresponds to a family of angular trajectories which can account for the substantial variability of the latter. A set of strategies for multijoint motor control following from this model is considered; one of them coincides with the frog wiping reflex model and resolves the kinematic inverse problem without inversion. The model trajectories exhibit certain properties observed in human multijoint reaching movements such as movement equifinality, straight end-point paths, bell-shaped tangential velocity profiles, speed-sensitive and speed-insensitive movement strategies, peculiarities of the response to double-step targets, and variations of angular trajectory without variations of the limb end-point trajectory in cartesian space. In humans, those properties are almost independent of limb configuration, target location, movement duration, and load. In the model, these properties are invariant to an affine transform of cartesian space. This implies that these properties are not a special goal of the motor control system but emerge from movement kinematics that reflect limb geometry, dynamics, and elementary principles of motor control used in planning. All the results are given analytically and, in order to compare the model with experimental results, by computer simulations. Received: 6 April 1994/Accepted in revised form: 25 April 1995     

Kinematics of the Coordination of Pointing during Locomotion

In natural motor behaviour arm movements, such as pointing or reaching, often need to be coordinated with locomotion. The underlying coordination patterns are largely unexplored, and require the integration of both rhythmic and discrete movement primitives. For the systematic and controlled study of such coordination patterns we have developed a paradigm that combines locomotion on a treadmill with time-controlled pointing to targets in the three-dimensional space, exploiting a virtual reality setup. Participants had to walk at a constant velocity on a treadmill. Synchronized with specific foot events, visual target stimuli were presented that appeared at different spatial locations in front of them. Participants were asked to reach these stimuli within a short time interval after a “go” signal. We analysed the variability patterns of the most relevant joint angles, as well as the time coupling between the time of pointing and different critical timing events in the foot movements. In addition, we applied a new technique for the extraction of movement primitives from kinematic data based on anechoic demixing. We found a modification of the walking pattern as consequence of the arm movement, as well as a modulation of the duration of the reaching movement in dependence of specific foot events. The extraction of kinematic movement primitives from the joint angle trajectories exploiting the new algorithm revealed the existence of two distinct main components accounting, respectively, for the rhythmic and discrete components of the coordinated movement pattern. Summarizing, our study shows a reciprocal pattern of influences between the coordination patterns of reaching and walking. This pattern might be explained by the dynamic interactions between central pattern generators that initiate rhythmic and discrete movements of the lower and upper limbs, and biomechanical factors such as the dynamic gait stability.     

Self-organized control of bipedal locomotion by neural oscillators in unpredictable environment

A new principle of sensorimotor control of legged locomotion in an unpredictable environment is proposed on the basis of neurophysiological knowledge and a theory of nonlinear dynamics. Stable and flexible locomotion is realized as a global limit cycle generated by a global entrainment between the rhythmic activities of a nervous system composed of coupled neural oscillators and the rhythmic movements of a musculo-skeletal system including interaction with its environment. Coordinated movements are generated not by slaving to an explicit representation of the precise trajectories of the movement of each part but by dynamic interactions among the nervous system, the musculo-skeletal system and the environment. The performance of a bipedal model based on the above principle was investigated by computer simulation. Walking movements stable to mechanical perturbations and to environmental changes were obtained. Moreover, the model generated not only the walking movement but also the running movement by changing a single parameter nonspecific to the movement. The transitions between the gait patterns occurred with hysteresis.     

The Temporal Structure of Vertical Arm Movements

The present study investigates how the CNS deals with the omnipresent force of gravity during arm motor planning. Previous studies have reported direction-dependent kinematic differences in the vertical plane; notably, acceleration duration was greater during a downward than an upward arm movement. Although the analysis of acceleration and deceleration phases has permitted to explore the integration of gravity force, further investigation is necessary to conclude whether feedforward or feedback control processes are at the origin of this incorporation. We considered that a more detailed analysis of the temporal features of vertical arm movements could provide additional information about gravity force integration into the motor planning. Eight subjects performed single joint vertical arm movements (45° rotation around the shoulder joint) in two opposite directions (upwards and downwards) and at three different speeds (slow, natural and fast). We calculated different parameters of hand acceleration profiles: movement duration (MD), duration to peak acceleration (D PA), duration from peak acceleration to peak velocity (D PA-PV), duration from peak velocity to peak deceleration (D PV-PD), duration from peak deceleration to the movement end (D PD-End), acceleration duration (AD), deceleration duration (DD), peak acceleration (PA), peak velocity (PV), and peak deceleration (PD). While movement durations and amplitudes were similar for upward and downward movements, the temporal structure of acceleration profiles differed between the two directions. More specifically, subjects performed upward movements faster than downward movements; these direction-dependent asymmetries appeared early in the movement (i.e., before PA) and lasted until the moment of PD. Additionally, PA and PV were greater for upward than downward movements. Movement speed also changed the temporal structure of acceleration profiles. The effect of speed and direction on the form of acceleration profiles is consistent with the premise that the CNS optimises motor commands with respect to both gravitational and inertial constraints.     

What moves us? How mobility and movement are at the center of human evolution

Movement is central to the survival of all free‐living organisms. Consequently, movement and what anthropologists often refer to as mobility, which is the sum of small‐scale movements tracked across larger geographic and temporal scales, are key targets of selection. Movement and mobility also underpin many of the key features that make us human and that allowed our lineage to adapt to changing environments across the globe. The most obvious example is the evolution of humans' singular mode of locomotion. Bipedalism is arguably the most important derived anatomical trait of the hominin lineage. The mechanisms and circumstances that gave rise to this novel mode of movement remain subjects of intense research.     

Rhythmic dynamics and synchronization via dimensionality reduction: application to human gait

Reliable characterization of locomotor dynamics of human walking is vital to understanding the neuromuscular control of human locomotion and disease diagnosis. However, the inherent oscillation and ubiquity of noise in such non-strictly periodic signals pose great challenges to current methodologies. To this end, we exploit the state-of-the-art technology in pattern recognition and, specifically, dimensionality reduction techniques, and propose to reconstruct and characterize the dynamics accurately on the cycle scale of the signal. This is achieved by deriving a low-dimensional representation of the cycles through global optimization, which effectively preserves the topology of the cycles that are embedded in a high-dimensional Euclidian space. Our approach demonstrates a clear advantage in capturing the intrinsic dynamics and probing the subtle synchronization patterns from uni/bivariate oscillatory signals over traditional methods. Application to human gait data for healthy subjects and diabetics reveals a significant difference in the dynamics of ankle movements and ankle-knee coordination, but not in knee movements. These results indicate that the impaired sensory feedback from the feet due to diabetes does not influence the knee movement in general, and that normal human walking is not critically dependent on the feedback from the peripheral nervous system.     

Rotational locomotion by the cockroach Blattella germanica

Locomotion on complex substrata can be expressed in a plane by two geometric components of body movement: linear locomotion and rotational locomotion. This study examined pure rotation by analysing the geometry of leg movements and stepping patterns during the courtship turns of male Blattella germanica. Strict rotation or translation by an insect requires that each side of the body cover equal distance with respect to the substrate. There are three mechanisms by which the legs can maintain this equality: frequency of stepping, magnitude of the leg arcs relative to the body and the degree to which legs flex and extend during locomotion. During the courtship behaviour of Blattella germanica selected males executed turns involving body rotation along with leg movements in which the legs on the outside of the turn swung through greater average arcs than those on the inside of the turn. This difference should have resulted in a translation component. However, legs on the inside of the turn compensated by flexion and extension movements which were greater than those of opposing legs. The net effect was that both sides of the body covered equal average ground. These cockroaches used a wide variety of stepping combinations to effect rotation. The frequency of these combinations was compared to an expected frequency distribution of stepping combinations and further to an expected frequency of these stepping combinations used for straight walking. These comparisons demonstrated a similarity between interleg coordination during straight walking and that during turning in place.     

Are interlimb interactions disturbed in patients with Parkinson’s disease? A study under unloading conditions

During natural human locomotion, neural connections are activated that are typical of regulation of the quadrupedal walking. The interaction between the neural networks generating rhythmic movements of the upper and lower limbs depends on tonic state of each of these networks regulated by motor signals from the brain. Distortion of these signals in patients with Parkinson’s disease (PD) may lead to disruption of the interlimb interactions. We examined the effect of movements of the limbs of one girdle on the parameters of the motor activity of another limb girdle at their joint cyclic movements under the conditions of arm and leg unloading in 17 patients with PD and 16 healthy subjects. We have shown that, in patients, the effect of voluntary and passive movements of arms, as well as the active movement of the distal parts of arms, on the voluntary movement of legs is weak, while in healthy subjects, the effect of arm movements on the parameters of voluntary stepping is significant. The effect of arm movements on the activation of the involuntary stepping by vibrational stimulation of-legs in patients was absent, while in healthy subjects, the motor activity of arms increased the possibility of involuntary rhythmic movements activation. Differences in the effect of leg movements on the rhythmic movements of arms were found in both patients and healthy subjects. The interlimb interaction appeared after drug administration. However, the effect of the drug was not sufficient for the recovery of normal state of the neural networks in patients. In PD patients, neural networks generating stepping rhythm have an increased tonic activity, which prevents the activation and appearance of involuntary rhythmic movements facilitating the effects of arms on legs.     

Behavioral states help translate dispersal movements into spatial distribution patterns of floaters

Within the field of spatial ecology, it is important to study animal movements in order to better understand population dynamics. Dispersal is a nonlinear process through which different behavioral mechanisms could affect movement patterns. One of the most common approaches to analyzing the trajectories of organisms within patches is to use random-walk models to describe movement features. These models express individual movements within a specific area in terms of random-walk parameters in an effort to relate movement patterns to the distributions of organisms in space. However, only using the movement trajectories of individuals to predict the spatial spread of animal populations may not fit the complex distribution of individuals across heterogeneous environments. When we empirically tested the results from a random-walk model (a residence index) used to predict the spatial equilibrium distribution of individuals, we found that the index severely underestimated the spatial spread of dispersing individuals. We believe this is because random-walk models only account for the effects of environmental conditions on individual movements, completely overlooking the crucial influence of behavior changes over time. In the future, both aspects should be accounted for when predicting general rules of (meta)population abundance, distribution, and dynamics from patterns of animal movements.     

Kinematic invariants during cyclical arm movements

It has been observed that the motion of the arm end-point (the hand, fingertip or the tip of a pen) is characterized by a number of regularities (kinematic invariants). Trajectory is usually straight, and the velocity profile has a bell shape during point-to-point movements. During drawing movements, a two-thirds power law predicts the dependence of the end-point velocity on the trajectory curvature. Although various principles of movement organization have been discussed as possible origins of these kinematic invariants, the nature of these movement trajectory characteristics remains an open question. A kinematic model of cyclical arm movements derived in the present study analytically demonstrates that all three kinematic invariants can be predicted from a two-joint approximation of the kinematic structure of the arm and from sinusoidal joint motions. With this approach, explicit expressions for two kinematic invariants, the two-thirds power law during drawing movements and the velocity profile during point-to-point movements are obtained as functions of arm segment lengths and joint motion parameters. Additionally, less recognized kinematic invariants are also derived from the model. The obtained analytical expressions are further validated with experimental data. The high accuracy of the predictions confirms practical utility of the model, showing that the model is relevant to human performance over a wide range of movements. The results create a basis for the consolidation of various existing interpretations of kinematic invariants. In particular, optimal control is discussed as a plausible source of invariant characteristics of joint motions and movement trajectories.     

Kinematic parameters of terrestrial locomotion in cursorial (ratites), swimming (ducks), and striding birds (quail and guinea fowl).

The importance of size, functional features and morphological features in adaptation for walking in birds were studied. The time and space kinematic parameters of locomotion were compared in two running birds, the ratites (rhea, kiwi, Paleognatiforms), in two swimming birds, (ducks) and two striding birds, (quail and Guinea fowl). The results showed that in the two phases, stance and swing, the time and space parameters worked in opposite ways: the duration of the swing was constant, but its length increased with speed. In contrast, the duration of the stance was correlated to speed, while its length was not (except in ducks). In all the birds, a higher speed was achieved by a decrease of the stance duration, and an increase of the swing length. The kinematic parameters were not used in the same way in all species: There is a size effect and large birds increase their speed mainly by increasing the frequency of their movements and the small species increase mainly their amplitude. Nevertheless, it is not the main factor and morphology, such as swimming adaptation features of the ducks, and behaviour, are important because they modify the mechanical constraints and influence the kinematics parameters.     

Body segment differences in surface area, skin temperature and 3D displacement and the estimation of heat balance during locomotion in hominins

The conventional method of estimating heat balance during locomotion in humans and other hominins treats the body as an undifferentiated mass. This is problematic because the segments of the body differ with respect to several variables that can affect thermoregulation. Here, we report a study that investigated the impact on heat balance during locomotion of inter-segment differences in three of these variables: surface area, skin temperature and rate of movement. The approach adopted in the study was to generate heat balance estimates with the conventional method and then compare them with heat balance estimates generated with a method that takes into account inter-segment differences in surface area, skin temperature and rate of movement. We reasoned that, if the hypothesis that inter-segment differences in surface area, skin temperature and rate of movement affect heat balance during locomotion is correct, the estimates yielded by the two methods should be statistically significantly different. Anthropometric data were collected on seven adult male volunteers. The volunteers then walked on a treadmill at 1.2 m/s while 3D motion capture cameras recorded their movements. Next, the conventional and segmented methods were used to estimate the volunteers' heat balance while walking in four ambient temperatures. Lastly, the estimates produced with the two methods were compared with the paired t-test. The estimates of heat balance during locomotion yielded by the two methods are significantly different. Those yielded by the segmented method are significantly lower than those produced by the conventional method. Accordingly, the study supports the hypothesis that inter-segment differences in surface area, skin temperature and rate of movement impact heat balance during locomotion. This has important implications not only for current understanding of heat balance during locomotion in hominins but also for how future research on this topic should be approached.     

Adaptive design of locomotion and foot form in prosobranch gastropods

A comparison of the locomotor types, speed, tenacity, and foot form of nearly 300 species in 52 families of marine prosobranchs has revealed that foot size and shape and even subtle variations of locomotion affect the speed and strength of adhesion to the substratum.Gastropods inhabiting soft substrata move primarily by pedal cilia or by discontinuous locomotion in which shell and foot move alternately. Both types of movement are accompanied by low tenacity. A specialized type of discontinuous locomotion, namely, leaping, surpasses all other methods of movement in speed. Species with ciliary locomotion have a very large foot while those with discontinuous movement have an exceedingly small foot relative to shell size.The majority of prosobranchs inhabit hard substrata, move by continuous pedal muscular gliding, and have moderately high tenacity during movement. Arhythmic pedal locomotion yields lower maximum speeds and tenacities than do rhythmic pedal waves. Foot size and shape relative to shell length in species with arhythmic locomotion vary from very short and broad to long and narrow. Studies of transects at several temperate and tropical marine littoral stations showed that these species are confined to low littoral or sublittoral habitats that are sheltered from heavy wave action. High speed and tenacity are simultaneously attained only by species with rhythmic pedal waves.Speed and tenacity do not represent competing selective pressures on the size and shape of the foot. Speed increases among species as the foot approaches or exceeds shell length and is highest if the foot is also broad; the greatest tenacities are attained by species with a long, broad foot whose dimensions do not exceed that of the shell. The optimal shape for both high tenacity and speed is a broad foot somewhat shorter than the shell; neither speed nor tenacity are much compromised by this form. In general, only species with rhythmic pedal waves whose foot size and shape approximate the optimal form for high tenacity and speed are found in habitats exposed to much wave action. Long rhythmic waves, moving a large proportion of foot area at once, are in theory energetically more economical than small, very rapid waves resulting in the same overall speed, but experiments showed that tenacity is significantly reduced in gastropods which increase speed by enlarging the waves. The optimal wave pattern of a species should be a balance between the demand for speed with the least expenditure of energy, favored by a pattern of many large waves at once, and the demand for tenacity, favored by a pattern of few and small waves.Retrograde ditaxic waves of elongation are the most common pattern encountered among prosobranchs, and are associated with a large range of foot sizes and shapes. Such waves are at least one third as long as the foot, while direct waves and other waves of compression are frequently much smaller. The range of foot forms of species with waves of compression is restricted, tending to be optimal for high tenacity or to be long and narrow. Waves of compression appear to be a specialization with the potential for maintaining high tenacity even at high speeds since the waves can be very small, and for giving superior speed since they can travel very rapidly.     

Motor primitives in vertebrates and invertebrates

In recent years different lines of evidence have led to the idea that motor actions and movements in both vertebrates and invertebrates are composed of elementary building blocks. The entire motor repertoire can be spanned by applying a well-defined set of operations and transformations to these primitives and by combining them in many different ways according to well-defined syntactic rules. Motor and movement primitives and modules might exist at the neural, dynamic and kinematic levels with complicated mapping among the elementary building blocks subserving these different levels of representation. Hence, while considerable progress has been made in recent years in unravelling the nature of these primitives, new experimental, computational and conceptual approaches are needed to further advance our understanding of motor compositionality.