|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1Neurological Sciences Institute, Oregon Health & Science University, Portland, Oregon; and 2Dipartimento di Elettronica, Informatica e Sistemistica, DEIS-University of Bologna, Bologna, Italy
Submitted 18 August 2004; accepted in final form 4 February 2006
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Peterka (2002)
showed that parameters of the proposed postural control model could be estimated from the anterior-posterior (AP) body-sway response to a pseudorandom stimulus that controlled the rotational motion of the SS or visual surround. Specifically, stimulus/response data were used to estimate experimental transfer functions (TFs) that characterized the postural control systems response sensitivity (gain) and timing (phase) as a function of the stimulus frequency. Then model parameter estimates were obtained from a curve fit of the models TF equation to the experimental TF data.
A key experimental result was that the TF gain decreased with increasing stimulus amplitude, but the overall shape of the TF gain and phase curves remained nearly unchanged. The model-based interpretation of this result was that a decrease in the relative contribution (sensory channel weight) of one sensory source (visual, proprioceptive, or graviceptive) was always accompanied by an increase in the contribution of other sensory sources and that other aspects of the system remained essentially unchanged. That is, a sensory reweighting mechanism is a dominant feature of postural regulation. This is in contrast to other potential interpretations where the contribution of only one sensory source might change and all other aspects of the control system remain constant, the overall stiffness (due to a combination of neurally mediated mechanisms and passive mechanics) might change to alter the amount of sway evoked by a given perturbation, or some nonlinear phenomenon might limit the amplitude of evoked sway.
The experimentally observed uniform TF gain changes imply that the amount of corrective torque generated in relation to a combined sensory error signal remains essentially constant across test conditions (Peterka 2002
). This relation between corrective torque and sensory error was referred to as the "torque normalization hypothesis, " and this hypothesis leads directly to what we refer to as the "sensory reweighting hypothesis" in the current study.
Specifically, under eyes-closed conditions with a SS perturbation, the model introduced in Peterka (2002)
indicated that the corrective torque was determined primarily by a "neural controller" acting on the value of a "sensory error signal", e, given by
![]() | (1) |
However, the second term, (Wp + Wg)·BS shows that a reduction in Wp would also cause a reduction in torque related to body sway relative to earth vertical. Modeling results demonstrate that such a reduction would produce resonant behavior and ultimately would result in instability if corrective torque was reduced below the level needed to resist gravity (Peterka and Loughlin 2004
). The disturbance caused by a tilt perturbation can be reduced and resonant behavior and instability avoided if a decrease in Wp is accompanied by a corresponding increase in Wg.
Although sudden changes in the availability of sensory information can induce transient resonant behavior (Peterka and Loughlin 2004
), experimental results obtained under steady-state conditions suggest that the postural control system is regulated in a manner that avoids resonant behavior (Peterka 2002
). This type of regulation can be represented by the equation Wp + Wg = constant. This equation is consistent with the torque normalization hypothesis when the neural controller properties and other passive mechanical factors contributing to corrective torque (associated with muscle and joint properties) do not change with stimulus conditions. Experimental results indicate that there are only small changes in neural controller and passive properties with changing stimulus amplitude (Peterka 2002
). This equation is also an explicit representation of the sensory reweighting hypothesis. Furthermore, there is no loss in generality (see METHODS) in representing the sensory reweighting hypothesis by the equation Wp + Wg = 1 where the sensory weighting factors represent the relative contributions of the two sensory systems.
A previous analysis of AP sway evoked by SS stimuli of different amplitudes yielded Wp estimates that decreased from 0.70 at the lowest SS stimulus amplitude tested (0.5° peak to peak) to 0.24 at the highest amplitude (8° peak to peak) in subjects with normal sensory function and with eyes closed (Peterka 2002
). Because this analysis assumed that the sensory reweighting hypothesis was true (i.e., Wg = 1 Wp), these results predict that Wg increased from 0.30 at the lowest stimulus amplitude to 0.76 at the highest amplitude. The primary goal of the current study is to quantitatively test that Wg increases with increasing SS amplitude in a manner predicted by the sensory reweighting hypothesis by comparing Wg measures derived from the TF analysis to Wg measures derived from an independent analysis method.
In developing an experimental test of the sensory reweighting hypothesis, we must first consider the possible sensory sources that contribute to Wg. An obvious contributor to graviception is the vestibular system, but others might include somatosensory sources related to the pressure distribution on the feet (Maurer et al. 2001
) and other nonvestibular sources (Mittelstaedt 1998
). Our previous results indicate that nonvestibular sensory sources do not contribute to graviception over the 0.05 to 1.2 Hz frequency bandwidth under our experimental conditions (Peterka 2002
).
If vestibular information is the sole source of graviceptive information under our experimental conditions, then a galvanic vestibular stimulus (GVS) can be used to test the prediction that Wg changes with SS stimulus amplitude. GVS is known to alter the activity of vestibular nerve afferents by acting directly on the nerve spike trigger zones (Goldberg et al. 1984
). Although GVS is a crude stimulus, in that it is uncertain what combination of vestibular receptor afferents are stimulated (Schneider et al. 2002
; Wardman and Fitzpatrick 2002
), it has long been known that robust postural perturbations are evoked by GVS applied through electrodes placed on the mastoid processes (Day et al. 1997
; Hlavacka and Njiokiktjien 1985
; Lund and Broberg 1983
; Nashner and Wolfson 1974
; Séverac Cauquil and Martinez 2000
). The uncertainty about which vestibular afferents are activated is not critical for the test of the sensory-reweighting hypothesis because we are primarily interested in determining whether or not the body-sway response to a given fixed intensity of the GVS current scales in accordance with the predictions of our model-based interpretation of SS responses. In particular, the sensory-reweighting hypothesis predicts that GVS responses should scale in proportion to estimates of Wg derived from responses to SS stimulation. In fact, we will go beyond this demonstration by showing that Wg estimates can be derived from GVS responses so that these GVS-derived Wg estimates can be directly compared with the Wg estimates derived from responses to SS stimulation.
The most robust responses to GVS are obtained with a bipolar-binaural electrode configuration where current is passed between the two ears through electrodes applied to the mastoid processes (Séverac Cauquil and Martinez 2000
). However, this bipolar-binaural stimulation evokes medial-lateral (ML) sway. Because measures of Wg obtained in previous experiments were based on SS stimuli that evoked AP body sway (Peterka 2002
), it was necessary to first demonstrate that analogous results were obtained from SS stimuli that evoke ML body sway.
The goal of the current study was to quantitatively test the prediction of the sensory reweighting hypothesis that a change in the proprioceptive contribution to postural control will be accompanied by an equal and opposite change in the graviceptive contribution under eyes closed conditions. To accomplish this, we first showed that ML body-sway responses to SS stimulation were analogous to our earlier AP results. We then analyzed body sway evoked by the simultaneous presentation of a pulsed GVS and a pseudorandom SS stimulus. The TF analysis method was applied to SS responses to derive estimates of Wp and thereby obtain model-based estimates of Wg = 1 Wp. A separate analysis of GVS responses provided a second estimate of Wg for comparison. Results showed a close correspondence between the two Wg measures.
| METHODS |
|---|
|
|
|---|
All subjects completed a medical screening prior to being tested to confirm that they had no history of balance problems or neurological defects and gave their informed consent prior to being tested using a protocol approved by the Institutional Review Board at Oregon Health & Science University and in accordance with the 1964 Helsinki Declaration.
Thirteen subjects (6 males and 7 females) participated in this study. The results from one subject were excluded due to extremely large amounts of oscillatory body-sway behavior, which differed considerably from the results of the other subjects. The remaining 12 subjects ages ranged from 21 to 48 yr [mean: 29 ± 7.5 (SD) yr].
Equipment
Subjects stood on a servo-controlled SS that could be rotated side to side using a servo-controlled motor. The SS rotational axis, orthogonal to the frontal plane, was centered between the subjects two feet and located 8 cm above the surface of the platform at a height approximately equal to ankle joint height.
To ensure that body movements remained similar to that of a single link inverted pendulum with sway only in a ML direction, most tests were performed with the subjects body secured to a backboard assembly with a padded headrest and straps that held the subject against padded supports at knee, hip, and shoulder levels. The backboard was free to rotate about an axis aligned with the rotational axis of the support surface. The backboard rested on its own set of bearings, so the subject did not support the backboard weight, although its mass and moment of inertia (11.4 kg and 12.9 kg · m2, respectively) were accounted for in the data analysis.
Eight subjects performed all tests with the backboard assembly. These subjects stood with their feet close together in a comfortable position, and this position was used throughout testing. The mean distance between the center of the two heels was 15 ± 2.4 cm (mean ± 1 SD). Four additional subjects performed a set of tests both with the backboard and freestanding. Their center-to-center heel distance was 15 cm on both backboard and freestanding trials.
An isolated galvanic current generator (A-M Systems analog stimulus isolator Model 2200, Carlsborg, WA) was used to pass current through two electrodes (Reusable Self-Adhering TENS/NMES/FES electrodes, 3.17 cm in diameter, # 626SS; UNI-PATCH) placed on the right and left mastoid processes in a bipolar binaural configuration. With the head facing forward on the body, this electrode configuration evoked ML body sway toward the anode (positive electrode) (Séverac Cauquil and Martinez 2000
). The impedance measured across the two electrodes was always <3.5 k
. Throughout the experiment, a voltage proportional to the galvanic current was monitored on an oscilloscope and recorded using an isolated current probe amplifier (Tektronix AM503) to ensure correct delivery of the galvanic current.
Stimulus delivery and data sampling were computer controlled (DEC LSI 11/73) using custom programs. The data-sampling frequency was 100/s. Sampled data included the SS angular position measured by a potentiometer, the backboard angular position measured by a potentiometer, backboard angular velocity measured by a rate sensor (Watson rate sensor, ARS-C141-1A, Watson Industries, Eau Claire, WI), and a voltage signal proportional to the galvanic current. Because the body was restrained to sway as a single-link inverted pendulum, the backboard angular position and velocity measures provided direct measures of the body center-of-mass (COM) angular sway and sway velocity. COM sway and sway velocity were considered to be the output variables of interest for data analysis.
In the freestanding tests, the method previously used to measure AP COM angular sway (Jeka et al. 2004
; Peterka 2002
) was applied to measure ML COM angular sway, which was considered to be the output variable of interest. Briefly, lateral body motion was measured using two horizontal rods that rotated in relation to ML displacement at shoulder and hip levels. Rotational positions of the rods were measured by two earth-fixed potentiometers. Trigonometric conversions were used to calculate lateral body displacements from the rotational positions of the rods. Based on results from a prior calibration trial, the COM lateral displacement was calculated from hip and shoulder displacements. Then the COM rotation angle was determined using estimates of the subjects COM height (based on anthropometric measures) above the ankle joint.
GVS and SS stimuli
Two different types of stimuli were used to evoke body sway. One was mechanical, the rotational motion of the SS, and the other was electrical, the GVS. SS rotation was used to provide a specific perturbation to identify input/output dynamics of postural control following the methods in Peterka (2002)
or to limit the contribution of proprioceptive cues to postural control using "sway-referencing" (Peterka and Loughlin 2004
). GVS was used to perturb the vestibular system, to determine the extent to which vestibular information was contributing to postural control in different test conditions.
We used a pseudorandom SS rotational stimulus based on a pseudorandom ternary sequence (Davies 1970
) that was similar to the one used previously to characterize the dynamics of AP-evoked sway (Peterka 2002
). The ternary sequence was mapped into a set of rotational velocities of +v, 0, or v°/s that were held throughout each of the n = 242 states of the pseudorandom sequence with a state duration of
t = 0.20 s. The period of the pseudorandom stimulus was N ·
t = 48.4 s. The velocity (v) was selected so that the mathematical integral of the pseudorandom stimulus (Fig. 1A) provided peak-to-peak rotations of 0 (control), 1, 2, 4, and 8° on different trials. The stimulus was symmetric, providing equal tilts to the right and left. Six consecutive cycles were presented during each test trial.
|
The GVS consisted of a succession of positive pulses alternating with negative pulses (Fig. 1B). The positive polarity was defined as a right-side anodal, left-side cathodal stimulation. The pulse duration was 0.5 s (50-ms rise and fall time). The pulse amplitude was 0 (control), 0.25, or 0.75 mA. The time period between the onset of each pulse was 6.05 s, whereas the time period between successive pulses of the same polarity was 12.1 s, which corresponded to one quarter of the period of the SS stimulus. A total of 24 positive and 24 negative pulses were presented during a test trial.
On some trials, the GVS and the pseudorandom SS stimulus were given simultaneously. However, these two stimuli were designed to be mathematically uncorrelated so that responses to the two stimuli could be isolated from one another and analyzed separately. Spectral analysis demonstrates that the GVS and SS stimuli are uncorrelated (Fig. 1, C and D). The pseudorandom SS stimulus has energy only at odd harmonic frequencies (dark vertical lines in Fig. 1, C and D). The pulsed GVS has energy only at even harmonic frequencies of the SS fundamental frequency (gray vertical lines in Fig. 1, C and D) where the pseudorandom SS stimulus has zero energy.
Separability of responses is guaranteed in the case of linear responses to the stimuli. For each SS stimulus amplitude, essentially linear responses to the SS stimulus were expected based on previous results from the characterization of AP sway dynamics (Peterka 2002
). Although less is known about the linearity of responses to GVS, there is some evidence that responses are quite linear for small current amplitudes (Hajos and Kirchner 1984
). However, even strong static nonlinearities, such as a rectification, in the GVS response would produce response components only at frequencies that are even integer multiples of the frequency components in the GVS. Therefore these higher frequency GVS harmonics, which are due to nonlinear responses, would still be at frequencies that are even harmonics of the fundamental frequency of the SS stimulus. Analysis of the experimental results showed that our separability assumption was justified.
Protocol
Eight of the 12 subjects performed tests consisting of various combinations of SS rotation and galvanic stimulation. All of these tests were performed using the backboard. These tests included six different SS conditions for each of the three different galvanic stimulus conditions. The six SS conditions included a fixed SS, four different amplitudes of the pseudorandom stimulus (1, 2, 4, and 8° peak to peak), and a sway-referenced SS. The GVS conditions included 0- (no-GVS), 0.25-, and 0.75-mA pulse stimulus amplitudes. Data were collected in two 2-h sessions on separate days. The test sequence was randomized to avoid habituation that might occur when similar trials were presented in succession and to reduce possible habituation to galvanic stimulation characterized by a reduction of body sway evoked by GVS. The eight subjects included in this portion of the study were able to complete all trials presented with the exception of one subject who was unable to complete the sway-referenced SS condition with 0.75-mA GVS even after many attempts.
Four subjects performed tests designed to determine if the backboard affected postural responses to the SS stimulus in the ML direction. These tests included the four different amplitudes of the pseudorandom SS stimulus (1, 2, 4, and 8° peak to peak) performed in both freestanding or backboard-supported conditions. No GVS was used in these experiments. Tests were completed in two sessions on two different days. All tests in each session were either freestanding or backboard supported, and the four SS amplitudes were presented in randomized order. Two of the subjects performed freestanding tests prior to backboard-supported tests, and the other two performed backboard-supported before freestanding tests.
All subjects were instructed to maintain a relaxed upright stance position. Subjects wore headphones and listened to audiotapes of short stories to mask equipment sounds and to maintain alertness.
Data analysis
PSEUDORANDOM SS STIMULUS/RESPONSE ANALYSIS.
The body-sway response to the SS stimulus was analyzed using methods described previously (Peterka 2002
). Spectral analysis methods were applied to sway data to compute a transfer function, relating the body-sway response to the SS stimulus, and a coherence function, indicating the degree of linearity between stimulus and response and the presence of noise in the system (Bendat and Piersol 2000
). Each transfer function was expressed as a gain and phase function over a range of stimulus frequencies. The power and cross-power spectra used to compute transfer and coherence functions were calculated using the discrete Fourier transform to decompose the SS stimulus (SS angular velocity) and COM response (backboard angular velocity) waveforms into their spectral components. For each of the last five stimulus cycles (1st cycle was excluded to avoid initial transient responses), spectral components were calculated at 150 frequencies ranging from 0.021 Hz (f1 = 1/48.4) to 3.1 Hz (f150 = 150/48.4) where 48.4 s is the period length of the SS stimulus. Because the SS stimulus has zero energy at all even harmonic frequencies of the fundamental, all even frequencies were discarded leaving 75 odd harmonic frequency samples. Power and cross-power spectra were computed by averaging the remaining 75 odd harmonic spectral components across the stimulus cycles. To further reduce the variability of the power spectra at higher frequencies, the higher frequency portion of the power spectra was smoothed by averaging spectral points at adjacent frequencies to obtain final power spectra at 17 frequencies ranging from 0.021 to 2.79 Hz and spaced at approximately equal intervals on a logarithmic frequency scale. Specifically, the lowest two odd-harmonic components were not averaged with adjacent frequencies, but adjacent odd-harmonic components at higher frequencies were averaged. The number of adjacent components included in an average increased with increasing frequency such that the highest frequency component was derived from an average of the highest 15 of the 75 odd harmonic components.
GVS ANALYSIS.
Responses to the GVS were analyzed by computing a pulse-triggered average response for both the positive and negative GVS pulses. A pulse-triggered average was computed for both COM body-sway angle and angular velocity data. Because the body-sway responses to the initial set of GVS pulses were often larger than responses later in the trial, consistent with habituation to the GVS as characterized in previous studies (Johansson et al. 1995
; Magnusson et al. 1990
), data from the first 48.4 s of each test trial, corresponding to the first cycle of the SS stimulus, were not included in the pulse-triggered average to avoid transient responses.
The pulse-triggered average was formed by averaging 6.05-s segments of body-sway or sway-velocity data after the beginning of each pulse, with separate averages formed for the responses to positive and negative pulses. Data from a total of 20 positive and 20 negative pulses were averaged to form the pulse-triggered average responses. The value of the first point in each 6.05-s segment was subtracted from all of the data samples in that segment prior to averaging so that the first sample in the pulse-triggered average was always zero. The peak amplitude and time to peak of each pulse-triggered average response was measured for both body-sway angle and angular velocity averages. The time-to-peak measures from GVS trials were used to define a window that was subsequently used to search for peak values on no-GVS control trials.
Model-based interpretation
A theoretical model of postural control was used to interpret the experimental results. The model represents a quantitative hypothesis whereby we explicitly state our assumptions about the organization of the postural control system. These assumptions can be later confirmed or refuted based on experimental results. The model in Fig. 2 is based on a modified version of the "independent channel model" described previously (Peterka 2002
) but includes a torque feedback mechanism and excludes the visual channel because eyes were closed during all trials.
|
The model assumes that the CNS is able to derive a veridical estimate of body in space from vestibular sensory receptors (otoliths and semicircular canals) and a veridical estimate of body sway relative to the feet from the various proprioceptive sensory receptors (Gardiner et al. 2000
). That is, the sensory systems represented in the Fig. 2 model do not include the complex dynamics of individual sensory receptors or the dynamics of different sensory subsystems, such as the semicircular canals, but rather represent the final outcome of the CNSs processing of motion information encoded by individual receptors. For the vestibular system, this assumption is consistent with experimental results showing that rotational information from the semicircular canals can be combined with linear acceleration and gravity signals from the otolith organs to accurately distinguish head tilt (rotation) with respect to gravity from head translation (linear motion) over a wide range of frequencies (Angelaki et al. 1999
, 2001
; Merfeld and Young 1995
; Merfeld et al. 1999
). Similarly, for the proprioceptive system, a huge array of sensory signals with various dynamic properties are available from individual sensors. However, the detailed studies of motion perception by Mergner and colleagues indicate that subjects have a veridical perception of body kinematics derived from proprioception over a wide range of frequencies (Mergner et al. 1983
, 1991
, 1993
). Additional neurophysiological evidence (Bosco and Poppele 2000
; Bosco et al. 2000
) is available indicating that sensory neurons in the CNS encode more global parameters of limb movement and posture rather than specific muscle and joint parameters.
There is some evidence that force sensors, such as Golgi tendon organs or pressure sensors in the feet, contribute to postural control (Maurer et al. 2001
), but this influence seems to be limited to low-frequency corrections (Mergner et al. 2003
; Peterka 2003
). The Fig. 2 model uses force sensors to detect the corrective torque, Tc, applied to the body. The model assumes that the torque signal is low-pass filtered (1st-order filter with gain, KT, and time constant,
T), fed back, and summed together with the graviceptive and proprioceptive signals to form the sensory error signal, e. The sign of the torque feedback is positive. As discussed in Peterka (2003)
, this form of positive torque feedback provides a mechanism that corrects for low-frequency deviations of body orientation away from vertical and accounts for the low-frequency phase lead and gain decline observed in experimental TFs. When torque feedback is included in the model, it is no longer necessary to include an integral control factor in the neural controller, as used in previous models (Peterka 2002
; Peterka and Loughlin 2004
), to account for low-frequency data. Others have suggested roles for positive force feedback in motor control systems (Duysens et al. 2000
; Geyer et al. 2003
; Prochazka et al. 1997a
,b
).
Inclusion of a torque feedback mechanism does not alter the basic concept behind the torque normalization hypothesis. Torque feedback influences low-frequency behavior and provides a relatively small contribution to e in comparison to the proprioceptive and graviceptive contributions. In contrast, the proprio- and graviceptive contributions influence the dynamic behavior of the postural control system over its entire bandwidth. The experimental results support the limited influence of torque feedback. Eyes-closed subjects with absent vestibular function fall without evidence for any corrective action when the SS is sway-referenced (Nashner et al. 1982
) suggest that force feedback is a minor contributor to stance control. In addition, the analysis of responses to pseudorandom SS rotations in eyes-closed vestibular-loss subjects identified Wp values of 1.0 (Peterka 2002
), suggesting that these subjects were relying entirely on the wide bandwidth information available from proprioceptive sensory cues to maintain stability and that the torque feedback only modulated the low-frequency behavior.
The weighted sum of vestibular, proprioceptive, and torque feedback signals, represented by e in Fig. 2, results in the generation of Ta via the action of a "neural controller. " The active torque generated by the neural controller includes two components, which are respectively proportional to e (gain constant KP) and to the rate of change of e (gain constant KD). A time delay,
d, accounts for neural transmission, processing, and muscle activation and torque generation delays.
For ML body sway with the feet placed at approximately hip width, the pelvis and legs move as a parallelogram with the pelvis remaining parallel to the SS, and the upper body, when constrained by a backboard assembly, tilts with respect to the pelvis when the body is not aligned with the SS. These mechanics are more complex than a simple inverted pendulum. Also in contrast to a simple inverted pendulum, where torque is applied at one joint, torque can be applied at a combination of ankle and hip joints to control lateral body sway. However, with regard to the control of COM motion, Winter et al. (1998)
showed that for small sway amplitudes, the equations of motion that govern a simple inverted pendulum are the same as those for the more complex lateral sway mechanics. Therefore for our study, the body was modeled as a linearized single-link inverted pendulum (Eq. 11, Peterka 2002
).
Estimates of the sensory channel weights as well as other postural control parameters were obtained by fitting the experimental transfer function data to the model equation given by
![]() | (2) |
ds represents the time delay, P = K + Bs models the passive mechanical properties, and T = KT/(
Ts + 1) models the torque feedback low-pass filter. The feet remain in contact with the support surface, therefore FS = SS. The biomechanics of the human body are represented by J, m, h, and g, which are, respectively, body moment of inertia about an axis at ankle joint level, body mass (excluding mass of the feet), COM height of the body above the ankle joints, and gravitational acceleration. The mass m was determined by directly weighing each subject, subtracting 2.6% of body mass to account for the mass of the feet (Diffrient et al. 1974
The curve fits of Eq. 2 to the experimental TF data were performed using a constrained nonlinear optimization implemented in the Matlab Optimization Toolbox (The MathWorks, Natick, MA) to minimize an error function based on the normalized difference between model and experimental TF data (Eq. 9 in Peterka 2002
). We found that Eq. 2 provided a good description of TF data in the low- and mid-frequency range of 0.021 to 1.47 Hz, but data above this frequency range were not always adequately described. Specifically, the higher frequency TF data showed greater phase lag than the fit. When TF data at all frequencies were included in the fit procedure, the fit typically showed resonant properties in the mid-frequency range that clearly deviated from TF data in this region. Therefore we chose to focus on the mid-frequency region that appears to be particularly affected by interactions associated with sensory integration. Accurate fits of Eq. 2 to the low- and mid-frequency TF data were obtained by excluding the highest four TF frequencies.
We tried variations of the model to improve the quality of the fit over the full range of frequencies. The most successful variation was to add an acceleration component to the neural controller. This addition accounted for the high-frequency phase but also resulted in unrealistically large time delay estimates. Because the inclusion of the acceleration component did not noticeably change the estimates of parameters other than the time delay, the conclusions of this study regarding sensory reweighting were not affected by the use of the simpler model without an acceleration component. It is likely that higher frequency dynamics are influenced by several mechanisms that are not fully accounted for in our model. These could include a more complex neural controller that includes an acceleration component, separate neural controllers for each sensory channel, different time delays in each sensory channel, shorter latency contributions from stretch reflexes (Mirbagheri et al. 2000
), muscle activation dynamics, and dynamics of the muscle/tendon system (Zajac 1989
). Future work is needed to develop a comprehensive model that can explain the full range of data.
Without loss of generality, we define the sensory channel weighting factors, Wg and Wp, so that they represent the relative contributions to postural control of the graviceptive and proprioceptive systems, respectively. That is
![]() | (3) |
The quality of the model fits was assessed by using model parameters derived from each test trials TF data to simulate (using Matlab Simulink) the steady-state time course of the response to the SS stimulus for that trial. The input to the model was the SS angular tilt recorded during the test trial. The output was the simulated body sway angle response. The "residual sway" over one stimulus cycle was calculated by subtracting the mean experimental body sway from the simulated body sway. The proportion of the variance accounted for by the model was determined by taking the difference between the mean squared value of the experimental body sway and the mean squared value of the residual sway and then normalizing by the mean squared value of the experimental body sway. The expectation was that on trials with combined GVS and SS stimuli, the residual sway traces should reveal the isolated responses to the GVS.
If the contribution of passive components (K and B) to the generation of corrective torque is small relative to contribution of active components and if the time delay and neural controller parameters do not change across stimulus conditions, then Eq. 2 shows that Wp is the major determiner of the overall gain of the transfer function and that changes in the value of Wp will only minimally influence the overall shape of the gain or phase functions. Therefore experimental conditions that yield higher transfer function gains without changing the transfer function shape are associated with greater utilization of proprioceptive cues for postural control (i.e., larger values of Wp). Because the sensory channel weights represent the relative contribution of each sensory system to postural control, a larger Wp value implies a smaller Wg value because Wg = 1 Wp in our hypothesized formulation of the sensory reweighting hypothesis (Eq. 3).
The Fig. 2 model includes a switch for turning on sway-referencing. When this switch is closed, the body-sway angle, BS, is used to control the angle of the support surface. For ideal sway-referencing, SS = BS, and therefore the body-on-foot angle, BF, which is sensed by the proprioceptive system, is always zero. Therefore information from the proprioceptive system should not contribute to the generation of corrective torque in the sway-referenced condition. However, information from the vestibular system still can contribute to the generation of corrective torque, and, in fact, subjects with normal vestibular function are able to maintain their balance during sway-referencing (Peterka and Black 1990
), but subjects with a bilateral vestibular loss cannot (Nashner et al. 1982
). Based on these considerations, we assume that to generate sufficient corrective torque, Wg must be increased in the sway-referenced condition compared with other conditions where proprioceptive cues contribute to postural control, to compensate for reduced availability of reliable proprioceptive orientation cues. To obtain good dynamic control over balance, it is reasonable to assume that Wg = 1 in the sway-referenced condition (Peterka 2002
; Peterka and Loughlin 2004
).
The model includes a GVS input that is assumed to perturb the afferent activity of vestibular receptors leading to an altered internal vestibular-based estimate of body-in-space position, bs + bsgalv. The model predicts that the body-sway response to the GVS should be largest in the sway-referenced condition where we assume Wg = 1 and should decrease in proportion to Wg in other conditions where proprioceptive cues are also contributing to postural control in which we assume that Wp > 0 and Wg = 1 Wp.
Based on the model predictions, a method was devised to calculate Wg from responses to GVS. Specifically, with the assumption that Wg = 1 in the sway-referenced condition and that neural controller parameters are constant across test conditions, an "uncorrected" GVS-Wg estimate in a particular test condition was calculated by dividing the peak GVS response in that condition, obtained from the pulse-triggered average, by the peak GVS response in the sway-referenced condition (Cenciarini 2001
). However, because passive mechanical properties do not contribute to the corrective torque in the sway-referenced condition, the uncorrected GVS Wg would be in error if the postural control system does not compensate for the loss of passive torque in the sway-referenced condition. Through simulation studies, we found that the difference between the actual Wg and uncorrected Wg, normalized by the actual Wg, varied systematically as a linear function of both K and B over the range of K and B values determined from the experimental data. Regression fits were used to calculate the coefficients of this linear function using data from simulation results. Separate regressions were used to fit data for Wg values obtained from peak position and peak velocity responses to GVS. These linear functions were solved to calculate the "corrected GVS-Wg" for the given values of K, B, and the uncorrected Wg obtained from each test trial. Comparisons were then made between both the uncorrected and the corrected GVS-Wg estimates and the TF-Wg estimates obtained from the same trial data. Favorable comparisons would support the Fig. 2 model and the assumptions made in determining values of Wg from both GVS and SS stimulation.
| RESULTS |
|---|
|
|
|---|
|
The bottom trace in each pair of traces in Fig. 3, B and C, shows the residual sway after subtracting the sway predicted from model simulations using parameters obtained from fits to TF data for each trial. The residual sway traces for trials shown in Fig. 3B demonstrate that the model accounted for most of the variability in the sway responses apart from some higher-frequency variability that likely remained because the TF fit did not fully account for the higher-frequency TF data. The average of residual error variance across all subjects and all no-GVS trials showed that the model accounted for 92 ± 4% of the sway variance.
The residual sway traces for trials shown in Fig. 3C demonstrate that the model accounted for most of the variability in the sway responses to the SS stimulus. The residual sway traces clearly revealed the individual responses to the GVS. The general trend is evident for an increasing GVS response amplitude with increasing SS stimulus amplitude.
Transfer functions characterizing responses to SS stimulation
EXAMPLE TRANSFER FUNCTIONS. Figure 4 shows two example transfer functions calculated from body-sway responses of one representative backboard-supported subject to a 2° peak-to-peak pseudorandom SS stimulus. The TF in Fig. 4A was calculated from a trial where no GVS was presented, whereas the TF in Fig. 4B was from a trial with a simultaneous presentation of a 0.75-mA pulsed GVS. There was no obvious distortion of the transfer function gain and phase data caused by the simultaneous presentation of the GVS, and no clear difference between TF gain and phase results obtained from trials both with and without GVS. The general pattern of the gain and phase changes as a function of stimulus frequency was similar to those shown in Fig. 4 for all subjects and test conditions.
|
). TF data > 2 Hz were sometimes not well explained by the model, suggesting that additional dynamic components, which are not included in the Fig. 2 model, are required to account for the higher frequency postural dynamics.
MEAN TRANSFER FUNCTIONSBACKBOARD VS. FREESTANDING.
Figure 5 shows two sets of the mean gain, phase, and coherence functions for the four subjects who performed both backboard-supported and freestanding trials in response to 4 SS amplitudes. These results are similar to results comparing TFs characterizing AP-evoked sway in freestanding and backboard-supported conditions (Fig. 7 in Peterka 2002
). The mean ratio of freestanding to backboard TF gains across all frequencies and for all stimulus amplitudes was 0.89, suggesting that subjects did not fully compensate for the added mass and moment of inertia of the backboard. However, the general pattern of TF gain decrease with increasing stimulus amplitude was similar in the backboard-supported and freestanding trials.
|
|
MEAN TRANSFER FUNCTIONS ACROSS GVS CONDITIONS.
Figure 6 shows a summary of the mean gain, phase, and coherence functions for the eight backboard-supported subjects in the three GVS conditions at four SS amplitudes. For all GVS conditions, the TF gain functions at frequencies below
1.6 Hz were largest for the 1° SS stimulus amplitude and progressively decreased with increasing stimulus amplitude. At frequencies above
1.6 Hz, TF gain values were indistinguishable from one another. For all GVS conditions, the family of phase functions converged at lower frequencies (below
0.1 Hz) but showed increasing separation with increasing stimulus frequency with phase functions from larger-amplitude SS trials showing less phase lag (less negative) at a given frequency.
|
Coherence function values of unity are indicative of a linear stimulus/response relationship with no noise in the system. The fact that coherence values were large and did not change appreciably with SS stimulus amplitude suggests that the postural response to SS stimulation was well characterized by linear system behavior for each SS stimulus amplitude. However, the overall behavior of the postural control system was nonlinear because gain declined and phase showed some changes with increasing SS amplitude. Furthermore because simultaneous presentation of the pulsed GVS did not alter the qualitative pattern of responses to SS stimulation, GVS can be used to test the sensory-reweighting hypothesis.
MODEL TRANSFER FUNCTION PARAMETERS. For the eight backboard-supported subjects who performed trials that included GVS, Fig. 7 shows the mean values of the TF parameters as function of the SS stimulus amplitude for each of the GVS conditions. The effects on each parameter of the SS stimulus amplitude and GVS condition were tested using repeated-measures ANOVA and the Tukey-Kramer multiple comparisons test.
Statistical analysis revealed that KP, K, B, and
d were not significantly affected by changes in the GVS condition (P > 0.2). The parameter KD was significantly lower in the 0.25-mA GVS condition (P < 0.05) compared with the 0.75-mA GVS condition, although the difference in mean KD values between GVS conditions was very small. The parameters KP, KD, and
d showed significant changes with increasing SS stimulus amplitude (P < 0.003). KP and KD increased, whereas
d decreased with increasing SS stimulus amplitude.
The median values of torque feedback parameters, KT and
T, are also shown in Fig. 7. The distributions of these parameters were highly skewed with half the values of KT between 0.02 and 0.09°/Nm and of
T between 2 and 12 s. These parameters exerted their largest influence on the TF dynamics at frequencies below
0.1 Hz where they produced a gain decline and phase advance. The limited amount of low-frequency TF data likely reduced the accuracy with which these parameters were identified, resulting in the skewed distributions.
The proprioceptive channel weighting factor, Wp, showed a large and significant decrease with increasing SS stimulus amplitude. Post hoc analysis revealed that Wp was significantly different at each SS amplitude from the values at the other SS amplitudes. This effect was not found for any of the other model parameters. The mean Wp (across all subjects and GVS conditions) was 0.76 at the 1° SS amplitude and decreased to 0.36 at the 8° amplitude. Also evident in Fig. 7, the mean Wp values obtained during simultaneous presentation of the 0.75-mA GVS were consistently larger than the mean Wp values in the no-GVS or the 0.25-mA GVS conditions, although the difference was significant only between the 0.25- and 0.75-mA GVS conditions (P < 0.05). There was no significant difference in Wp values between the no-GVS and 0.25-mA GVS conditions. Across all SS amplitudes, the mean value of Wp = 0.62 in the 0.75-mA GVS condition was on average 10% greater than in the no-GVS or 0.25-mA GVS condition.
Figure 8 shows mean Wg values based on Wp estimates from TF curve fits and the relation Wg = Wp 1. The sensory-reweighting hypothesis will be supported if GVS responses and corresponding GVS-Wg estimates show a similar dependency on SS amplitude.
|
Figure 9 shows examples of responses to positive GVS pulses determined by calculating a pulse-triggered average from sway position data (Fig. 9A) and from sway-velocity data (Fig. 9B) from a representative subject. The 0.75-mA pulsed GVS (Fig. 9C) was delivered concurrently with a 2° pseudorandom SS stimulus. In all cases, the position pulse responses to the GVS had a characteristic peak that was directed toward the anodal electrode. A second peak, smaller and in the opposite direction, sometimes followed the first peak in agreement with previous studies (Hajos and Kirchner 1984
). The velocity pulse response had a biphasic waveform with the first peak always toward the side of the anode. Responses to positive and negative pulses were similar to each other with only the direction differing.
|
The time to the first peak was, on average, 1.26 ± 0.14 and 0.83 ± 0.08 (SD) s for position and velocity responses, respectively. The peak amplitude of GVS responses increased with increasing SS amplitude. Body-sway responses from a representative subject to a 0.75-mA positive galvanic pulse are shown in Fig. 10 during all of the SS conditions tested. The body-sway response to the galvanic stimulus was smallest for a fixed surface, increased with increasing pseudorandom SS amplitude, and was largest in the sway-referenced SS condition. The peak amplitudes measured from body sway-velocity responses showed an analogous trend with changes of SS condition.
|
The summary of GVS peak position and GVS peak velocity responses is shown in Fig. 11 for responses averaged across subjects at each test condition. The peak responses were smallest when subjects stood on a fixed SS (0°), and both position and velocity peak responses increased with increasing SS stimulus amplitude. Peak responses were consistently larger in the sway-referenced SS condition than in any other condition for a given GVS amplitude.
|
|
Comparisons of both corrected and uncorrected GVS-Wg estimates with corresponding TF-Wg estimates are shown in Fig. 13. Paired comparisons of both the corrected and uncorrected GVS- and TF-Wg estimates showed high correlations ranging from 0.58 to 0.90 with the correlations being the same or slightly greater for the uncorrected compared with the corrected GVS-Wg data (Fig. 13). Correlations were lower for the comparisons based on 0.25-mA GVS responses than for 0.75-mA GVS responses, consistent with higher variability in GVS responses to lower-amplitude GVS pulses. Correlations were higher for GVS-Wg estimates based on peak velocity than for estimates based on peak position.
|
Paired t-test showed no significant differences (P > 0.19) between uncorrected 0.25-mA GVS-Wg values and TF-Wg values for the GVS-Wg values derived from both peak position and peak velocity measures. However, there were significant differences (P < 0.003) between uncorrected 0.75-mA GVS-Wg and TF-Wg values. The GVS-Wg values, based on peak position and peak velocity measures, were on average 0.07 larger than the TF-Wg values.
Linear regression analyses were applied to uncorrected GVS-Wg and TF-Wg estimates using GVS-Wg estimates derived from both GVS amplitude conditions (0.25 and 0.75 mA) and from sway-position (Fig. 13A) and sway-velocity data (Fig. 13B). The regression analyses were performed using the paired-Wg estimates from each individual subject for each SS stimulus condition even though Fig. 13 shows only the mean (±SE) uncorrected GVS-Wg and mean (±SE) TF-Wg estimates obtained at the four pseudorandom SS amplitudes. For three of the four sets of data shown in Fig. 13, the regression lines that best fit the data had intercepts that were not significantly different from 0 (P > 0.12), and the 95% of confidence limits on the linear regression slopes included the slope of 1. The exception was the Wg comparison based on the 0.25-mA GVS sway-position data. Considering the results from both the paired comparisons and linear regressions, the closest match between TF-Wg and uncorrected GVS-Wg estimates was obtained for the 0.25-mA GVS peak-velocity data.
| DISCUSSION |
|---|
|
|
|---|
A number of assumptions and practical considerations were required to develop a protocol for testing the sensory-reweighting hypothesis. These include the applicability of the sensory-reweighting hypothesis to the control of ML sway for which GVS evokes a robust sway response, the separability of sway responses to concurrently applied GVS and support surface stimuli, and the validity of the assumptions made in the calculation of graviceptive weights from responses to GVS and SS stimuli. These considerations are discussed in the following text and are followed by additional discussion of interaction effects between GVS and SS stimuli, and a comparison of our results with previous evidence of sensory reweighting.
ML postural dynamics resemble AP dynamics
Results shown in Figs. 38, compared with the AP results from a previous study (Figs. 3, 5, 6, 10, and 14 in Peterka 2002
), confirmed that ML sway evoked by SS tilts had properties similar to AP-evoked sway. Specifically, ML body-sway responses to increasing amplitudes of SS stimuli showed a saturation phenomenon (Fig. 3). Like the AP responses, this saturation was not likely due to some static nonlinearity, which would be expected to cause a reduction of coherence function values at increasing SS amplitudes. Results showed that the coherence functions remained essentially unchanged or even increased slightly with increasing SS amplitude (Fig. 6). This implies that the stimulus/response behavior was essentially linear at each SS amplitude tested even though the overall behavior was clearly nonlinear because the overall responsiveness to the SS stimulus declined with increasing SS amplitude. Additionally, the overall shape of ML transfer function gain and phase curves (Figs. 46) were similar to those derived from AP-evoked sway and showed only small changes in shape with changing SS amplitude.
Curve fits of Eq. 2 provided a very good description of the data over a frequency range of 0.021 to 1.47 Hz (Fig. 4). The values of the model parameters derived from these curve fits showed systematic changes with SS amplitude (Fig. 7) that were very similar to those previously shown for AP-evoked sway (Fig. 10 in Peterka 2002
). A minor difference between ML and AP results was that the mean value of Wp was
37% higher for ML than AP results, indicating that subjects tend to make greater use of proprioceptive cues for controlling ML sway than for AP sway under the conditions tested. A greater reliance on proprioceptive information is potentially consistent with proprioceptive cues being more sensitive and accurate for ML compared with AP sway due to the involvement of a greater number of joints and muscles during ML sway. Day et al. (1993)
showed that ML sway involves both ankle and hip joints, and the hip joint motion increases with increasing stance width. They speculated that the reduction in GVS-induced sway observed at greater stance widths could be due to increased reliance on proprioceptive information because this information becomes more accurate when more joints and larger joint motions are involved in the sway.
The similarities between the results of ML and AP sway evoked by pseudorandom SS stimulation suggest that there are similar underlying mechanisms involved in the control of balance in both directions of sway. Hence we concluded that the original sensory-reweighting hypothesis developed from AP body-sway experiments also applies to postural mechanisms controlling ML sway.
Separability of SS and GVS responses
Testing of the sensory-reweighting hypothesis relied on our ability to separate responses to simultaneously presented GVS and SS stimuli. The stimuli were specifically designed to be mathematically uncorrelated (Fig. 1) to facilitate later separation of the responses to each stimulus component. For uncorrelated stimuli, the superposition principle (Bendat and Piersol 2000
) guarantees that responses to the two stimuli are separable if the response to each stimulus component is linearly related to that component.
Our results and the results of others suggest that there is good reason to expect linear responses from GVS and SS stimuli. A previous study found an approximately linear dependence on galvanic current amplitude and the amount of GVS-evoked body sway (Hajos and Kirchner 1984
). The high coherence function values derived from ML-evoked sway responses to SS stimuli (Figs. 5 and 6) suggest that the linearity assumption is justified for ML sway evoked by each amplitude of the SS stimulus. Results of the analyses of sway responses to simultaneously presented GVS and SS stimuli and comparisons to control-trial results (no-GVS condition) suggest that the responses to GVS and SS stimuli were separable. There were no obvious distortions of TF estimates due to simultaneous presentation of the GVS (Figs. 4 and 6). When the pulse-triggered averaging for GVS analysis was applied to SS stimulation trials performed without GVS, there was no false pulse response detected even when the SS amplitude was large (control trials responses in Figs. 10 and 11).
Estimating Wg from GVS responses
Previous experimental results suggest that the graviceptive contribution to postural control increases in the sway-referenced condition to ensure that sufficient corrective torque is generated to maintain balance (Peterka and Loughlin 2004
). If we assume that Wg = 1 during sway-referencing, the model predicts that the amplitude of GVS-evoked sway in the sway-reference condition should be larger than in any other condition and should be indicative of the full contribution to body sway caused by the GVS-induced perturbation to vestibular afferent firing rate. Then in any other condition, the reduction in the amplitude of GVS-evoked sway is indicative of a reduction in Wg, and "uncorrected" Wg estimates can be obtained by taking the ratio of the GVS pulse-response amplitude in that condition to the pulse-response amplitude in the sway-referenced condition.
The potential error in this method for estimating Wg is that muscle/joint mechanics contribute passive torque in all conditions except the sway-referenced condition. In the sway-referenced condition, the surface tracks the rotational body-sway motion, and there is no change in body orientation relative to the surface over time. Because BF = 0 in the model, the model predicts that the muscle/joint mechanics make no passive torque contribution. If the nervous system does not compensate for the loss of the passive torque component during sway-referencing, the GVS response amplitude in the sway-referenced condition will be too large, and the uncorrected GVS-Wg estimates will be too small.
Alternatively, if the nervous system does fully compensate for the loss of passive torque in the sway-referenced condition, then the "uncorrected" GVS-Wg estimates would provide an accurate measure of Wg, whereas the corrected measure will be too large as also shown in Fig. 13. Such a compensation would be consistent with the torque normalization hypothesis (see INTRODUCTION) and could occur by increasing the neural controller gain. We currently have no method for measuring the neural controller parameters in the sway-referenced condition to determine whether there was a compensation for the loss of passive torque. However, the results in Fig. 7 show small increases in the neural controller gain factors KP and KD with increasing SS amplitude, suggesting that the nervous system is able to adjust the neural controller gain. Furthermore, the uncorrected GVS-Wg estimates were more closely matched to the TF-Wg estimates than were the GVS-Wg estimates corrected for the assumed loss of passive torque in the sway-referenced condition (Fig. 13). A parsimonious interpretation of these results is that the torque normalization does occur in the sway-referenced condition, and therefore the uncorrected Wg estimates are more accurate than the corrected Wg estimates.
Even though the uncorrected GVS-Wg estimates were closer to the TF-Wg estimates than the corrected GVS-Wg estimates, there were still some systematic differences between the uncorrected 0.75-mA GVS-Wg estimates and TF-Wg estimates (Fig. 13, right). Specifically, the 0.75-mA GVS-Wg estimates were slightly larger than the corresponding TF-Wg estimates. Such a discrepancy could occur if the peak response to the 0.75-mA GVS obtained on the sway-referenced SS was biased toward a lower value. This may have occurred for responses to 0.75-mA, but not 0.25-mA pulses, because responses to the 0.75-mA pulses in the sway-referenced condition brought subjects close to the limits of stability. Indeed one subject was unable to complete a full sway-referenced trial with the 0.75-mA GVS. Such a large response may have triggered the participation of other nonvestibular sensory systems, such as somatosensory cues signaling excessive sway based on pressure cues from the feet or loading/unloading cues in the legs (Dietz 1998
; Maurer et al. 2001
) and therefore reduced the response to the 0.75-mA GVS. Sway evoked by the 0.25-mA GVS did not bring subjects close to their stability limits and therefore would not have triggered the participation of nonvestibular sensory systems.
An additional assumption in the GVS-Wg calculation is that the GVS-induced perturbation of the vestibular afferent firing rate remains the same for a given GVS pulse amplitude so that changes in responsiveness to GVS are associated with centrally mediated reweighting and not with habituation in the peripheral vestibular nerve. This assumption is supported by neurophysiological experiments that show only minimal evidence for habituation of afferent discharge rate changes during application of GVS (Courjon et al. 1987
; Goldberg et al. 1982
, 1984
).
The close correspondence between Wg estimates from responses to GVS and from the TF analysis supports the idea that all of the assumptions required for the estimate of Wg from responses to GVS and from TF analysis are, at least to a good approximation, correct. Of course, we cannot rule out the possibility that there are offsetting errors in these two analyses. However, a parsimonious conclusion is that our model description (Fig. 2) and the associated assumptions about how sensory channel weights change to facilitate stance control in changing conditions do explain important features of the human postural control system.
Interactions between SS and GVS responses
The primary interaction of interest in our study was the effect of SS stimulus amplitude on the response to the GVS. The qualitative prediction of the sensory-reweighting hypothesis was that the body-sway response to the pulsed GVS should increase with increasing SS amplitude. This prediction was confirmed (Figs. 10 and 11). Additionally, the model-based interpretation of experimental results allowed a quantitative test of the sensory-reweighting hypothesis by comparing estimates of Wg from the TF analysis of responses to SS stimulation with Wg estimates from the GVS analysis. The comparison of TF Wg estimates with GVS Wg estimates (Fig. 13) showed a close correspondence and therefore quantitatively confirms the sensory-reweighting hypothesis.
A secondary interaction between SS condition and GVS responses was also observed. The estimates of Wp derived from TFs during concurrent application of the 0.75-mA GVS were larger than Wp estimates from the no-GVS or 0.25-mA GVS conditions (Fig. 7). No other model parameter showed any systematic influence of concurrent GVS presentation. The larger Wp values mean that subjects relied more on proprioceptive cues and less on graviceptive cues during trials that concurrently presented the larger amplitude GVS. Such a result indicates that the GVS stimulus produced a change in the responses to SS stimulation. A change toward increased reliance of proprioceptive information is consistent with the general predictions of a recently developed optimal control model of postural control (van der Kooij et al. 1999
, 2001
) and is consistent with other motor-control models (van Beers et al. 1999
, 2002
). These models predict that in obtaining an optimal internal estimate of important system variables, such as the position and velocity of body or limb movement, the nervous system weights sensory information in inverse proportion to the amount of variability (noise) in the sensory-channel signals. Therefore adding significant "noise" to the vestibular channel by applying the larger amplitude GVS caused a detectable shift toward greater use of proprioceptive information. Similarly, increasing the SS stimulus amplitude effectively increases the variability in the proprioceptive channel, producing a reduction in use of proprioceptive cues for postural control with increasing SS amplitude.
Comparison with previous studies
The idea that sensory reweighting and sensory substitution contribute to human postural control is not new. However, consideration of the feedback model structure (Fig. 2) and the torque normalization hypothesis provides a quantitative understanding of the benefits of sensory reweighting. Specifically, sensory reweighting affords the ability to reduce the effect of an external perturbation while preserving stability by maintaining a constant amount of corrective torque in relation to sway about the upright position (see Eq. 1).
Previous studies have identified many qualitative examples that are consistent with the idea that the nervous system alters the sources of sensory orientation information used for balance control as environmental, biomechanical, and pathological conditions change. Studies have shown increases in GVS-evoked responses primarily under three conditions: tasks performed with limited visual and/or proprioceptive cues (Day and Bonato 1995
; Day et al. 1993
; Fitzpatrick et al. 1994
; Nashner and Wolfson 1974
; Welgampola and Colebatch 2001
), when subjects stand on a moving surface (Inglis et al. 1995
), and when pathology limits the somatosensory contribution to postural control (Day and Cole 2002
; Horak and Hlavacka 2001
).
It is clear from several studies that vestibular information is most relevant to the postural control system when it is critical for maintaining balance (Day et al. 2002
; Fitzpatrick et al. 1994
; Peterka 2002
; Welgampola and Colebatch 2001
). That is, in conditions where information from other sensory systems is compromised, a sensory reweighting occurs, causing an increased reliance on vestibular information while decreasing reliance on other sensory information. The increased reliance on vestibular information is supported by the increased responsiveness to GVS (Fitzpatrick et al. 1994
; Welgampola and Colebatch 2001
) and by decreased responsiveness to other sensory cues (Ivanenko et al. 1999
). In fact, when all three primary sensory inputs (visual, proprioceptive, and vestibular) are available and are providing veridical cues regarding body orientation in space, it is generally accepted that the vestibular system contributes less than the visual or proprioceptive systems. For example, using very-small-amplitude perturbations from SS or visual surround motions that were typically not perceived by the subjects, Peterka (2002)
showed that the proprioceptive, visual, and vestibular systems contributed
50, 35, and 15%, respectively, to quiet stance control. Peterka (2002)
also demonstrated experimental results consistent with sensory reweighting for a wide range of conditions that manipulated visual, proprioceptive, and combinations of these two orientation cues.
Others have considered models and mechanisms to account for sensory interactions for postural control. Various related models have been used to understand the influence of different perturbations (visual, light touch on a moving surface, and combinations of the two) on stance control (Dijkstra et al. 1994
; Giese et al. 1996
; Jeka et al. 1998
, 2000
; Kiemel et al. 2002
; Oie et al. 2002
). These models differ from the Fig. 2 model in that they do not consider the influence of time delay in the feedback system, and they characterize the effect of a perturbation on a system that is a priori stable. In contrast, the stability of the Fig. 2 model depends on the regulation of the sensory integration process represented by the selection of sensory channel weights.
Building on earlier studies of the perception of head/neck movements (Mergner et al. 1983
, 1991
, 1993
), Mergner and colleagues have developed a model that accounts for body sway evoked by sinusoidal external force perturbations in different surface conditions (Mergner et al. 2003
). This model includes multiple nonlinear threshold elements that alter the contribution of different sensory systems to corrective torque generation as a function of stimulus conditions and body sway. Effectively, these thresholds provide a mechanism for changing weights, and therefore there is some resemblance to the Fig. 2 model.
Optimal control models have also been developed recently that include high dimensionality and predictive control (van der Kooij et al. 1999
, 2001
). In contrast, the Fig. 2 model is of low dimensionality and does not include prediction. Prediction could play a role in the responses to our pseudorandom stimulus because the stimulus is, in fact, deterministic. Even if it was not, the postural control system could still use prediction to take advantage of the statistical properties of the stimulus and internal knowledge about body dynamics to improve performance and to compensate for feedback time delay. However, the Fig. 2 model provided a good description of the experimental data without including any predictive mechanism.
| SUMMARY |
|---|
|
|
|---|
| GRANTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Address for reprint requests and other correspondence: R. J. Peterka, Neurological Sciences Institute, OHSU West Campus, 505 NW 185th Ave., Beaverton, OR 97006, USA (E-mail: peterkar{at}ohsu.edu)
| REFERENCES |
|---|
|
|
|---|
Angelaki DE, Wei M, and Merfeld DM. Vestibular discrimination of gravity and translational acceleration. Ann NY Acad Sci 942: 114127, 2001.[Web of Science][Medline]
Bendat JS and Piersol AG. Random Data: Analysis and Measurement Procedures. New York: John Wiley & Sons, 2000.
Bosco G and Poppele RE. Reference frames for spinal proprioception: kinematics based or kinetics based? J Neurophysiol 83: 29462955, 2000.
Bosco G, Poppele RE, and Eian J. Reference frames for spinal proprioception: limb endpoint based or joint-level based? J Neurophysiol 83: 29312945, 2000.
Cenciarini M. Stimulus Dependent Changes in the Vestibular Contribution to Human Postural Control (Masters thesis). Bologna, Italy: University of Bologna, 2001.
Courjon JH, Precht W, and Sirkin DW. Vestibular nerve and nuclei unit responses and eye movement responses to repetitive galvanic stimulation of the labyrinth of the rat. Exp Brain Res 66: 4148, 1987.[CrossRef][Web of Science][Medline]
Davies WDT. System Identification for Self-Adaptive Control. London: Wiley Interscience, 1970.
Day BL and Bonato C. Modification of the galvanic sway response by visual conditions. In: Multisensory Control of Posture, edited by Mergner T and Hlavacka F. New York: Plenum Press, 1995, p. 169172.
Day BL and Cole J. Vestibular-evoked postural responses in the absence of somatosensory information. Brain 125: 20812088, 2002.
Day BL, Guerraz M, and Cole J. Sensory interactions for human balance control revealed by galvanic vestibular stimulation. In: Sensorimotor Control of Movement and Posture, edited by Gandevia SC, Proske U, and Stuart DG. New York: Kluwer Academic/Plenum Publishers, 2002, p. 129137.
Day BL, Severac Cauquil A, Bartolomei L, Pastor MA, and Lyon IN. Human body-segment tilts induced by galvanic stimulation: a vestibularly driven balance protection mechanism. J Physiol 500: 661672, 1997.
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: 479499, 1993.
Dietz V. Evidence for a load receptor contribution to the control of posture and locomotion. Neurosci Biobehav Rev 22: 495499, 1998.[CrossRef][Web of Science][Medline]
Diffrient N, Tilley AR, and Bardagjy JC. Humanscale 1/2/3. Cambridge, Massachusetts: MIT Press, 1974.
Dijkstra TMH, Schöner G, Giese MA, and Gielen CCAM. Frequency dependence of the action-perception cycle for postural control in a moving visual environment: relative phase dynamics. Biol Cybern 71: 489501, 1994.[Web of Science][Medline]
Duysens J, Clarac F, and Cruse H. Load-regulating mechanisms in gait and posture: comparative aspects. Physiol Rev 80: 83133, 2000.
Fitzpatrick R, Burke D, and Gandevia SC. Task-dependent reflex responses and movement illusions evoked by galvanic vestibular stimulation in standing humans. J Physiol 478: 363372, 1994.
Gardiner EP, Martin JH, and Jessell TM. The bodily senses. In: Principles of Neural Science (4th ed.), edited by Kandel ER, Schwartz JH, and Jessell TM. New York: McGraw-Hill, 2000, p. 430450.
Geyer H, Seyfarth A, and Blickhan R. Positive force feedback in bouncing gaits? Proc R Soc Lond B Biol Sci 270: 21732183, 2003.[Medline]
Giese MA, Dijkstra TMH, Schöner G, and Gielen CCAM. Identification of the nonlinear state-space dynamics of the action-perception cycle for visually induced postural sway. Biol Cybern 74: 427437, 1996.[Web of Science][Medline]
Goldberg JM, Fernandez C, and Smith CE. Responses of vestibular-nerve afferents in the squirrel monkey to externally applied galvanic currents. Brain Res 252: 156160, 1982.[CrossRef][Web of Science][Medline]
Goldberg JM, Smith CE, and Fernandez C. Relation between discharge regularity and responses to externally applied galvanic currents in vestibular nerve afferents of the squirrel monkey. J Neurophysiol 51: 12361256, 1984.
Hajos A and Kirchner W. Körperlageregelung des menschen bei elekrischer reizung des vestibularsystems. In: Sensory Experience, Adaptation, and Perception, edited by Spillman L and Wooten BR. London: Lawrence Erlbaum, 1984, p. 255280.
Hlavacka F and Njiokiktjien C. Postural responses evoked by sinusoidal galvanic stimulation of the labyrinth. Acta Otolaryngol 99: 107112, 1985.[Medline]
Horak FB and Hlavacka F. Somatosensory loss increases vestibular sensitivity. J Neurophysiol 86: 575585, 2001.
Inglis JT, Shupert CL, Hlavacka F, and Horak FB. Effect of galvanic vestibular stimulation on human postural responses during support surface translations. J Neurophysiol 73: 896901, 1995.
Ishida A and Imai S. Responses of the posture-control system to pseudorandom acceleration disturbances. Med Biol Eng Comput 18: 433438, 1980.[CrossRef][Web of Science][Medline]
Ishida A and Miyazaki S. Maximum likelihood identification of a posture control system. IEEE Trans Biomed Eng 34: 15, 1987.[Web of Science][Medline]
Ivanenko YP, Talis VL, and Kazennikov OV. Support stability influences postural responses to muscle vibration in humans. Eur J Neurosci 11: 647654, 1999.[CrossRef][Web of Science][Medline]
Jeka J, Kiemel T, Creath R, Horak F, and Peterka R. Controlling human upright posture: velocity information is more accurate than position or acceleration. J Neurophysiol 92: 23682379, 2004.
Jeka J, Oie K, Schöner G, Dijkstra T, and Henson E. Position and velocity coupling of postural sway to somatosensory drive. J Neurophysiol 79: 16611674, 1998.
Jeka J, Oie KS, and Kiemel T. Multisensory information for human postural control: integrating touch and vision. Exp Brain Res 134: 107125, 2000.[CrossRef][Web of Science][Medline]
Johansson R, Magnusson M, and Akesson M. Identification of human postural dynamics. IEEE Trans Biomed Eng 35: 858869, 1988.[CrossRef][Web of Science][Medline]
Johansson R, Magnusson M, and Fransson PA. Galvanic vestibular stimulation for analysis of postural adaptation and stability. IEEE Trans Biomed Eng 42: 282292, 1995.[CrossRef][Web of Science][Medline]
Kiemel T, Oie KS, and Jeka JJ. Multisensory fusion and the stochastic structure of postural sway. Biol Cybern 87: 262277, 2002.[CrossRef][Web of Science][Medline]
Lund S and Broberg C. Effects of different head positions on postural sway in man induced by a reproducible vestibular error signal. Acta Physiol Scand 117: 307309, 1983.[Web of Science][Medline]
Magnusson M, Enbom H, Johansson R, and Wiklund J. Significance of pressor input from the human feet in lateral postural control. The effect of hypothermia on galvanically induced body sway. Acta Otolaryngol 110: 321327, 1990.[Medline]
Maurer C, Mergner T, Bolha B, and Hlavacka F. Human balance control during cutaneous stimulation of the plantar soles. Neurosci Lett 302: 4548, 2001.[CrossRef][Web of Science][Medline]
Merfeld DM and Young LR. The vestibulo-ocular reflex of the squirrel monkey during eccentric rotation and roll tilt. Exp Brain Res 106: 111122, 1995.[Web of Science][Medline]
Merfeld DM, Zupan L, and Peterka RJ. Humans use internal models to estimate gravity and linear acceleration. Nature 398: 615618, 1999.[CrossRef][Medline]
Mergner T, Hlavacka F, and Schweigart G. Interaction of vestibular and proprioceptive inputs. J Vest Res 3: 4157, 1993.[Medline]
Mergner T, Maurer C, and Peterka RJ. A Multisensory posture control model of human upright stance. Prog Brain Res 142: 189201, 2003.[Web of Science][Medline]
Mergner T, Nardi GL, Becker W, and Deecke L. The role of canal-neck interaction for the perception of horizontal trunk and head rotation. Exp Brain Res 49: 198208, 1983.[Web of Science][Medline]
Mergner T, Siebold C, Schweigart G, and Becker W. Human perception of horizontal trunk and head rotation in space during vestibular and neck stimulation. Exp Brain Res 85: 389404, 1991.[Web of Science][Medline]
Mirbagheri MM, Barbeau H, and Kearney RE. Intrinsic and reflex contributions to human ankle stiffness: variation with activation level and position. Exp Brain Res 135: 423436, 2000.[CrossRef][Web of Science][Medline]
Mittelstaedt H. Origin and processing of postural information. Neurosci Biobehav Rev 22: 473478, 1998.[CrossRef][Web of Science][Medline]
Nashner LM, Black FO, and Wall III C. Adaptation to altered support and visual conditions during stance: patients with vestibular deficits. J Neurosci 2: 536544, 1982.[Abstract]
Nashner LM and Wolfson P. Influence of head position and proprioceptive cues on short latency postural reflexes evoked by galvanic stimulation of the human labyrinth. Brain Res 67: 255268, 1974.[CrossRef][Web of Science][Medline]
Oie KS, Kiemel T, and Jeka JJ. Multisensory fusion: simultaneous re-weighting of vision and touch for the control of human posture. Cogn Brain Res 14: 164176, 2002.[CrossRef][Medline]
Peterka RJ. Sensorimotor integration in human postural control. J Neurophysiol 88: 10971118, 2002.
Peterka RJ. Simplifying the complexities of maintaining balance. IEEE Eng Med Biol Mag 22: 6368, 2003.[CrossRef][Web of Science][Medline]
Peterka RJ and Black FO. Age-related changes in human posture control: sensory organization tests. J Vestib Res 1: 7385, 1990.[Medline]
Peterka RJ and Loughlin PJ. Dynamic regulation of sensorimotor integration in human postural control. J Neurophysiol 91: 410423, 2004.
Prochazka A, Gillard D, and Bennett DJ. Implications of positive feedback in the control of movement. J Neurophysiol 77: 32373251, 1997a.
Prochazka A, Gillard D, and Bennett DJ. Positive force feedback control of muscles. J Neurophysiol 77: 32263236, 1997b.
Schneider E, Glasauer S, and Dieterich M. Comparison of human ocular torsion patterns during natural and galvanic vestibular stimulation. J Neurophysiol 87: 20642073, 2002.
Séverac Cauquil A and Martinez P. Orientation of the body response to galvanic stimulation as a function of the inter-vestibular imbalance. Exp Brain Res 133: 501505, 2000.[CrossRef][Web of Science][Medline]
van Beers RJ, Sittig AC, and Gon JJ. Integration of proprioceptive and visual position-information: an experimentally supported model. J Neurophysiol 81: 13551364, 1999.
van Beers RJ, Wolpert DM, and Haggard P. When feeling is more important than seeing in sensorimotor adaptation. Curr Biol 12: 834837, 2002.[CrossRef][Web of Science][Medline]
van der Kooij H, Jacobs R, Koopman B, and Grootenboer H. A multisensory integration model of human stance control. Biol Cybern 80: 299308, 1999.[CrossRef][Web of Science][Medline]
van der Kooij H, Jacobs R, Koopman B, and van der Helm F. An adaptive model of sensory integration in a dynamic environment applied to human stance control. Biol Cybern 84: 103115, 2001.[CrossRef][Web of Science][Medline]
Wardman DL and Fitzpatrick RC. What does galvanic vestibular stimulation stimulate. In: Sensorimotor Control of Movement and Posture, edited by Gandevia SC, Proske U and Stuart DG. New York: Kluwer Academic/Plenum, 2002, p. 119128.
Welgampola MS and Colebatch JG. Vestibulospinal reflexes: quantitative effects of sensory feedback and postural task. Exp Brain Res 139: 345353, 2001.[CrossRef][Web of Science][Medline]
Winter DA. Biomechanics and Motor Control of Human Movement. New York: John Wiley & Sons, 2004.
Winter DA, Patla AE, Prince F, Ishac M, and Gielo-Perczak K. Stiffness control of balance in quiet standing. J Neurophysiol 80: 12111221, 1998.
Zajac FE. Muscle and tendon: properties, models, scaling, and application to biomechanics and motor control. In: Critical Reviews in Biomedical Engineering, edited by Bourne JR. Boca Raton: CRC Press, 1989, p. 359411.
This article has been cited by other articles:
![]() |
A. D. Goodworth and R. J. Peterka Contribution of Sensorimotor Integration to Spinal Stabilization in Humans J Neurophysiol, July 1, 2009; 102(1): 496 - 512. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Masani, A. H. Vette, N. Kawashima, and M. R. Popovic Neuromusculoskeletal Torque-Generation Process Has a Large Destabilizing Effect on the Control Mechanism of Quiet Standing J Neurophysiol, September 1, 2008; 100(3): 1465 - 1475. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Schmid, A. Nardone, A. M. De Nunzio, M. Schmid, and M. Schieppati Equilibrium during static and dynamic tasks in blind subjects: no evidence of cross-modal plasticity Brain, August 1, 2007; 130(8): 2097 - 2107. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. L. Day and M. Guerraz Feedforward versus feedback modulation of human vestibular-evoked balance responses by visual self-motion information J. Physiol., July 1, 2007; 582(1): 153 - 161. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Lakie and I. D. Loram Manually controlled human balancing using visual, vestibular and proprioceptive senses involves a common, low frequency neural process J. Physiol., November 15, 2006; 577(1): 403 - 416. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |