JN Watch the video to learn how APS reaches out to developing nations.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Neurophysiol 90: 3774-3782, 2003. First published August 27, 2003; doi:10.1152/jn.00730.2002
0022-3077/03 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
90/6/3774    most recent
00730.2002v2
00730.2002v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (38)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Masani, K.
Right arrow Articles by Nozaki, D.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Masani, K.
Right arrow Articles by Nozaki, D.

Importance of Body Sway Velocity Information in Controlling Ankle Extensor Activities During Quiet Stance

Kei Masani1,2, Milos R. Popovic2, Kimitaka Nakazawa3, Motoki Kouzaki1 and Daichi Nozaki3

1Department of Life Sciences, Graduate School of Arts and Sciences, The University of Tokyo, Tokyo 153-8902, Japan; 2Rehabilitation Engineering Laboratory, Institute of Biomaterials and Biomedical Engineering, University of Toronto, Toronto, Ontario, M5S3G9, Canada; and 3Department of Motor Dysfunction, Research Institute of National Rehabilitation Center for the Disabled, Saitama 359-8555, Japan

Submitted 25 August 2002; accepted in final form 25 August 2003


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 DYNAMICS OF QUIET STANCE
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 ACKNOWLEDGMENTS
 REFERENCES
 
In literature, it has been suggested that the CNS anticipates spontaneous change in body position during quiet stance and continuously modulates ankle extensor muscle activity to compensate for the change. The purpose of this study was to investigate whether velocity feedback contributes by modulating ankle extensor activities in an anticipatory fashion, facilitating effective control of quiet stance. Both theoretical analysis and experiments were carried out to investigate to what extent velocity feedback contributes to controlling quiet stance. The experiments were carried out with 16 healthy subjects who were asked to stand quietly with their eyes open or closed. During the experiments, the center of pressure (COP) displacement (COPdis), the center of mass (COM) displacement (COMdis), and COM velocity (COMvel) in the anteroposterior direction were measured. Rectified electromyograms (EMGs) were used to measure muscle activity in the right soleus muscle, the medial gastrocnemius muscle, and the lateral gastrocnemius muscle. The simulations were performed using an inverted pendulum model that described the anteroposterior kinematics and dynamics of quiet stance. In the simulations, an assumption was made that the COMdis of the body would be regulated using a proportional-derivative (PD) controller. Two different PD controllers were evaluated in these simulations: 1) a controller with the high-derivative/velocity gain (HDG) and 2) a controller with the low-derivative/velocity gain (LDG). Cross-correlation analysis was applied to investigate the relationships between time series obtained in experiments 1) COMdis and EMGs and 2) COMvel and EMGs. Identical cross-correlation analysis was applied to investigate the relationships between time series obtained in simulations 3) COMdis and ankle torque and 4) COMvel and ankle torque. The results of these analyses showed that the COMdis was positively correlated with all three EMGs and that the EMGs temporally preceded the COMdis. These findings agree with the previously published studies in which it was shown that the lateral gastrocnemius muscle is actively modulated in anticipation of the body's COM position change. The COMvel and all three EMGs were also correlated and the cross-correlation function (CCF) had two peaks: one that was positive and another that was negative. The positive peaks were statistically significant, unlike the negative ones; they were larger than the negative peaks; and their time shifts were much shorter compared with the time shifts of the negative peaks. When these results were compared with the CCF results obtained for simulated time series, it was discovered that the cross-correlation results for the HDG controller closely matched cross-correlation results for the experimental time series. On the other hand, the simulation result obtained for LDG controller did not match the experimental results. These findings suggest that the actual postural control system during quiet stance adopts a control strategy that relies notably on velocity information and that such a controller can modulate muscle activity in anticipatory manner without using a feed-forward mechanism.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 DYNAMICS OF QUIET STANCE
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 ACKNOWLEDGMENTS
 REFERENCES
 
Human bipedal stance is inherently unstable because a large body mass is kept in erect posture with its center of mass (COM) located high above a relatively small base of support. The mechanism responsible for equilibrium control of quiet stance has attracted the attention of many researchers in the field. However, the nature of the control mechanism is still an object of controversy.

One approach that has proved useful for investigating human equilibrium control is the analysis of the responses of a quietly standing body to various perturbations, such as perturbations to the proprioceptive subsystem (Fitzpatrick et al. 1992aGo,bGo, 1996Go; Gurfinkel et al. 1995Go; Horak and Nashner 1986Go; Nashner 1976Go, 1977Go; Woollacott et al. 1988Go), perturbations to the visual subsystem (Dijkstra et al. 1994aGo,bGo; Peterka and Benolken 1995Go; Schöner 1991Go), and perturbations to the vestibular subsystem (Day et al. 1997Go; Fitzpatrick et al. 1996Go; Johansson et al. 1995Go; Pavlik et al. 1999Go). These studies have played a significant role in elucidating the contribution made by different subsystems to overall equilibrium control (Dietz 1992Go; Horak and Macpherson 1996Go).

Another approach to investigating human equilibrium control is to analyze spontaneous body sway: trajectory of the center of pressure (COP) (Collins and DeLuca 1993Go, 1994Go, 1995Go; Day et al. 1993Go; Gatev et al. 1999Go; Lacour et al. 1997Go; Mauritz and Dietz 1980Go; Panzer et al. 1995Go), trajectory of the COM (Gatev et al. 1999Go; Panzer et al. 1995Go; Winter et al. 1998Go, 2001Go), trajectory of the ankle joint angle (Fitzpatrick et al. 1994aGo,bGo), and trajectory of other body points (Accornero et al. 1997Go; Aramaki et al. 2001Go). These studies investigated the quantitative and qualitative properties of the spontaneous body sway by comparing the system's behavior under different physiological conditions, such as aging (Accornero et al. 1997Go; Panzer et al. 1995Go), disruption or alteration of proprioception feedback (Fitzpatrick et al. 1994bGo; Mauritz and Dietz 1980Go), disruption or obstruction of vision feedback (Accornero et al. 1997Go; Collins and DeLuca 1995Go; Fitzpatrick et al. 1994bGo), loss or alteration of vestibular sense (Lacour et al. 1997Go), and varying stance width (Day et al. 1993Go; Winter et al. 1998Go).

