1Jenks Vestibular Physiology Laboratory, Massachusetts Eye and Ear Infirmary, Department of Otology and Laryngology, Harvard Medical School, Boston, Massachusetts 02114; and 2Neurological Sciences Institute, Oregon Health Sciences University, Portland, Oregon 97209
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Zupan, L. H., R. J. Peterka, and D. M. Merfeld. Neural Processing of Gravito-Inertial Cues in Humans. I. Influence of the Semicircular Canals Following Post-Rotatory Tilt. J. Neurophysiol. 84: 2001-2015, 2000. Sensory systems often provide ambiguous information. Integration of various sensory cues is required for the CNS to resolve sensory ambiguity and elicit appropriate responses. The vestibular system includes two types of sensors: the semicircular canals, which measure head rotation, and the otolith organs, which measure gravito-inertial force (GIF), the sum of gravitational force and inertial force due to linear acceleration. According to Einstein's equivalence principle, gravitational force is indistinguishable from inertial force due to linear acceleration. As a consequence, otolith measurements must be supplemented with other sensory information for the CNS to distinguish tilt from translation. The GIF resolution hypothesis states that the CNS estimates gravity and linear acceleration, so that the difference between estimates of gravity and linear acceleration matches the measured GIF. Both otolith and semicircular canal cues influence this estimation of gravity and linear acceleration. The GIF resolution hypothesis predicts that inaccurate estimates of both gravity and linear acceleration can occur due to central interactions of sensory cues. The existence of specific patterns of vestibuloocular reflexes (VOR) related to these inaccurate estimates can be used to test the GIF resolution hypothesis. To investigate this hypothesis, we measured eye movements during two different protocols. In one experiment, eight subjects were rotated at a constant velocity about an earth-vertical axis and then tilted 90° in darkness to one of eight different evenly spaced final orientations, a so-called "dumping" protocol. Three speeds (200, 100, and 50°/s) and two directions, clockwise (CW) and counterclockwise (CCW), of rotation were tested. In another experiment, four subjects were rotated at a constant velocity (200°/s, CW and CCW) about an earth-horizontal axis and stopped in two different final orientations (nose-up and nose-down), a so-called "barbecue" protocol. The GIF resolution hypothesis predicts that post-rotatory horizontal VOR eye movements for both protocols should include an "induced" VOR component, compensatory to an interaural estimate of linear acceleration, even though no true interaural linear acceleration is present. The GIF resolution hypothesis accurately predicted VOR and induced VOR dependence on rotation direction, rotation speed, and head orientation. Alternative hypotheses stating that frequency segregation may discriminate tilt from translation or that the post-rotatory VOR time constant is dependent on head orientation with respect to the GIF direction did not predict the observed VOR for either experimental protocol.
![]() |
INTRODUCTION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The vestibular system includes two types of
sensors: the semicircular canals and the otolith organs. The
semicircular canals behave as integrating angular accelerometers
measuring head angular velocity (Wilson and Melvill Jones
1979). The otolith organs behave as linear accelerometers
measuring the specific gravito-inertial force (GIF), which is the sum
of gravitational force and inertial force due to linear acceleration.
Using information from these sensors, the human CNS develops
perceptions of spatial orientation and generates VORs that help
stabilize gaze in response to head rotation and translation.
A problem arises when the CNS must distinguish head tilt with respect
to gravity from linear translation (acceleration) of the head. For
example, a compensatory response for tilting the head toward the left
shoulder is a torsional eye movement while acceleration toward the
right induces a compensatory horizontal eye movement. In both cases,
the interaural shear force measured by the otolith organs might be
identical, but the compensatory eye responses are quite different.
Ideally, having devices that independently measure gravity and linear
acceleration could solve the tilt-translation dilemma. However,
Einstein's equivalence principle states that no physical device can
distinguish gravitational force from inertial force due to linear
acceleration. This implies that the CNS is unable to solve the
tilt-translation problem using only information from the otolith organs
or any other combination of linear accelerometers. Indeed, a given GIF
(f) measured by the otolith organs can be
produced by an infinite number of combinations of gravity
(g), and linear acceleration (a). For example,
while in a supine orientation, a linear acceleration of the head to the
left down a slight slope (Fig.
1A) can lead to the exact same
GIF as a slight yaw tilt () to the right (Fig. 1B).
|
In many situations, the eye movements generated in response to head
tilts and translations are appropriately compensatory (Angelaki
et al. 1999; Merfeld and Young 1995
), implying
that the CNS is able to somehow distinguish gravity from linear
acceleration. At least three hypotheses have been proposed to explain
how the CNS makes use of GIF information from the otolith organs to
influence eye movements, peripheral processing, frequency segregation,
and GIF resolution.
The peripheral processing hypothesis states that some of the separation
of the GIF into estimates of gravity and linear acceleration occurs
peripherally with the phasic irregular otolith afferent signal
interpreted primarily as linear acceleration, while the tonic regular
otolith afferent signal might represent gravity (Mayne
1974; Young and Meiry 1967
). The frequency
segregation hypothesis states that the CNS might resolve the GIF
ambiguity using a form of central processing in which low-frequency
cues elicit tilt responses and high-frequency cues elicit translation responses (Mayne 1974
; Paige and Tomko
1991
). This hypothesis is in accordance with low-frequency
characteristics of tilt perception during centrifugation (Clark
and Graybiel 1951
; Glasauer 1992
). It also
matches VOR data during interaural translation, which has both
horizontal translation-related eye movements with high-frequency characteristics and torsional tilt-related eye movements with low-frequency characteristics (Paige and Tomko 1991
).
However, as observed by Angelaki (1998)
, frequency
segregation during interaural translation is affected by the way
translation and tilt sensitivities are expressed (Paige and
Tomko 1991
).
The GIF resolution hypothesis (Merfeld and Young 1995;
Merfeld et al. 1993a
), an explicit refinement of the
"multisensory integration" hypothesis (Guedry 1974
;
Mayne 1974
; Oman 1982
; Young 1984
), states that additional sensory information is required for the CNS to resolve the GIF ambiguity. Specifically, the CNS systematically separates the otolith GIF measurement
(f) into estimates of gravity
(
) and linear acceleration (â) using
multi-sensory convergence, so that the difference between these
estimates approximately matches the measured GIF (f =
â). [The specific gravito-inertial force measured by
the otolith organs (f) is defined as the
sum of gravitational force (g) plus an inertial force per
unit mass (
a) acting on the otolith organs and exactly opposing the direction of linear acceleration (a). We choose to adopt the notation of Young (1984)
for representing
the effect that gravity and linear acceleration have on the otolith
organs. Physical variables are mathematically represented by
three-dimensional vectors, f, g, and
a.] The GIF resolution hypothesis appears
indistinguishable (other than notational differences) from another
recently formulated approach used to explain how the CNS uses
semicircular canal cues to distinguish sinusoidal tilt from sinusoidal
translation (Angelaki et al. 1999
).
When visual cues are absent or cannot be used, the semicircular canal
and otolith cues may provide the CNS with conflicting sensory
information. Consider the sensory situation experienced by a subject
statically positioned in a nose-up orientation (Fig. 2A) following a sustained
counterclockwise rotation (toward the subject's left) in darkness
about an earth-horizontal axis. First, for this experimental protocol,
the otolith organs measure a constant GIF due to gravity alone
(f = g). In the nose-up orientation, this force is aligned with the subject's nasooccipital axis. At the same time, the semicircular canals, because of their dynamics (Wilson and Melvill Jones 1979), have a
post-rotatory response indicating an on-going rotation
(
) in a clockwise (toward the subject's right)
direction even though the subject is actually at rest. Rotational cues
are known to influence the perceived orientation of gravity
(Stockwell and Guedry 1970
), often leading to illusory
tilt (Dichgans et al. 1972
; Merfeld et al.
1999
; von Holst and Grisebach 1951
). This
evidence suggests that the yaw rotational cue following post-rotatory
tilt (
) influences the estimated orientation of
gravity (
) such that the estimate of gravity
(
) rotates in the same direction as if there was an
actual rotation. We hypothesize that the CNS computes estimates of
gravity (
) and linear acceleration (â) such that their difference matches the measured
GIF, which equals gravity in this example (f = g =
â).
Therefore a nonzero estimate of linear acceleration (â =
g) would be generated whenever the
estimate of gravity (
) does not match true gravity (g), as demonstrated in Fig. 2A. The estimate of
interaural linear acceleration
(ây) may induce a horizontal VOR
component (Fig. 2B) that should be similar to a horizontal
VOR response compensatory to an actual interaural linear acceleration;
we refer to this VOR component as an induced VOR. This induced VOR
component should combine more or less linearly with the horizontal
angular VOR component as previously demonstrated in squirrel monkeys
(Sargent and Paige 1991
). Since the tilt direction of
the estimate of gravity (
) depends on rotation
direction, the induced VOR component must depend on rotation direction.
Similarly, since the estimate of interaural linear acceleration
(ây) depends on subject
orientation (
, Fig. 2C), the induced VOR component must
depend on subject orientation.
|
To gain insight into how human subjects distinguish tilt from translation, we designed two experiments that created a conflicting sensory situation between semicircular canal and otolith cues, and we measured reflexive slow phase eye movements. In one experiment, subjects were tilted 90° from an upright orientation after they were rotated for 150 s at a constant velocity about an earth-vertical axis ("dumping" protocol). In another experiment, subjects were rotated for 150 s at a constant velocity about an earth-horizontal axis and then decelerated to a stop ("barbecue" protocol). For both protocols, we focused on the reflexive eye movements following deceleration to a stop, referred to as the post-rotatory responses. Both of these protocols are described in detail in the following section.
Horizontal responses alone have been presented in a previous report for
200°/s dumping protocols (Merfeld et al. 1999). The present report extends this previous study by including VOR responses for 50 and 100°/s dumping protocols to investigate the influence of
the magnitude of the rotational stimulus on the GIF resolution. In
addition, vertical VOR responses are investigated to identify properties of the axis-shift. Moreover responses following
barbecue-spit rotation are presented to evaluate the dynamic influence
of post-rotational tilt on the measured responses.
![]() |
METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Experimental setup and eye movement recording
Informed consent was obtained in accordance with institutional procedures, and subjects were instructed about potential risks, including motion sickness, prior to each testing session.
All protocols were conducted on a two-axis rotation device. An inner gimbal powered by a 160-N-m DC motor (velocity servo control) provided full-circle rotations about the subject's yaw axis. An outer pitch/roll gimbal powered by a 2,300-N-m hydraulic actuator (position servo control) provided rotations about an earth-horizontal axis. Each subject was seated in a kneeling position on a car-race-type seat mounted in the device's inner yaw gimbal. The subject's body was restrained using a 5-point seat-belt system, lateral shoulder supports, a wide waist belt, and knee restraints. Foam pads were added as needed to ensure maximum stability. During testing the subject gripped a pair of handlebars that provided additional stabilization. The subject's head was secured in an adjustable foam-lined head restraint. The subject's interaural axis was aligned with the outer gimbal earth-horizontal axis and the interaural axis midpoint was aligned with the inner gimbal yaw axis.
Binocular eye movements were recorded (Panasonic AG-DS850 SVHS VCR) using small video cameras (Machine Vision Hyper CCD Cameras CV-36SH) mounted on a bite-bar assembly. Infrared light-emitting diodes (LEDs) provided lighting for the video cameras. A mold of each subject's mouth was formed on the bite-bar using a dental-impression compound (3 M Express, 3 M Dental Products, St. Paul, MN). The weight of the camera assembly was supported by elastic bands attached to the yaw gimbal.
An off-line application of a custom image-analysis program provided measurements of the horizontal and vertical coordinates of the pupil center for each field of the video image (59.94 video fields/s). The video system was calibrated by having the subject sequentially direct gaze at 9 horizontal, 9 vertical, and 12 oblique LED target lights. Target viewing evoked horizontal and/or vertical eye movements up to approximately ±20°. Measurements of the three-dimensional (3-D) locations in space of both the eyes and target LEDs were used to determine the horizontal and vertical angular orientations (Fick coordinates) of the eyes with respect to the head for each gaze orientation. The calibration procedure used these known gaze orientations to determine polynomial equations that related the horizontal and vertical pupil center coordinates to the horizontal and vertical eye orientation angles. These polynomial equations were applied to video data obtained during subsequent test runs to determine horizontal and vertical eye orientation. The overall accuracy of the video measurements is provided by the root mean square error between the absolute eye position prediction from the polynomial fit and the presumed actual eye position. The rms error was less than 0.5° for both horizontal and vertical eye orientations. The repeatability of the video measurements is provided by results of redigitizing and reanalyzing a single video image. This procedure produced a standard deviation of 0.05° for both horizontal and vertical eye position measures. All reported eye movement responses are at least one order of magnitude greater than the video-system sensitivity. Horizontal and vertical eye position data were digitally filtered and differentiated to yield horizontal and vertical eye velocity, respectively. Torsional eye position was not calculated. Fast phases were automatically removed using a computer algorithm based on peak acceleration detection, with manual editing by experienced personnel, leaving the slow phase eye velocity (SPV).
Experimental protocols
dumping protocol. Eight healthy subjects age 26-47 (6 males and 2 females) with no history of peripheral or central vestibular disorders volunteered for this study. Clinical testing (including but not limited to rotating chair test battery, Hallpike maneuvers, computerized dynamic posturography, and caloric testing) was performed on seven of the eight subjects and indicated no abnormalities. (The data from the subject not clinically examined were similar to the data from the other subjects.)
The test subjects, with head upright, were accelerated in 2 s to a clockwise (CW) or counter-clockwise (CCW) constant yaw angular velocity about an earth-vertical axis. (A clockwise rotation is defined as a yaw head rotation toward the subject's right shoulder and a counterclockwise rotation as a yaw head rotation toward the subject's left shoulder.) As a convention, the post-rotatory period of a trial is designated by the direction of the preceding rotation. Each subject was tested using three different constant velocities: 50, 100, and 200°/s. The constant velocity rotations were maintained for 150 s and were followed by a 2-s angular deceleration to a stop. To control "viewing" distance, a light was turned on, 10 s prior to the stop. During this 5-s "lights-on" period, the subject was instructed to look straight ahead at a poster fixed to the chair 40 cm in front of the subject's eyes. (Any residual nystagmus was suppressed during this 5-s interval, but a small nystagmus returned by the time the subject was brought to a stop.) Immediately after stopping, the subject was tilted 90° about an earth-horizontal axis through the center of the head (at ear level) in 1.5 s. This passive post-rotatory tilt positioned the subject in one of eight evenly spaced final orientations: nose-up (NU), nose-down (ND), right-down (RD), left-down (LD), and each of the four orientations midway between these principal orientations (NU-RD, NU-LD, ND-RD, ND-LD). The trial order was randomized and secret; the direction of rotation (CW or CCW), and the direction of tilt were counter-balanced across eight subjects for each speed. The same testing order was used for each subject for each speed. The order with which the speeds were presented was counter-balanced across six subjects, and the two remaining subjects were tested using two previously tested speed sequences (200, 100, 50 and 200, 50, 100°/s). The speed was always the same within a test session, and at least two nights separated test sessions. For each speed, there were three testing sessions: one session with all NU and ND trials, one session with all LD and RD trials, and one session with all in-between orientation trials. As controls, each subject was also tested with no tilt following the upright rotation for each velocity (in the 1st 2 sessions). Each subject was also passively tilted 90° to one of eight orientations (in a separate session). To minimize order effects, lights were turned on for 2 min between trials while the subjects were upright and stationary. During data collection, subjects were instructed to keep their eyes open and to look straight ahead but not to focus on any point, real or imagined, and were challenged with mental arithmetic to maintain alertness.barbecue spit protocol. Four healthy male subjects, ages 29-49, with no history of peripheral and central vestibular disorders participated in this portion of the study (2 of the 4 also participated in the dumping protocol). The number of subjects performing this protocol was limited because of severe motion sickness that prevented most subjects from completing the test protocols.
Subjects were initially positioned in a supine orientation and were then accelerated in 2 s to a CW or CCW constant yaw angular velocity of 200°/s about an earth-horizontal axis. The constant velocity rotation was maintained for 150 s before decelerating (2 s) to a stop in either the nose-up (NU) or nose-down (ND) orientations. The remainder of the protocol was identical to the dumping protocol. (Ongoing nystagmus was suppressed during the 5-s light interval preceding the deceleration but returned to normal by the time the subject was brought to a stop.) In addition, two subjects who also participated in the dumping study were tested with the barbecue protocol using all eight stop orientations. The trial order for these two subjects was the same as for their dumping tests.Data analysis
head-fixed reference frame. All vector coordinates (physical variables and internal estimates) were expressed in an orthogonal head-fixed frame of reference with x, y, and z axes corresponding to the subject's nasooccipital, interaural, and rostrocaudal axes respectively. The positive axes were directed nasally (x), toward the left ear (y), and toward the top of the skull (z).
fick angles. To describe eye position and eye angular velocity, we used Fick angles and angular rates defined by the right-hand rule with positive x, y, and z coordinates corresponding to CW torsional, downward vertical, and leftward horizontal movements, respectively.
sinusoidal fit. To characterize the variations of the horizontal induced VOR component, horizontal VOR time constant, and vertical VOR as a function of head orientation, we fit spatial sinusoids to these different data sets for each direction of rotation. For the horizontal-induced VOR component, a sinusoidal variation with head orientation is expected based on the following considerations. The induced VOR component should be proportional to the interaural projection of the estimate of linear acceleration (ây, Fig. 2B). If the estimate of linear acceleration is constant in a space-fixed frame of reference, its interaural projection and the associated compensatory induced VOR component should vary sinusoidally with head orientation in a head-fixed frame of reference (Fig. 6). This sinusoidal variation was quantified by measuring the mean horizontal induced VOR component between 3 and 4 s after the tilt and fitting a sine function to these data (Figs. 6 and 8).
Similarly for the vertical VOR, if the VOR axis-shift component stays constant in a space-fixed frame of reference, its projection on the vertical axis should vary sinusoidally with head orientation in a head-fixed frame of reference. The sinusoidal variation in the vertical VOR with head orientation was quantified by measuring the mean vertical VOR between 10 and 11 s after the tilt (Fig. 10) and fitting a sine function to these data (Fig. 10). We waited 10 s to allow adequate development of the axis-shift, which occurs gradually in humans (Fetter et al. 1996
![]() |
(1) |
time constants.
To determine time constants, a constrained nonlinear optimization
algorithm was used (function const in Matlab 5.2 for
Macintosh, The Mathworks). The slow phase velocity was fit with three
different analytical functions D(t)
![]() |
(2) |
![]() |
(3) |
![]() |
(4) |
statistics. Analysis of variance (ANOVA) was used to determine the statistical significance of velocities and time constants of VOR, angular VOR, and induced VOR components as a function of speed of rotation, direction of rotation, and head final orientation. All statistical analyses were performed with Systat 7.0 (SPSS).
Multivariate ANOVA (MANOVA) methods were applied to the statistical analysis of sinusoidal fit results. Each sinusoidal curve, defined by an amplitude M and a phase ![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
This report focuses primarily on the responses following
deceleration, i.e., slow phase velocity of the post-rotatory nystagmus. However, the per-rotatory responses to yaw rotations about an earth-vertical axis during dumping protocols had average amplitudes and
average decay time constants within the normal range (Table 1) for all speeds (Hess et al.
1985; Honrubia et al. 1984
). The per-rotatory
horizontal VOR time constants decreased with increasing velocity as
previously reported (Baloh et al. 1979
; Paige
1989
). For barbecue spit trials (Table
2), the per-rotatory responses showed the
usual peak amplitudes and sinusoidal modulations (Wall and
Furman 1989
).
|
|
Post-rotatory horizontal eye movements following earth-vertical axis rotations ("dumping")
HORIZONTAL VOR. The horizontal post-rotatory VOR slow-phase velocity was determined in eight subjects for three speeds (200, 100, 50°/s), two directions (CW and CCW) of rotation, and eight final tilt orientations. During the first 20 s following the post-rotatory tilt, the post-rotatory VOR exhibited four very different patterns dependent on rotation direction and head orientations. First, for each of the three speeds tested, the VOR magnitude following a tilt in the NU orientation was greater than the VOR magnitude following a tilt in the ND orientation after identical CCW and CW rotations (Fig. 3, left). We are focusing here on the 20 s following the post-rotatory tilt period that begins 3.5 s after the VOR peaked (see graphs to the right of the vertical dash-dotted line on each plot in Fig. 3). Second, the VOR magnitude was greater for LD trials than for RD trials following identical CCW rotations, but the VOR magnitude was smaller for LD trials than for RD trials following identical CW rotations for each of the three speeds tested (Fig. 3, 2nd column). Third, for NU-RD and ND-LD orientations, little difference was observed following CCW rotation, but the VOR magnitude was greater for NU-RD trials than for ND-LD trials following CW rotation (Fig. 3, 3rd column). Fourth, for NU-LD and ND-RD trials, little difference was observed following CW rotation, but the VOR magnitude was greater for NU-LD trials than for ND-RD trials following CCW rotation (Fig. 3, 4th column). These response patterns were one of our principal findings and were predicted by the GIF resolution hypothesis as discussed later.
|
|
|
ANGULAR- AND ORIENTATION-DEPENDENT COMPONENTS OF THE VOR.
We hypothesized that subject orientations separated by 180° (e.g., NU
and ND) result in an opposite estimate of interaural linear
acceleration (ây),
leading to oppositely directed induced VOR components as well (see
DISCUSSION for more details). At the same time, rotational
cues from the semicircular canals are identical, consistent with
experimental psychophysical reports showing that the duration of
rotation sensation following 90° tilts does not depend on subject
orientation (Benson and Bodin 1966). Therefore if the
VOR slow phase velocities for two head orientations separated by 180°
are noted VOR1 and VOR2, we
calculated, as a simple approximation, an estimate of the angular
VOR (AVOR) component as
![]() |
(5) |
![]() |
(6) |
|
|
|
Post-rotatory horizontal eye movements following an earth-horizontal axis rotation ("barbecue")
HORIZONTAL EYE MOVEMENTS. Four subjects completed barbecue rotations with constant yaw angular velocity of 200°/s about an earth-horizontal axis followed by a stop in either the NU or ND orientations. The barbecue post-rotatory VOR (Fig. 8A) and induced VOR (Fig. 8B) were very similar to the results for 200°/s dumping trials (Figs. 3 and 5). During the 35 s following the tilt, for both CW and CCW rotations, the VOR magnitude following a stop in the NU orientation was greater than following a stop in the ND orientation. For both CW and CCW rotations, the VOR time constant was greater for the NU orientation (12.2 and 12.3 s, respectively) than for the ND orientation (8.7 and 8.9 s, respectively), but due to the small number (n = 4) of subjects tested, the difference was not significant (P = 0.07). These results confirm that the dynamic tilt rotation measured by the semicircular canals in the dumping protocol did not substantially affect or contaminate the horizontal VOR eye movements recorded in the post-rotatory period.
|
Post-rotatory vertical eye movements following earth-vertical axis rotations (dumping)
In addition to the postulated induced VOR component
discussed previously, previous results in humans (Fetter et al.
1996; Harris and Barnes 1987
), rhesus monkeys
(Angelaki and Hess 1994
), and squirrel monkeys
(Merfeld et al. 1993b
) demonstrate that a shift in the
axis of eye rotation is also observed following post-rotatory tilt.
This axis-shift, which is not specifically predicted (nor precluded) by
the GIF resolution hypothesis, indicates the tendency of the VOR
rotation axis to align with gravity, as described by any one of several
spatial orientation hypotheses (Angelaki and Hess 1994
;
Cohen et al. 1999
; Merfeld et al. 1993b
).
While both AVOR and induced VOR components should be horizontal eye
movements, any axis-shift of the VOR rotation axis toward alignment
with gravity should induce either vertical or torsional eye movements, depending on the subject's orientation. The axis-shift of the VOR
rotation axis is not directly related to the induced VOR component, although the induced VOR component would influence the magnitude of the
axis-shift, since the induced VOR component modifies the horizontal VOR.
We evaluated the spatial reorientation of the post-rotatory VOR induced
by head movements in the roll plane by plotting post-rotatory vertical
VOR versus horizontal VOR (Fig. 9) for
both LD or RD orientations, for all speeds, and for both directions of
rotation. The upward vertical components were on average larger than
the downward vertical components. The axis-shift was initially small (<10° for all speeds) at the time when the induced VOR component peaked (about 3 s after the post-rotatory tilt, Fig. 5), and then increased and peaked at an amplitude of roughly 15-30° between 10 and 20 s after the tilt, depending on the rotation angular velocity. A slow buildup in the axis-shift has been previously observed
following post-rotatory tilt (Fetter et al. 1996) and following off-vertical axis rotation in humans (Furman and
Koizuka 1994
; Harris and Barnes 1987
). Figure
10 shows plots of the mean vertical VOR
between 10 and 11 s after completion of the tilt as a function of
head orientation.
|
|
For each rotation direction and speed, the peak amplitude of the
vertical VOR calculated between 10 and 11 s after completion of
the tilt varied sinusoidally as a function of head orientation for each
individual subject (not shown). As described in METHODS, each individual sinusoidal fit can be represented by a pair of real
number [X = M cos , Y = M sin
]
where M is the amplitude and
the phase of the sinusoidal
fit. For two speeds (100 and 200°/s), the group of real pairs
[X1,Y1]
for CW rotations were significantly different (P = 0.02 and P = 0.03, respectively, for Hotelling
T2-test) from the group of real pairs
[X2,Y2]
for the CCW rotations. Since response amplitudes were not significantly
different, we concluded that there was a significant phase difference
between the results from CW and CCW rotations. The mean phase shift for vertical VOR following CW rotations (63 and 66°, respectively) was
significantly different from following CCW rotations (
58 and
64°, respectively).
The mean peak amplitude across subjects of the vertical VOR calculated between 10 and 11 s after completion of the tilt also varied as a function of head orientation (Fig. 10). The CW and CCW responses were well fit by sinusoids shifted with respect to one another. The amplitude, the bias, and the phase shift absolute value of the sinusoidal fits to the mean vertical VOR data decreased with decreasing speed for both CW and CCW rotations.
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
For both dumping and barbecue protocols, the post-rotatory horizontal VOR demonstrated a clear dependence on both head orientation and direction of the preceding rotation, indicating that interactions between otolith and semicircular canal cues influenced the resulting eye movements. In addition, we have shown that the dynamic tilt of the dumping protocol cannot explain the various asymmetries observed in horizontal eye movements since eye responses following yaw rotation about an earth-horizontal axis demonstrated the same asymmetries (Fig. 8). We first compare our results to previous experimental findings and then consider whether any of the existing hypotheses (VOR spatial orientation, frequency segregation, and GIF resolution) that characterize interactions between otolith and semicircular canal cues are consistent with experimentally determined features of dumping- and barbecue-induced post-rotatory VOR.
Horizontal VOR decay following post-rotatory tilts
Previous human studies reported that an eye movement response
attenuation (dumping) occurs when graviceptor cues are inconsistent with rotation cues, e.g., after yaw (Benson and Bodin
1966) and roll (Udo de Haes and Schöne
1970
) rotation about an earth-horizontal axis and after active
(Schrader et al. 1985
) and passive (Benson and
Bodin 1966
) post-rotational tilts following yaw rotation about an earth-vertical axis. Our data confirm these findings, with the decay
time constant of post-rotatory VOR significantly smaller for trials
with post-rotatory tilt than for trials without post-rotatory tilt.
Previous studies also demonstrated asymmetries in VOR responses for
orientations separated by 180° after both passive (Benson and
Bodin 1966
) and active (Schrader et al. 1985
)
post-rotatory tilts. Other paradigms also appear to induce related
orientation-dependent asymmetries. For example, human horizontal
optokinetic nystagmus and afternystagmus induced by yaw rotation of an
optokinetic surround, with subjects and optokinetic surround aligned
with earth-horizontal, have been shown to be greater for NU than for ND
orientations (Wall et al. 1999
). A similar asymmetry has
been observed during caloric stimulation of canal plugged squirrel
monkeys (Minor and Goldberg 1990
; Paige
1985
). These extensive sets of data allow for detailed
comparisons of predictions of different hypotheses regarding the
interactions of otolith and semicircular canal cues on the generation
of VOR eye movements.
VOR spatial orientation hypothesis
The VOR spatial orientation hypothesis has been proposed to
explain orientation dependent asymmetries in the VOR (Cohen et al. 1999). First, this hypothesis states that the angular VOR time constant (10-20 s) is prolonged when compared with primary semicircular canal afferent (5-6 s) by a velocity-storage mechanism (Raphan et al. 1979
). Second, the velocity storage
mechanism is dependent on head orientation with respect to the
direction of the GIF as demonstrated in rhesus and cynomolgus monkeys
(Dai et al. 1991
) and humans (Gizzi et al.
1994
). Third, a "cross-coupling" in the velocity-storage
mechanism occurs from the primary horizontal nystagmus to a secondary
vertical nystagmus after a roll tilt or to a secondary torsional
component after a pitch tilt (Raphan and Cohen 1988
).
Because of this cross-coupling, the rotation axis of reflexive eye
movements demonstrates a tendency to shift toward alignment with the
GIF direction. In several species of monkeys (rhesus, cynomolgus, and
squirrel monkeys), this so-called axis-shift has been observed for
optokinetic nystagmus (OKN) and afternystagmus (OKAN) (Dai et
al. 1991
), for VOR during centrifugation (Merfeld and
Young 1995
; Wearne et al. 1999
), for VOR
following off-vertical axis rotation (Raphan et al.
1992
), and for VOR following a post-rotatory tilt
(Angelaki and Hess 1994
; Merfeld et al.
1993b
). For similar protocols in humans (Fetter et al.
1996
; Gizzi et al. 1994
; Harris and
Barnes 1987
; Merfeld et al. 1998
), the rotation axis of reflexive eye responses also tends to align with the GIF direction but to a much smaller extent.
According to the VOR spatial orientation hypothesis, a directional gain
asymmetry in the secondary nystagmus induced by "cross-coupling" might lead to a time constant asymmetry in the horizontal primary nystagmus (Raphan and Sturm 1991). For example, upward
OKN and OKAN are larger than downward OKN and OKAN in rhesus and
cynomolgus monkeys (Matsuo et al. 1979
). Therefore
horizontal OKN and OKAN responses following a roll tilt in either LD or
RD orientations demonstrate asymmetries that reverse with the direction
of the rotational cues in rhesus and cynomolgus monkeys (Raphan
and Cohen 1988
). Assuming that a similar upward/downward
response asymmetry exists for human VOR responses, the VOR spatial
orientation hypothesis might be consistent with our observed LD/RD
asymmetries in the horizontal VOR and their dependence on rotation direction.
While upward/downward response asymmetries are consistently observed in
rhesus and cynomolgus monkeys, the presence of an upward/downward
response asymmetry in humans is rather controversial and is certainly
not consistent across subjects. An OKN/OKAN upward/downward asymmetry
was reported in humans in some studies (Murasugi and Howard
1989; Wei et al. 1994
) but not in other studies
where asymmetries were reported for only a minority of human subjects
(Baloh et al. 1986
; Stiefel 1962
). No
consistent vertical VOR asymmetry has been reported in humans during
pitch rotation about an earth-vertical axis (Allum et al.
1988
; Baloh et al. 1983
) or about an
earth-horizontal axis (Baloh and Demer 1991
).
It is certainly safe to conclude upward/downward response asymmetries
are inconsistent and at best very weak in humans. Furthermore no
asymmetries in torsional responses have ever been reported in humans
for either optokinetic or vestibular stimulation (Morrow and
Sharpe 1993; Peterka 1992
; Seidman and
Leigh 1989
). Therefore the VOR spatial orientation hypothesis
cannot explain the observed NU/ND asymmetries in the horizontal VOR.
Even if an asymmetry in torsional responses existed, the VOR spatial
orientation hypothesis would predict that the NU/ND asymmetry in the
post-rotatory horizontal VOR should reverse with opposite rotation
direction for both dumping and barbecue protocols. This prediction is
refuted by our dumping and barbecue protocol data that demonstrate a
similar NU/ND asymmetry following both CW and CCW yaw rotations.
Therefore the VOR spatial orientation hypothesis is unable to predict
many characteristics found in the measured human responses. This may be
because this hypothesis has primarily been developed to explain monkey
but not human responses. Furthermore the fact that models of reflexive
eye response spatial orientation (Raphan and Sturm 1991;
Wearne et al. 1999
) do not implement the influence of
rotational cues on the estimation of gravity and linear acceleration may explain some of these shortcomings.
Frequency segregation hypothesis
The frequency segregation hypothesis states that the frequency
content of the GIF measured by the otolith organs can be used to
separate gravity from linear acceleration (Paige and Tomko 1991; Seidman et al. 1998
). Specifically,
low-frequency content might contribute to an estimate of gravity while
high-frequency content might contribute to an estimate of linear
acceleration. Our experimental results are not consistent with this hypothesis.
Rapid tilts yield high-frequency forces measured by the otolith organs. According to the frequency segregation hypothesis, these high-frequency cues should lead to an estimate of linear acceleration in the plane of the post-rotatory tilt and therefore should generate a horizontal linear VOR component when the interaural projection of this estimate of linear acceleration is nonzero. First, there should be a horizontal linear VOR component following a dynamic 90° tilt in either LD or RD orientation with no rotation preceding the tilt. These predictions are inconsistent with the absence of horizontal eye movements for both of these passive tilts. Second, there should be no horizontal linear VOR component when the subject is pitched into either a NU or ND orientation. These predictions are inconsistent with the large asymmetry observed between the NU and ND orientations (Fig. 3). In addition, we have shown that the post-rotatory tilt alone cannot explain the various asymmetries observed in horizontal VOR eye movements since eye responses following yaw rotation about an earth-horizontal axis demonstrated the same asymmetries (Fig. 8). Some additional central processing, probably involving sensory interactions, must be assumed to explain our data.
GIF resolution hypothesis
The GIF resolution hypothesis predicts that the difference between
the estimate of gravity () and the estimate of linear
acceleration (â) matches the measured GIF
(f =
â). Therefore the GIF resolution hypothesis predicts a
nonzero estimate of linear acceleration (â), whenever
the measured GIF, in this case gravity (g), and the estimate
of gravity (
) do not match. After a post-rotatory tilt following a yaw rotation about an earth-horizontal axis, the
post-rotatory semicircular canal cue induces an illusory tilt (von Holst and Grisebach 1951
). Therefore the estimate
of gravity (
) does not match gravity (g)
and the GIF resolution hypothesis predicts a nonzero estimate of linear
acceleration (â, Fig. 2A). If present,
the interaural estimate of linear acceleration should induce a
horizontal VOR. This induced VOR component is similar to the VOR evoked
by actual linear acceleration (Paige and Tomko 1991
;
Schwarz and Miles 1991
; Schwarz et al. 1989
) but is evoked by an estimate of linear acceleration even in the absence of true linear acceleration. The induced VOR component should depend on both rotation direction and head orientation because
the rotation direction influences the direction of illusory tilt and
therefore the direction of the interaural estimate of linear
acceleration and head orientation alters the projection of
the estimated linear acceleration on the interaural axis
(ây, Fig. 2, B and
C). For example, the post-rotational response
(
) for a CCW rotation is CW. For the NU
head orientation, the estimate of gravity (
) tilts
toward the subject's right (Fig. 3, 1st column, top). Since
the estimate of gravity does not match gravity, a nonzero interaural
estimate of linear acceleration (ây) is therefore induced toward
the subject's right. The compensatory eye response to this interaural
estimate of linear acceleration is a horizontal leftward induced VOR
component, which adds to the AVOR component. In contrast, for the ND
orientation following a CCW rotation, both the predicted tilt and
interaural estimate of linear acceleration are toward the subject's
left. This leftward interaural estimate of linear acceleration induces
a rightward horizontal induced VOR component, which now deducts from
the AVOR component. Consistent with these predictions, in the 20 s
following the post-rotatory tilt, the magnitude of the VOR following a
tilt in the NU orientation was greater than the VOR following a tilt in
the ND orientation after identical CCW rotations for each of the three
speeds tested (Fig. 3). When the direction of the rotation preceding
the post-rotatory tilt is CW, the predicted AVOR and induced VOR
components both reversed. Again the induced VOR component adds to the
AVOR component for the NU orientation but deducts for the ND
orientation. Therefore the VOR asymmetry is similar for both rotation
directions (Fig. 3).
For the LD orientation, the interaural estimate of linear acceleration
(ây) is toward the subject's
right for both directions of rotation (Fig. 3, 2nd column),
inducing a leftward horizontal induced VOR component. This induced VOR
component should add to the AVOR component following CCW rotation and
deduct from the angular VOR component following CW rotation. For the RD
orientation, both the interaural estimate of linear acceleration and
the induced VOR component are inverted. Therefore the induced VOR
component should deduct from the AVOR component following CCW rotation
and add following CW rotation. Consistent with previous data
(Schrader et al. 1985), the LD response magnitude was
greater than the RD response magnitude following identical CCW
rotations and less than the RD response magnitude following identical
CW rotations for each of the three speeds tested (Fig. 3).
For orientations midway between the primary tilt orientations, the
observed patterns are also consistent with predictions of the GIF
resolution hypothesis. If we assume that the estimate of linear
acceleration is tilted approximately 45° to the right of vertical
following CCW rotations, the interaural estimate of linear acceleration
is near-zero for NU-RD and ND-LD orientations and maximal for NU-LD and
ND-RD orientations following CCW rotations (Fig. 3). When the rotation
direction is opposite (CW), the estimate of linear acceleration tilts
45° to the left of vertical. Therefore the asymmetry patterns are
inverted. Indeed the maximal interaural estimates of acceleration for
CW trials occur now for NU-RD and ND-LD orientations (Fig. 3), while
near-zero estimates of linear acceleration for CW trials occur for
NU-LD and ND-RD orientations (Fig. 3). The dependence of the horizontal
VOR time courses on head orientation after a post-rotatory tilt (Fig.
3) are predicted by previous VOR models that implement the GIF
resolution hypothesis (Merfeld et al. 1993a;
Zupan 1995
). These models include mechanisms that allow
the combination of otolith and canal information to contribute to the
estimates of both gravity and head angular velocity.
The GIF resolution hypothesis is further supported by sudden increases
observed in the horizontal VOR immediately following the post-rotatory
tilt and by the induced VOR component reversal pattern. First, the VOR
actually increases immediately after the post-rotatory tilt in some
orientations (Fig. 3): this increase is easily noticeable for CW RD, CW
NU-RD, and CCW NU-LD trials at all speeds. This increase is
inconsistent with any type of response attenuation but is consistent
with additive influence of an induced VOR component. Second, for
200°/s rotations, the induced VOR component reverses after about
30 s for all orientations except for LD and RD orientations (Fig.
5, 3rd row). It is well known that semicircular canal
afferents typically show a secondary phase reversal approximately
30 s after deceleration following constant velocity rotation
(Goldberg and Fernandez 1971) and a reversal of VOR eye
movements consistent with this pattern of afferent discharge is
commonly recorded (Correia and Guedry 1966
). As the
semicircular canal cue diminishes and reverses (Fig.
11), the angle between the estimates of
gravity (
) and gravity (g) should diminish
(t2) and then reverse
(t3). Therefore the GIF resolution
hypothesis predicts a reversal of the estimate of linear acceleration
as well as a reversal of the induced VOR component for most
orientations (as illustrated for the NU orientation in Fig.
11B). Such reversals are evident in the horizontal induced VOR component (Fig. 5) following 200°/s yaw rotations. These
reversals are also present following 50 and 100°/s yaw rotations but
are harder to observe.
|
However, in the RD and LD positions, the projection of the estimated linear acceleration on the interaural axis is always in a single direction. For example, in the RD orientation (Fig. 11C), this projection is always toward the subject's left following both CW and CCW rotations. Therefore for the RD and LD orientations, reversal in the estimated tilt does not reverse the interaural estimate of linear acceleration nor does it reverse the induced VOR component even though the VOR for these conditions reverse. Consistent with this prediction, the horizontal induced VOR component does not show a reversal for RD and LD orientations for any of three speeds tested.
According to the GIF resolution hypothesis (Fig. 2), the induced VOR
component is proportional to the interaural estimate of linear
acceleration (ây) and therefore
should vary sinusoidally as a function of head orientation. The induced
VOR component should be maximum and positive when the subject's
interaural axis is aligned with the estimate of linear acceleration and
when this estimate of linear acceleration is toward the right ear.
Therefore the amplitude (M) of the sinusoidal fit to induced
VOR data (Fig. 6, left) should be proportional to the norm
of the estimate of linear acceleration (â). In
addition, the phase shift () of the sinusoidal fits should be
related to the tilt angle that the estimate of linear acceleration
makes with the earth horizontal (
, Fig. 2B).
Specifically, the phase shift (
) should be equal to
for
CCW rotations and equal to 180°
for CW rotations. In
agreement with this prediction, the induced VOR component varied
sinusoidally with subject orientation at each of the three velocities
tested (Fig. 6, left). Moreover, these sinusoidal fits to
the 50, 100, and 200°/s data closely matched the variations of the
projection of the estimate of linear acceleration (Fig. 6,
right) onto the subject's interaural axis
(ây) for angles of the estimate
of linear acceleration with respect to earth-horizontal (
,
Fig. 2B) of 22.5, 32, and 45°, respectively.
In the simple formulation of the GIF resolution hypothesis discussed in
this paper, the estimate of gravity () has a constant
norm [i.e., central processing of canal and graviceptor information
leads to an estimate of gravity (
) such that the
magnitude of
is always a fixed value, presumably 1 g]. Therefore the trajectory of the estimate of linear
acceleration (â) is theoretically an arc of a circle
with a radius equal to the magnitude of
(Figs. 7,
right, and 11). Using simple trigonometry and planar
geometry (Fig. 2B), the theoretical trajectory of the
estimate of linear acceleration (Fig. 7, right) is defined
by the parametric equations
(t) =
â(t)
cos
(t) and
(t) =
â(t)
sin
(t). In
these equations, t is the time after the post-rotatory tilt
and
is the tilt of the estimate of linear acceleration with
earth-horizontal (Fig. 2B). The amplitude of the estimate of
linear acceleration (
â
) is theoretically
proportional to the amplitude (M) of the induced VOR
sinusoidal fit. In addition, the tilt of the estimate of linear acceleration with earth-horizontal for CCW rotations (
, Fig. 2B) and for CW rotations (180°
) is
theoretically equal to the phase shift (
) of the sinusoidal fits to
the induced VOR component. Therefore if the GIF resolution hypothesis
is valid, the induced VOR experimental trajectories (Fig. 7,
left) defined by the parametric equations
X(t) = M(t) cos
(t) and Y(t) = M(t) sin
(t) should closely match
the theoretical trajectories of the estimate of linear acceleration
(Fig. 7, right). Indeed, the experimental trajectories
appear consistent with these predictions.
Influence of speed on the GIF resolution hypothesis
Due to semicircular canal dynamics (Wilson and Melvill
Jones 1979), the semicircular canal post-rotatory cue decreases
with decreasing angular velocity. Since the tilt of the estimate of gravity (
) relative to the earth-vertical axis is
driven by the post-rotatory semicircular canal cue, this illusory tilt should decrease with decreasing speed. Therefore both the amplitude of
the estimate of linear acceleration (â) and the angle (
) of the estimate of linear acceleration with
earth-horizontal (i.e., the phase shift) should decrease with
decreasing speed (Fig. 6, right).
Using trigonometry and planar geometry (Fig. 2B), we can
calculate that the magnitude of the estimate of linear acceleration (â
) when tilted with respect to
earth-horizontal is
â
= 2
g
sin
. This magnitude decreases with decreasing velocity since the tilt angle (
) of the estimate of
acceleration decreases with decreasing velocity. Since the amplitude
(M) of the sinusoidal fits to the induced VOR data are hypothesized to be proportional to the magnitude of the estimate of
linear acceleration (â), this amplitude (M)
should decrease with decreasing velocity. This expectation is
consistent with experimental results showing a decrease in M
with decreasing velocity (Fig. 6).
The decrease in the magnitude of the induced VOR component could simply
be due to the scaling factor of the rotation angular velocity. However,
a simple scaling factor explanation cannot explain the phase shift
variation between sinusoidal fits to CW and CCW data with speed (Fig.
6). In fact, according to the GIF resolution hypothesis, the phase
shift () of the sinusoidal fit to the induced VOR data are related
to the angle (
) of the estimate of linear acceleration
(â) with earth-horizontal. Specifically, the
orientation angle (
, see Fig. 2C) where the induced VOR
component is at its positive maximum is equal, by definition, to the
phase shift (
). This angle also corresponds to the orientation where the subject's interaural axis is aligned with the estimate of linear
acceleration (â), and this estimate of linear
acceleration is directed toward the right ear. As a result,
=
=
for CCW rotations and
= 180°
=
for CW rotations. Therefore with decreasing speed,
the peak of the sinusoidal fits to the induced VOR data should shift
toward 0° for CCW rotations and toward 180° for CW rotations. The
experimental results confirm this theoretical prediction (Fig. 6).
Vertical responses
In addition to horizontal responses, vertical eye movements were also evoked following post-rotatory tilts (Figs. 9 and 10). These vertical VOR eye movements were not compensatory to either the horizontal post-rotatory semicircular canal cues or the post-rotatory tilt itself. Therefore these vertical responses must be due to the interaction of these sensory cues.
As previously reported after post-rotatory tilts following yaw
rotations in rhesus monkeys (Angelaki and Hess 1994) and
squirrel monkeys (Merfeld and Young 1992
; Merfeld
et al. 1993b
), the vertical eye movements observed in the
dumping protocol illustrate the tendency of the VOR axis of rotation to
align with either gravity (g) or an estimate of gravity
(
). Consider a subject tilted to the NU orientation
following CW or CCW yaw rotations. If the VOR displays a tendency to
align its rotation axis with the direction of gravity (g),
the axis-shift component should be exclusively torsional (Fig.
12A,
A-shift1). Alternatively, if the VOR exhibits a
tendency to align its rotation axis with the direction of the estimate
of gravity (
), the corresponding axis-shift component
(Fig. 12A, A-shift2) should include a
secondary vertical axis-shift component in addition to the torsional component. For both CW and CCW directions of rotation, the
post-rotatory secondary vertical VOR component should be positive
(downward) for the NU orientation (Fig. 12A) and negative
(upward) for the ND orientation (Fig. 12B).
|
As shown (Fig. 10), the sinusoidal fits to the vertical VOR data for CW
and CCW rotations are phase shifted with respect to one another, with a
phase shift smaller than 180°. This experimental phase shift is
better matched by the predicted phase shift when the VOR rotation axis
aligns with the estimate of gravity (, Fig.
12D) rather than gravity (g, Fig.
12C). In particular, if the VOR rotation axis aligns with
gravity (g), no vertical eye movement should be present for
NU and ND orientations (Fig. 12C, orientation angle 0 and
180°, respectively). However, experimental results show that a small
vertical VOR is present for both NU and ND orientations, consistent
with an axis-shift toward alignment with the estimate of gravity
(
), as opposed to alignment with gravity.
Conclusion
The GIF resolution hypothesis accurately predicts that post-rotatory VOR eye movements recorded during dumping and barbecue protocols include an "induced" VOR component. The induced VOR component is consistent with the idea of a response compensatory to an internally generated neural estimate of linear acceleration even though no physical linear acceleration is present. The GIF resolution hypothesis accurately predicts the VOR and induced-VOR dependency with head orientation, direction, and speed of rotation. Alternative hypotheses, stating that the post-rotatory VOR time constant depends on head orientation relative to the GIF or that gravity can be separated from linear acceleration by frequency segregation are not consistent with the various characteristics of the measured post-rotatory VOR during dumping and barbecue protocols.
![]() |
ACKNOWLEDGMENTS |
---|
We thank S. Clark-Donovan and J. Roth for assistance throughout the data acquisition and analysis, P. Dimitri for advice on statistical analysis, and Dr. R. Lewis for comments on the manuscript.
This work was supported by a grant from the European Space Agency External Fellowship Program to L. H. Zupan, National Institute on Deafness and Other Communication Disorders/National Aeronautics and Space Administration Grant P60-DC-02072 to R. J. Peterka, and NIDCD Grants R29-DC-03066 and R01-DC-04158 to D. M. Merfeld.
![]() |
FOOTNOTES |
---|
Address for reprint requests: L. H. Zupan, Jenks Vestibular Physiology Laboratory, Suite 421, Massachusetts Eye and Ear Infirmary, 243 Charles St., Boston, MA 02114 (E-mail: lionel_zupan{at}meei.harvard.edu).
1
Indeed, the tilt angle () of the estimate of
linear acceleration with respect to earth-horizontal can be
approximated as
= arctan
(iVORLD-RD/iVORNU-ND),
where iVORLD-RD is the induced VOR component
for the orientation pair LD-RD and iVORNU-ND
is the induced VOR component for the orientation pair NU-ND for a given
rotation direction.
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.
Received 29 December 1999; accepted in final form 16 June 2000.
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|