A group of investigators who adopted the latter approach, Gatev et al. (1999Go) recently reported significant correlation between spontaneous body sway [COM displacement (COMdis) as a representative parameter of body sway] and the activity of the lateral gastrocnemius muscle (LG). They also discovered that the LG activity temporally preceded COMdis. This result could not be interpreted as a linear response of the LG to the positional change of the body. Instead, Gatev et al. (1999Go) proposed that the CNS applies feed-forward control, which anticipates the body position change and activates the LG in advance to regulate balance during quiet stance, rather than a feedback control. Other experimental (Fitzpatrick et al. 1996Go) and theoretical (Morasso and Schieppati 1999Go) studies have also proposed feed-forward control as compensation for the inevitable transmission delay in the neural process, which is adequate to destabilize quiet stance. However, one difficulty with the feed-forward hypothesis is: how does the CNS predict body position?

Over the last 10 years, a number of independent studies have reported that certain postural reactions were adapted according to velocity information (Dijkstra et al. 1994bGo; Jeka et al. 1998Go; Schöner 1991Go). The importance of velocity information in controlling balance during quiet stance was also suggested in simulations carried out by Morasso and Schieppatti in 1999. Velocity feedback can play a significant role in anticipating body position change because it carries information about the subsequent state of the body, i.e., a change in COM velocity (COMvel) indicates the direction and intensity with which the current COMdis will be changed in the following time instant. In general, the velocity feedback in addition to the position feedback, called the proportional + derivative (PD) control, can potentially predict the future condition of the system and can stabilize it more effectively than only a position/proportional controller. Several simulation studies, which applied a single joint inverted pendulum model to simulate human quiet stance, revealed that the PD controller can facilitate stable control of the proposed model (Morasso and Schieppatti 1999Go). However, there is no experimental study without perturbations that investigates the contribution of velocity information in controlling the body during quiet stance, and the structure of the PD controller, i.e., the ratio of position and velocity information, remains unclear.

The purpose of this study was to investigate whether a velocity feedback mechanism makes a significant contribution in an anticipatory modulation of ankle extensor activities during quiet stance. Therefore theoretical and experimental studies were carried out to test this hypothesis.


    DYNAMICS OF QUIET STANCE
 TOP
 ABSTRACT
 INTRODUCTION
 DYNAMICS OF QUIET STANCE
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 ACKNOWLEDGMENTS
 REFERENCES
 
The ankle joint torque needed to stabilize the body during quiet stance can be generated actively and passively.

Passive torque components are the result of tension/stiffness produced by muscle tonus and by the stiffness of the surrounding tissue, such as ligaments and tendons. However, the stabilization of quiet stance by passive torque alone is a very challenging task, and an active component is required to maintain stability (Morasso and Schieppati 1999Go). The active torque component is produced by the CNS, which modulates/controls muscle contractions based on the overall body kinematics and dynamics of spontaneous body sway that are influenced by external disturbances (Morasso and Schieppati 1999Go; Winter 1990Go; Winter et al. 1996Go). The fact that LG activity was modulated in accordance with COMdis in the study by Gatev et al. (1999Go) suggests that this active mechanism is operating.

Herein, dynamics of quiet stance is discussed, as well as generation of the active torque component. For the purpose of this analysis, the human body was approximated as a single segment, single joint inverted pendulum (Fig. 1) that rotates about the ankle joint (Morasso and Schieppati 1999Go; Peterka and Benolken 1995Go; Winter et al. 1998Go). The dynamic equation of the inverted pendulum model is

(1)
where {theta} is the sway angle, theta; is the sway acceleration, m is the mass of the body, I is the moment of inertia of the body (excluding feet), h is the distance of COM from the ankle, g is the gravitational acceleration, T is the total ankle torque, and {epsilon} is the torque disturbance, which is sufficiently small compared with other torque contributions. Notably, ankle torque dominates body movement in this equation. Because COM is located in front of the ankle joint, backward ankle torque is continuously applied to the body to prevent it from falling forward (Smith 1957Go). Because ankle flexor activities are rare and ankle extensors are considerably activated (Gatev et al. 1999Go; Joseph and Nightingale 1952Go; Loram and Lakie 2002Go; Panzer et al. 1995Go), it can be said that ankle extensors contribute the most toward control of the ankle joint torque and therefore the body posture during quiet stance.



View larger version (14K):
[in this window]
[in a new window]
 
FIG. 1. An inverted pendulum model of quiet stance, where y is the center of mass (COM) position, u is the center of pressure (COP) position, f is the ground reaction force divided with fh and fv components, {theta} is the sway angle, g is the acceleration of gravity, T is the total muscle torque about the ankle, and h is the distance of COM from ankle.

 

The ankle torque should satisfy the following equation concerning to the foot segment

(2)
where fv is the vertical component of ground reaction force and u is COP position (COPdis). If we take into account that fv {approx} mg in quiet stance, this equation shows that changes in ankle torque are immediately and linearly translated into changes of COP position. From Eqs. 1 and 2 one can derive

(3)
where y is the COM position. The COP position and COM position must coincide (u = y) under the static equilibrium condition (the inertial term = 0). However, due to excess ankle torque, the COP frequently departs from this instant equilibrium point (Zatsiorsky and Duarte 2000Go) (u != y) and the excess ankle torque generates COM acceleration (theta; != 0). Thus reduction of excess ankle torque to facilitate alignment of the COP position and the COM projection on the standing surface represents the active mechanism for controlling quiet stance.

The high correlation between COPdis and muscle activity (Gatev et al. 1999Go; Schieppati et al. 1994Go) can be accounted for by Eq. 2, which shows a direct relationship between COPdis and the torque generated by muscles. On the other hand, the high correlation between COMdis and LG (Gatev et al. 1999Go) suggests that the ankle torque generated by ankle extensors are almost equal to the gravity torque of the body (mgh sin{theta} {approx} T, in Eq. 1), and that the active mechanism stabilizes the body well.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 DYNAMICS OF QUIET STANCE
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 ACKNOWLEDGMENTS
 REFERENCES
 
To investigate the impact of velocity information on balance control during quiet stance, the following approach was adopted. At first, we obtained kinematic parameters, such as COPdis, COMdis, and COMvel, in the anterioposterior direction and individual EMGs of the ankle extensors during quiet stance. Cross-correlation analysis was applied to investigate the relationship between kinematic parameters and rectified EMGs. We then conducted a theoretical analysis of quiet stance with an inverted pendulum model. In this study, an assumption was made that the COMdis of the pendulum is regulated using a PD controller. Two PD controllers were evaluated in these simulations: 1) a controller with the high-derivative/velocity gain (HDG) and 2) a controller with very low-derivative/velocity gain (LDG). We also obtained the relationship between the kinematic parameters of the inverted pendulum and ankle torque using the cross-correlation analysis. Finally, cross-correlation results of both studies were compared with examination of the hypothesis that the velocity information contributed significantly to the control of quiet standing.

Experiments

Sixteen healthy men (mean ± SD age, 23.8 ± 3.9 yr; mean ± SD height, 169 ± 6.6 cm) participated in this study. All subjects gave informed consent, which was approved by the ethical committee of our research institute.

Each subject was requested to keep a quiet stance posture standing barefoot on a force platform (type 9281B, Kistler, Zürich, Switzerland) with their eyes open (EO condition) or closed (EC condition). The subjects had their arms hanging along the sides of their body, their feet were parallel, and the distance between their heels was 15 cm. The duration of each trial was approximately 40 s, and data from the latter 30 s were subjected to the subsequent analyses. Five trials were conducted for each condition, and sufficient resting time was allowed between trials. In this paper, we have only focused on the anteroposterior body sway, because the ankle extensor muscles are the main contributors to stabilizing body sway in this direction.

The COPdis was obtained using a force platform measurement. The horizontal position of the waist point was measured as an approximation of the COMdis using a CCD (charge coupled device) laser displacement sensor (1 µm resolution; LK-2500, Keyence, Osaka, Japan). This parameter was used to represent body sway. The COMvel was calculated by numerically differentiating the COMdis as a function of time. Electromyograms (EMGs) were recorded by Ag/AgCl surface electrodes with a diameter of 5 mm, which were connected to a preamplifier and a differential amplifier having a bandwidth of 5 Hz to 1 kHz (1253A, NEC Medical Systems, Tokyo, Japan). After careful abrasion of the skin, the electrodes were placed longitudinally over the right soleus muscle (SOL), medial gastrocnemius muscle (MG), LG, and tibialis anterior muscle, with an inter-electrode distance of 20 mm. A 16-bit A/D converter with 1-kHz sampling frequency was used to measure the above data. The measured data were stored on a personal computer.

Measured EMGs were numerically rectified. Both rectified EMGs and the kinematic data were low-pass filtered using the fourth-ordered, zero-phase-lag Butterworth filter (Winter 1990Go). Since this study mainly addressed the concordance of low-frequency body movements (a major part of the signal was in the frequency range below 1 Hz (Fitzpatrick et al. 1992bGo; Gurfinkel 1973Go; Masani et al. 2001Go), the cutoff frequency of the Butterworth filter was set to 4 Hz. The rectified and smoothed EMGs were considered to represent the level of muscle activity during quiet stance. Since the activity of the tibialis anterior muscle during the experiments was marginal, a decision was made not to consider it further in our analysis. These filtered time series were used in the following cross-correlation analysis.

Theoretical analysis

Figure 2 shows the closed-loop controller and the plant/body regulated by the PD controller adapted from Peterka (2000Go). The body dynamics and kinematics during quiet stance were described using an inverted pendulum model with parameters set to values of the typical adult male (m = 76 kg, I = 66 kgm2, and h = 0.87 m; as shown in Fig. 1). The input to the body model was the torque exerted about the ankle joint, which consisted of two components. One was due to a random disturbance torque (Td) that generated body sway patterns similar to those observed experimentally, and the other one was due to the command torque (Tc) generated by the CNS in response to body motion and disturbances. The "sensory" time delay ({tau}1 = 0.05 s), which represents cumulative time loss due to neural transmission from the ankle proprioceptors to the CNS and the sensory information processing by the CNS, was introduced. The "command" time delay ({tau}2 = 0.05 s), which represents cumulative time loss due to the CNS decision-making process and the neural transmission from the CNS to the ankle extensors, was introduced. The simulations were performed using Simulink software (MathWorks). A Gaussian noise time series with zero mean and unit variance was used to model Td. This noise was first-order, low-pass filtered with a time constant of {tau}f= 1 s, which produced waveforms with sway characteristics similar to those seen in experiments. The system's dynamics was typically simulated no longer than 40 s and the sampling frequency was 100 Hz.



View larger version (11K):
[in this window]
[in a new window]
 
FIG. 2. Close-loop control scheme of quiet stance used in the theoretical analysis.

 

As discussed above, a PD controller that used COMdis and COMvel as feedback variables was applied in simulations. The PD controller was defined with proportional and derivative gain factors, KP and KD, respectively. The command torque was then calculated according to the following equation

(4)
In our simulations, two PD controllers were considered. One was the LDG controller, and the other one was the HDG controller. In the LDG case, KP was arbitrarily set to 20N · m · deg-1 and KD was set to 4N · m · deg-1, where ratio KD/KP was 0.2. In the HDG case, KP was set to 20N · m · deg-1, the same as in the LDG case and the KD was set to 10N · m · deg-1, where ratio KD/KP was 0.5. It should be noted that the pure proportional controller (KD = 0 or KD << KP) applied to the inverted pendulum system was found to be unstable.

Before the simulations with the proposed model were carried out, a Nyquist analysis was performed to demonstrate that both the LDG and HDG controllers generated stable behavior of the given system and that they were both able to compensate for Td disturbances despite time delays {tau}1 and {tau}2. Following that, 50 simulations were carried out for each of the proposed controllers. The cross-correlation analysis was applied to the simulation outputs as discussed in the following text.

Cross-correlation analysis

Cross-correlation was applied to the following time series: 1) COPdis and each EMG, 2) COMdis and each EMG, 3) COMvel and each EMG, 4) COMdis and Tc, and 5) COMvel and Tc. Please note that an assumption was made that the measured EMG signals in experiments and the Tc obtained in simulations represent the control command of the CNS in response to system's perturbations.

To evaluate the cross-correlation and time shift between two time series x and y that had zero means, the CCF (Rxy({tau})) was applied. Rxy({tau}) was defined as follows

(5)
where {tau} denotes the time lag of y with respect to x, and the overbar denotes an average over time t. Since the direct calculation of Eq. 5 needs considerable time, the fast Fourier transform (FFT) was used instead (Bloomfield 2000Go). First, the cross-power spectral density between x(t) and y(t) was calculated using a 213-point FFT. The was then obtained by applying the inverse FFT to this cross-power spectral density function (Bloomfield 2000Go).

The following is the actual procedure that was applied to calculate the CCFs. The 30-s-long data sets were first divided into seven segments that were 213 points long, i.e., 8.192 s long. Please note that almost half of the selected 8.192-s-long segments were overlapped with the adjacent segments. Then, 13-bit FFT algorithm was applied to these segments to generate segments' periodograms. Next, an ensemble-averaged CCF of these periodograms was calculated as a CCF for each trial. Therefore the cross-correlation coefficients were calculated for 8.192-s-long time segments, and the CCF time resolution was 0.002 s (Nyquist criterion for sampling period of 0.001 s). The value and time shift of the highest peaks of the ensemble-averaged CCF of five trials for each eye condition (EO and EC) were calculated for each subject individually. Later, the group mean value was obtained for each eye condition. For theoretical data, the group mean value of 50 simulations was obtained.

The peak value of CCF and the time shift were used as variables of interest in the Results and the Discussion. Fisher's Z-transform was applied to each peak CCF value to normalize the data for subsequent statistical analysis. However, the graphs presented in the document are provided in the original form, i.e., before the Z-transform was applied, for the reader's convenience. The significance of the group mean peak CCF value and the difference of the group mean time shift with respect to zero value were tested by t-test. P < 0.05 was used as a level of significance to prevent excessive false-positive results.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 DYNAMICS OF QUIET STANCE
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 ACKNOWLEDGMENTS
 REFERENCES
 
Cross-correlation analysis for position parameters and muscle activities

Figure 3 illustrates a typical time series of the COPdis, COMdis, and COMvel, and EMG time series for SOL, MG, and LG, for the EC condition. The slow component of COPdis coincides with COMdis, while the faster component of COPdis oscillates around the COMdis. Figure 4, A and B show the CCF between COPdis and EMGs and between COMdis and EMGs, respectively, for EC condition for the same subject as in Fig. 3. The horizontal broken line indicates the r value of ±0.195 at which r is statistically different from zero (P < 0.05, n = {infty}). There is a clear positive significant peak for each EMG, which-has a negative time shift in both figures. A positive correlation indicates that, as EMG increased (decreased), the COPdis and COMdis moved forward (backward). A negative time shift indicates that the position parameters (COPdis and COMdis) lagged behind the EMGs.



View larger version (37K):
[in this window]
[in a new window]
 
FIG. 3. Representative examples of time series for a single trial in one subject during EC condition: A: COP displacement (COPdis; thin line) and COM diplacement (COMdis; bold line); B: COM velocity (COMvel); C: soleus (SOL) electromyogram (EMG); D: medial gastrocnemius (MG) EMG; E: lateral gastrocnemius (LG) EMG.

 


View larger version (19K):
[in this window]
[in a new window]
 
FIG. 4. Examples of the normalized cross-correlation function (CCF) ensemble averaged for 5 trials during eye closed (EC) condition for the same subject as shown in Fig. 3: A: CCF between COPdis and EMGs; B: CCF between COMdis and EMGs; C: CCF between COMvel and EMGs. In all 3 figures the thin line indicates CCF for SOL, the bold line indicates CCF for MG, and the dashed line indicates CCF for LG. The horizontal broken line indicates an r value of ±0.195 at which r is different from zero (P < 0.05, n = {infty}).

 

Figure 5 summarizes the peak values and time shifts of CCFs between position parameters and EMGs. Figure 5, A and B show peak values and time shifts of CCF between COPdis and EMGs, respectively. Figure 5, C and D show peak values and time shifts of CCF between COMdis and EMGs, respectively. All peak values were significant, and time shifts were significantly different from zero. These results indicate that changes of position parameters occurred with considerably long time shift, 0.147 to 0.198 s, after the corresponding muscle activities. In other words, forward displacements of COP and COM followed increasing muscle activity, and backward displacements of COP and COM followed decreasing muscle activity.



View larger version (23K):
[in this window]
[in a new window]
 
FIG. 5. Group means of the cross-correlation analysis between COP and COM measurements and EMGs: A: peak CCF values between the COPdis and EMGs; B: time shifts between the COPdis and EMGs at peak CCF values; C: peak CCF values between the COMdis and EMGs; and D: time shifts between the COMdis and EMGs at peak CCF values. For each graph, vs. SOL, vs.MG, and vs. LG indicate CCFs with muscle activity for the respective muscles. Open bar indicates CCF in eye open (EO) condition, and closed bar indicates CCF in EC condition. Data are group means ± SE (n = 16). The horizontal broken line indicates an r value of ±0.195 at which r is different from zero (P < 0.05, n = {infty}).

 

Cross-correlation analysis for COMvel and muscle activities

Figure 4C shows the CCF between COMvel and EMGs during the EC condition for the subject discussed in Fig. 3. There were two clear peaks for each EMG; one had a positive value and the other one had a negative value. The time shifts of the positive peaks were close to zero, while those of the negative peaks had larger negative values.

Figures 6A and 5B summarize the peak values and time shifts of CCFs between COMvel and EMGs, respectively. The group means of the positive peak correlations of 1) MG for both eye conditions, 2) SOL for EC condition, and 3) LG for EC condition were all significant. All group means of peak correlations for negative peaks were not significant. It is important to note that the positive peaks were significantly larger than the negative peaks for all muscles and for both eye conditions. In addition to this, we should note that the time shifts for positive peaks were much shorter than those for the negative peaks. Although all group means of time shifts for positive peaks were positive, only group means of time shifts of MG in both conditions were significantly different from zero. All group means of time shifts for negative peaks were significantly negative.



View larger version (16K):
[in this window]
[in a new window]
 
FIG. 6. Group means of the cross-correlation analysis between COMvel measurements and EMGs: A: peak CCF values between the COMvel and EMGs; B: time shifts between the COMvel and EMGs at peak CCF values. For CCF between COMvel and muscle activities, there were two peaks. One was a positive peak and the other was a negative peak. The graphs show these two peaks, respectively. For each graph, vs. SOL, vs. MG, and vs. LG indicate CCFs with muscle activity for the respective muscles. Open bar indicates CCF in EO condition, and closed bar indicates CCF in EC condition. Data are group means ± SE (n = 16). The horizontal broken line indicates an r value of ±0.195 at which r is different from zero (P < 0.05, n = {infty}).

 

Cross-correlation analysis for simulation study

The typical simulated time series of Tc and COMdis in the case of HDG are shown in Fig. 7A. One can observe that the system was stabilized by this simple controller and that the obtained COMdis time series had similar properties compared with measured COMdis time series. The CCF between COMdis and Tc was also bell-shaped and appeared to be similar to the CCF between COMdis and EMGs. Peak values of CCF between COMdis and Tc were 0.938 and 0.994 for the HDG and LDG controllers, respectively.



View larger version (24K):
[in this window]
[in a new window]
 
FIG. 7. Results of the numerical simulation experiment: A: representative examples of time series of the command torque from the controller (Tc) and the COMdis in the high-derivative gain (HDG) case; B: time shift at peak of the normalized CCF between COMdis and Tc; C: peak value of CCF between COMvel and Tc; and D: time shift at peak of CCF between COMvel and Tc. In B–D, the closed bar indicates CCF in the HDG case, and the open bar indicates CCF in the low-derivative gain (LDG) case. In C and D, the graphs are shown for the positive and the negative peaks, respectively.

 

Figure 7B shows the time shift of these CCFs. In the LDG case, the time shift was 0.059 s, which was much smaller than the actual measured value (from 0.147 to 0.198 s). However, in the HDG case, the time shift was 0.121 s, which was close to the actual measured values. Figure 7C shows peak values of CCF between COMvel and Tc, which had a two-peaks shape similar to the experimentally obtained data. In the LDG case, its positive and negative peaks had almost the same values. However, in the HDG case, the positive peak was larger than the negative peak, similar to the experimental results (see Fig. 6A). In addition, the time shift result of HDG was also similar to the experimentally obtained time shift results (see Fig. 6B). Figure 7D shows the time shifts of CCF between COMvel and Tc. In the LDG case, the time shifts of the positive and negative peaks were similar in value. However, in the HDG case, the time shift of the positive peak had a value similar to the one obtained in experiments, and the time shift of the negative peak was significantly larger compared with the positive peak.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 DYNAMICS OF QUIET STANCE
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 ACKNOWLEDGMENTS
 REFERENCES
 
Anticipatory modulation of muscle activities to COMdis

The ankle extensor activities were positively correlated with, and temporally preceded, both COPdis and COMdis. Since the finding regarding the COPdis can be accounted for by Eq. 2, only the findings pertaining to COMdis will be discussed.

The cross-correlations between COMdis and EMGs indicate that the body sway is closely related to muscle activities, meaning that the body sway reflects certain aspects of the active control mechanism of posture control (Morasso and Schieppati 1999Go). The time shift results indicate that the ankle extensors are controlled in anticipation of the change in the COMdis, which confirms findings reported by Gatev et al. (1999Go). Morasso and Schieppati (1999Go) in their theoretical analysis suggested the necessity of an anticipatory activity in the calf muscle with respect to COM position change to compensate for the neural transmission delay. Fitzpatrick et al. (1996Go) demonstrated experimentally that anticipatory action is necessary to counteract body position change. Our result and the results of Gatev et al. (1999Go) present additional evidence for these two findings. However, until now, the mechanism of how the CNS anticipates future COM position remained unclear.

Contribution of velocity information to controlling quiet stance

The numerical simulations carried out with a single joint inverted pendulum model, the balance of which was regulated with a PD controller, generated similar time shift results between COMdis and Tc as the time shift between COMdis and EMGs measured during quiet stance experiments. In particular, similar time shift results were obtained for the HDG type of controller. In the case of the LDG controller, the time shift was very different from the one observed in the experiments. This indicates that a feed-forward mechanism is not needed to regulate balance control and that the velocity provides suffi-cient information to anticipate the COM position change. In addition, the parameters of the CCF for COMvel and Tc time series obtained in simulations using the HDG controller were similar to the ones obtained during quiet standing experiments (Figs. 6, A and B and 7, C and D). The CCF between COMvel and EMGs showed that the positive peaks were larger than the negative ones and that the time shifts for positive peaks were much shorter compared with the negative ones. Similar CCF parameters were also achieved in simulation with the HDG controller, while the results obtained with the LDG controller substantially differed from the ones observed in experiments. These results indicate that the actual postural control system adopts a control strategy that relies notably on velocity information and that such a controller can modulate the muscle activity in an anticipatory manner similar to the one observed in quiet standing experiments. Furthermore, this finding shows that a feed-forward mechanism for controlling balance is not needed and that the CNS may be able to anticipate changes in COMdis if it is provided with sufficient velocity information.

Change in COMvel indicates the direction and intensity with which the current COM position will be changed in the next time instant. Thus the COMvel information is critical for anticipating how COM position will change and what corrective measures need to be taken by the CNS to compensate for these disturbances. Morasso and Schieppati (1999Go) suggested that the process in the CNS that integrates multisensory information to obtain position and velocity information of the COM (state vector) is needed to stabilize the body. In their computational simulation, the state vector was used as a feedback variable to control an inverted pendulum. Peterka (2000Go) also reported that a simple PID (proportional + integral + derivative) control, which uses the velocity information as a feedback variable, could stabilize an inverted pendulum and simulate the random walk property of the COP trajectory. These theoretical studies predicted that velocity information, in addition to position information, would contribute to the control of quiet stance. The results of this study present evidence that velocity information makes a crucial contribution to the control of quiet stance.

Because there is no sensory system that directly measures the COMvel, it is speculated that integrating multisensory information at the CNS (Horak and Macpherson 1996Go) could contribute to the velocity feedback mechanism. The strong coupling of visual information with body sway has been previously reported (Dijkstra et al. 1994aGo,bGo; Peterka and Benolken 1995Go; Schöner 1991Go). Although visual sensation is sensitive to the velocity of the visual stimulus (Dijkstra et al. 1994bGo; Schöner 1991Go), it would not be a main sensory source in the velocity feedback mechanism because it was equally observed in the EC condition. Proprioceptive sensation (Fitzpatrick and McCloskey 1994), as in group I (Griffin et al. 1990Go; Weiss and White 1986Go) and group II (Schieppati and Nardone 1997Go), is conceivable as a contributor to the velocity feedback mechanism because it plays a significant role in controlling quiet stance. Several studies reported the significant contribution of plantar cutaneous receptors in the control of quiet stance (Kavounoudias et al. 1998Go; Magnusson et al. 1990Go). Morasso and Schieppati (1999Go) pointed out the potential crucial role of these receptors in the estimation of state vector. These receptors are adequate to detect vertical and horizontal components of ground reaction force and COP position. This kinetic and kinematic information is important in calculating COM position and COMvel (Morasso and Schieppati 1999Go). Thus integrating this multisensory information, especially proprioceptive and plantar cutaneous sensations, would play a significant role in the velocity feedback mechanism. However, further investigation is needed to clarify the source and process of this mechanism.

In conclusion, we confirm the previously reported finding that activities of ankle extensors are actively modulated in anticipation of the body's position changes. By comparing the experimental results with the simulation results, we conclude that the actual postural control system during quiet stance adopts a control strategy that relies notably on velocity information and that such a controller can modulate muscle activity in an anticipatory manner without using a feed-forward mechanism.


    DISCLOSURES
 TOP
 ABSTRACT
 INTRODUCTION
 DYNAMICS OF QUIET STANCE
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 ACKNOWLEDGMENTS
 REFERENCES
 
This work was supported by the Mizuno Sports Promotion Foundation.


    ACKNOWLEDGMENTS
 TOP
 ABSTRACT
 INTRODUCTION
 DYNAMICS OF QUIET STANCE
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 ACKNOWLEDGMENTS
 REFERENCES
 
We thank T. Kimura for technical assistance with the experiments. We also thank Professors Tetsuo Fukunaga and Hiroaki Kanehisa and Dr. A. Thrasher for their advice and assistance with manuscript preparation.


    FOOTNOTES
 
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Address for reprint requests and other correspondence: K. Masani, Department of Life Sciences, Graduate School of Arts and Sciences, The University of Tokyo, 3–8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan (E-mail: masani{at}idaten.c.u-tokyo.ac.jp).


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 DYNAMICS OF QUIET STANCE
 METHODS
 RESULTS
 DISCUSSION
 DISCLOSURES
 ACKNOWLEDGMENTS
 REFERENCES
 
Accornero N, Capozza M, and Manfredi GW. Clinical multisegmental posturography: age-related changes in stance control. Electroenceph Clin Neurophysiol 105: 213-219, 1997.[Medline]

Aramaki Y, Nozaki D, Masani K, Sato T, Nakazawa K, and Yano H. Reciprocal angular acceleration of the ankle and hip joints during quiet standing in humans. Exp Brain Res 136: 463-473, 2001.[Medline]

Bloomfield P. Fourier Analysis of Time Series. Toronto: Wiley, 2000.

Collins JJ and DeLuca CJ. Open-loop and closed-loop control of posture: a random-walk analysis of center-of-pressure trajectories. Exp Brain Res 95: 308-318, 1993.[ISI][Medline]

Collins JJ and DeLuca CJ. Random walking during quiet standing. Phys Rev Lett 73: 764-767, 1994.[ISI][Medline]

Collins JJ and DeLuca CJ. The effects of visual input on open-loop and closed-loop postural control mechanisms. Exp Brain Res 103: 151-163, 1995.[Medline]

Day BL, Cauquil AS, Bartolomei L, Pastor MA, and Lyon IN. Human body-segment tilts induced by galvanic stimulation: a vestibularly driven balance protection mechanism. J Physiol 500: 661-672, 1997.[ISI][Medline]

Day BL, Steiger MJ, Thompson PD, and Marsden CD. Effect of vision and stance width on human body motion when standing: implications for afferent control of lateral sway. J Physiol 469: 479-499, 1993.[Abstract/Free Full Text]

Dietz V. Human neuronal control of automatic functional movements: interaction between central programs and afferent input. Physiol Rev 72: 33-69, 1992.[Free Full Text]

Dijkstra TM, Schöner G, and Gielen CC. Temporal stability of the action–perception cycle for postural control in a moving visual environment. Exp Brain Res 97: 477-486, 1994a.[ISI][Medline]

Dijkstra TM, Schöner G, Giese MA, and Gielen CC. Frequency dependence of the action-perception cycle for postural control in a moving visual environment: relative phase dynamics. Biol Cybern 71: 489-501, 1994b.[ISI][Medline]

Fitzpatrick R, Burke D, and Gandevia C. Loop gain of reflexes controlling human standing measured with the use of postural and vestibular disturbances. J Neurophysiol 76: 3994-4008, 1996.[Abstract/Free Full Text]

Fitzpatrick R and McCloskey D. Proprioceptive, visual and vestibular thresholds for the perception of sway during standing in humans. J Physiol 478: 173-186, 1994a.[ISI][Medline]

Fitzpatrick R, Rogers DK, and McCloskey DI. Stable human standing with lower-limb muscle afferents providing the only sensory input. J Physiol 480: 395-403, 1994b.[Medline]

Fitzpatrick R, Taylor JL, and McCloskey DI. Ankle stiffness of standing humans in response to imperceptible perturbation: reflex and task-dependent components. J Physiol 454: 533-547, 1992a.[Abstract/Free Full Text]

Fitzpatrick RC, Gorman RB, Burke D, and Gandevia SC. Postural proprioceptive reflexes in standing human subjects: bandwidth of response and transmission characteristics. J Physiol 458: 69-83, 1992b.[Abstract/Free Full Text]

Gatev P, Thomas S, Kepple T, and Halett M. Feedforward ankle strategy of balance during quiet stance in adults. J Physiol 514: 915-928, 1999.[Abstract/Free Full Text]

Griffin JW, Cornblath DR, Alexander E, Campbell J, Low PA, Bird S, and Feldman EL. Ataxic sensory neuropathy and dorsal root ganglionitis associated with Sjogren's syndrome. Ann Neurol 27: 304-315, 1990.[ISI][Medline]

Gurfinkel VS. Physical foundations of the stabilography. Agressologie 14C: 9-14, 1973.[Medline]

Gurfinkel VS, Ivanenko YP, Levik YS, and Babakova IA. Kinesthetic reference for human orthograde posture. Neuroscience 68: 229-243, 1995.[ISI][Medline]

Horak FB and Macpherson JM. Postural orientation and equilibrium. In: Handbook of Physiology. Exercise: Regulation and Integration of Multiple Systems. Control of Respiratory and Cardiovascular Systems: Bethesda, MD: Am. Physiol. Soc., 1996, sect. 12, pt. II, chapt. 7, p. 255-292.

Horak FB and Nashner LM. Central programming of postural movements: adaptation to altered support-surface configurations. J Neurophysiol 55: 1369-1381, 1986.[Abstract/Free Full Text]

Jeka JJ, Oie K, Schöner G, Dijkstra T, and Henson E. Position and velocity coupling of postural sway to somatosensory drive. J Neurophysiol 79: 1661-1674, 1998.[Abstract/Free Full Text]

Johansson R, Magnusson M, and Fransson PA. Galvanic vestibular stimulation for analysis of postural adaptation and stability. IEEE Trans Biomed Eng 42: 282-292, 1995.[ISI][Medline]

Joseph J and Nightingale A. Electromyography of muscles of posture: leg muscles in males. J Physiol 117: 484-491, 1952.[Free Full Text]

Kavounoudias A, Roll R, and Roll J-P. The planter sole is a `dynamometric map' for human balance control. Neuroreport 9: 3247-3252, 1998.[ISI][Medline]

Lacour M, Barthelemy J, Borel L, Magnan J, Xerri C, Chays A, and Ouaknine M. Sensory strategies in human postural control before and after unilateral vestibular neurotomy. Exp Brain Res 115: 300-310, 1997.[ISI][Medline]

Loram ID and Lakie M. Human balancing of an inverted pendulum: position control by small, ballistic-like, throw and catch movements. J Physiol 540: 1111-1124, 2002.[Abstract/Free Full Text]

Magnusson M, Enbom H, Johansson R, and Pyykko I. Significance of pressor input from the human feet in anterior–posterior postural control: the effect of hypothermia on vibration-induced body-sway. Acta Otolaryngol 110: 182-188, 1990.[Medline]

Masani K, Nakazawa K, and Kouzaki M. Two frequency components in ankle extensor activity during human quiet standing. Soc Neurosci Abstr 31: 305.2, 2001.

Mauritz K-H and Dietz V. Characteristics of postural instability induced by ischemic blocking of leg afferents. Exp Brain Res 38: 117-119, 1980.[ISI][Medline]

Morasso PG and Schieppati M. Can muscle stiffness alone stabilize upright standing? J Neurophysiol 83: 1622-1626, 1999.

Nashner LM. Adapting reflexes controlling the human posture. Exp Brain Res 26: 59-72, 1976.[ISI][Medline]

Nashner LM. Fixed patterns of rapid postural responses among leg muscles during stance. Exp Brain Res 30: 13-24, 1977.[ISI][Medline]

Panzer VP, Bandinelli S, and Hallett M. Biomechanical assessment of quiet standing and changes associated with aging. Arch Phys Med Rehabil 76: 151-157, 1995.[ISI][Medline]

Pavlik AE, Inglis JT, Lauk M, Oddsson L, and Collins JJ. The effects of stochastic galvanic vestibular stimulation on human postural sway. Exp Brain Res 124: 273-280, 1999.[ISI][Medline]

Peterka RJ. Postural control model interpretation of stabilogram diffusion analysis. Biol Cybern 82: 335-343, 2000.[ISI][Medline]

Peterka RJ and Benolken MS. Role of somatosensory and vestibular cues in attenuating visually induced human postural sway. Exp Brain Res 105: 101-110, 1995.[ISI][Medline]

Schieppati M, Hugon M, Grasso M, Nardone A, and Galante M. The limits of equilibrium in young and elderly normal subjects and parkinsonians. Electroenceph Clin Neurophysiol 93: 286-298, 1994.[ISI][Medline]

Schieppati M and Nardone A. Medium-latency stretch reflexes of foot and leg muscles analyzed by cooling the lower limb in standing humans. J Physiol 503: 691-698, 1997.[ISI][Medline]

Schöner G. Dynamic theory of action–perception patterns: the "moving room" paradigm. Biol Cybern 64: 455-462, 1991.[ISI][Medline]

Smith JW. The forces operating at the human ankle joint during standing. J Anat 91: 545-564, 1957.[Medline]

Weiss JA and White JC. Correlation of 1a afferent conduction with the ataxia of Fisher syndrome. Muscle Nerve 9: 327-332, 1986.[ISI][Medline]

Winter DA. Biomechanics and Motor Control of Human Movement. Toronto: Wiley, 1990.

Winter DA, Patla AE, Prince F, Ishac M, and Gielo-Perczak K. Stiffness control of balance in quiet standing. J Neurophysiol 80: 1211-1221, 1998.[Abstract/Free Full Text]

Winter DA, Patla AE, Rietdyk S, and Ishac MG. Ankle muscle stiffness in the control of balance during quiet standing. J Neurophysiol 85: 2630-2633, 2001.[Abstract/Free Full Text]

Winter DA, Prince F, Frank JS, Powell C, and Zabjek KF. Unified theory regarding A/P and M/L balance in quiet stance. J Neurophysiol 75: 2334-2343, 1996.[Abstract/Free Full Text]

Woollacott MH, von Hosten C, and Rosblad B. Relation between muscle response onset and body segmental movements during postural perturbations in humans. Exp Brain Res 72: 593-604, 1988.[ISI][Medline]

Zatsiorsky VM and Duarte M. Rambling and trembling in quiet standing. Motor Control 4: 185-200, 2000.[ISI][Medline]




This article has been cited by other articles:


Home page
J. Appl. Physiol.Home page
C. D. Tokuno, S. J. Garland, M. G. Carpenter, A. Thorstensson, and A. G. Cresswell
Sway-dependent modulation of the triceps surae H-reflex during standing
J Appl Physiol, May 1, 2008; 104(5): 1359 - 1365.
[Abstract] [Full Text] [PDF]


Home page
SURG INNOVHome page
G. Lee, T. Lee, D. Dexter, R. Klein, and A. Park
Methodological Infrastructure in Surgical Ergonomics: A Review of Tasks, Models, and Measurement Systems
Surgical Innovation, September 1, 2007; 14(3): 153 - 167.
[Abstract] [PDF]


Home page
Br. J. Sports. Med.Home page
T. Paillard, R. Bizid, and P. Dupui
Do sensorial manipulations affect subjects differently depending on their postural abilities?
Br. J. Sports Med., July 1, 2007; 41(7): 435 - 438.
[Abstract] [Full Text] [PDF]


Home page
J. Neurophysiol.Home page
G. Mochizuki, T. D. Ivanova, and S. J. Garland
Factors Affecting the Common Modulation of Bilateral Motor Unit Discharge in Human Soleus Muscles
J Neurophysiol, June 1, 2007; 97(6): 3917 - 3925.
[Abstract] [Full Text] [PDF]


Home page
ptjournalHome page
K. P Westlake, Y. Wu, and E. G Culham
Sensory-Specific Balance Training in Older Adults: Effect on Position, Movement, and Velocity Sense at the Ankle
Physical Therapy, May 1, 2007; 87(5): 560 - 568.
[Abstract] [Full Text] [PDF]


Home page
J. Physiol.Home page
I. D. Loram, P. J. Gawthrop, and M. Lakie
The frequency of human, manual adjustments in balancing an inverted pendulum is constrained by intrinsic physiological factors
J. Physiol., November 15, 2006; 577(1): 417 - 432.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.