The vestibular receptor of cupula acts an important role in maintaining body balance. However, the cupula buried in the semicircular canals (SCCs) will be destroyed if it is detached from the relevant environment. The mechanical properties of human cupula still remain ambiguous. In this paper, we explored the cupula’s elastic modulus changing with temperature by experiments and numerical simulation of SCCs model. We obtained 3 volunteers’ nystagmus induced by constant angular acceleration when the temperature of volunteers’ SCCs was 36℃ and 37℃ respectively. The slow-phase velocity of 3 volunteers decreased by approximately 3°/s when the temperature of SCCs reduced by 1℃, which corresponded to the reduction of cupula deformation by 0.3–0.8 µm in the numerical model. Furthermore, we investigated the effects of the variation of endolymphatic properties induced by temperature reduction on cupula deformation through numerical simulation. We found that the decrease of cupula deformation was not caused by the change of endolymphatic properties, but probably by the increase of cupula’s elastic modulus. With the temperature reducing by 1℃, the cupula’s elastic modulus may increase by 6%-20%, suggesting that the stiffness of cupula is enhanced. This exploration of temperature characteristic of human cupula promotes the research of alleviating vestibular diseases.
Research Article
Exploring the Mechanical Properties of Human Cupula by Numerical Analysis of Temperature Experiments
https://doi.org/10.21203/rs.3.rs-139506/v1
This work is licensed under a CC BY 4.0 License
Journal Publication
published 15 Apr, 2021
You are reading this latest preprint version
Semicircular canals
Nystagmus
Slow-phase velocity
Numerical model
Cupula deformation
The semicircular canals (SCCs) in the vestibular system detect the angular motion of the head1. 3 SCCs are approximately mutually orthogonal, each of which terminal dilated cavity contains a gelatinous structure called cupula (see Fig. 1). When the head experiences angular motion, the endolymph filled in the SCCs interacts with the cupula due to the inertia. Then, the cupula and the embedded hair cell bundles are deflected, transmitting the neural signals to the brain and inducing the nystagmus to maintain visual stability [the vestibulo-ocular refex (VOR)]2,3. The cupula acts a crucial role in the maintenance of body balance. However, the SCCs are small and complex structures buried in the inner ear. Moreover, the physical properties of cupula will be changed if it is detached from the relevant environment of the SCCs. Thus, it is inappropriate to thoroughly elucidate the material properties and mechanical behaviors of the cupula by dissecting human SCCs. The vestibular diseases involved in the cupula are important factors inducing vertigo, such as motion sickness, benign paroxysmal positional vertigo, and Meniere’s disease, which cause serious disturbance to people's daily life. The research of alleviating vestibular diseases will be promoted effectively if the mechanical properties of the cupula are clarified.
In the present studies, Silver et al.4 investigated the oyster toadfish by anatomy experiment, and found that the cupula of oyster toadfish was a gel containing collagen connective tissue fibers which acted an important role in strengthening the stiffness of cupula. Mchenry and van Netten5 experimentally measured the cupular Young's modulus of zebrafish. These experiments on animals provide a useful basis for researching the mechanical properties of human cupula. However, the anatomic methods used for animal are not suitable for humans. The current technology is not appropriate to directly measure the cupula in human SCCs. Hence, some researchers studied the mechanical behavior of the cupula in SCCs by establishing mathematical or numerical model2,6−11. For instance, Selva et al.6 built a finite-element model of the cupula and estimated Young’s modulus of human cupula. The cupula is considered as a material of collagen gel because it is a gelatinous flexible structure containing collagen connective tissue fibers4,6,12. Collagen is an important factor affecting the stiffness of cupula, and the elastic modulus of collagen gel4,13. With the rise of collagen concentration, the elastic modulus of collagen gel increases13,14. Furthermore, the variation of elastic modulus of different concentration in collagen gel changes diversely with temperature14. The investigation of the cupula’s elastic modulus changing with temperature may identify the concentration of collagen gel with similar temperature characteristics to the cupula, which is also helpful to determine the approximate range of collagen concentration in the cupula. It plays an important role in clarifying the mechanical properties of cupula and further researching the treatment for vestibular diseases.
Based on the principle of VOR, the cupula is deflected and the physiological response of nystagmus is induced, when the SCCs experience angular motion. The slow-phase velocity (SPV) of nystagmus is an external characteristic of cupula deformation, which rises with the cupula deformation increasing15. Additionally, the temperature of the segment of horizontal semicircular canal closest to the temporal bone can be reduced by 1℃ when the external auditory canal is provided irrigation of cold air or water that is 7℃ below nominal body temperature (37℃)16,17. In this study, we explored the elastic modulus of cupula changing with the temperature by the VOR experiment and numerical simulation of SCCs model. Firstly, 3 volunteers’ nystagmus under the condition of SCCs’ temperature at 36℃ and 37℃ respectively was induced by constant angular acceleration in the VOR experiment. We compared the SPV of volunteers’ nystagmus when their SCCs’ temperature was 36℃ and 37℃, respectively. Then, the variation of cupula deformation caused by temperature reduction was quantitatively identified in the numerical model of SCCs. Moreover, we explained that the change of the cupula’s elastic modulus might be the dominant factor affecting cupula deformation by describing intricate fluid–structure interaction in the numerical model of SCCs. Finally, we quantitatively investigated the variation of cupula’s elastic modulus if it was caused by the temperature reduction of 1℃ in SCCs.
Volunteers’ SPV.
Figure 2 shows the trajectory of a volunteer’s nystagmus when the temperature of SCCs was 37℃. The trajectory with the negative slope represents the slow phases of nystagmus. We found that the SPV in the first 3 s increased gradually, and became stable in the last 4 s. This phenomenon was consistent with the report of Bockisch et al.18. Then, we discarded the volunteers’ SPV in the first 3 s and calculated the average of the SPV in the last 4 s (see Table 1). Table 1 shows the average of the SPV of 3 volunteers’ horizontal nystagmus when the temperature of SCCs was 36℃ and 37℃ respectively under the constant angular acceleration of 30°/s2. When the temperature of SCCs was 36℃, the volunteers’ SPV was less than that of the SCCs at 37℃. The results revealed that the volunteers’ SPV induced by constant angular acceleration of 30°/s2 was reduced by approximately 3°/s with the temperature of SCCs reducing by 1℃, which indicated that the cupula deformation decreased. In addition, we found differences in SPV for volunteers when they were simulated by same angular acceleration. The variation of SPV was also different for volunteers when the temperature of SCCs reduced. These phenomena may be caused by individual differences.
|
Nystagmus SPV (mean, unit: °/s) |
||
|---|---|---|
|
Temperature |
36℃ |
37℃ |
|
First volunteer |
28.92 |
32.08 |
|
Second volunteer |
28.27 |
31.67 |
|
Third volunteer |
26.36 |
29.17 |
Responses of SCCs.
As the temperature of SCCs reduced by 1℃, the volunteers’ SPV decreased by approximately 3°/s. According to the relationship between cupula deformation and SPV, the SPV rises slowly with the cupula deformation increasing where approximate 0.1-0.2667 µm of the maximal cupula deformation can induce 1°/s of horizontal nystagmus SPV considering individual differences15. When the SPV decreased by 3°/s, the cupula deformation could be considered as a linear change with the SPV because both the variation of cupula deformation and SPV were very small. Meanwhile, considering that the SPV of volunteers was different due to individual differences, the maximal cupula deformation in the numerical model of SCCs reduced by 0.3–0.8 µm when the SPV decreased by 3°/s. In order to explain the reason for the decrease of cupula deformation, we quantitatively investigated the effects of the changes of endolymphatic properties including density and viscosity on cupula deformation due to the temperature reduction of SCCs. The transcupular pressure, cupula deformation and strain was compared with the same elastic modulus of cupula (5 Pa) when the temperature of SCCs was 36℃ and 37℃ respectively.
When the SCCs rotate with the head, the endolymph is pulled by SCCs due to the action of inertia. The endolymphatic absolute velocity has similar trend with the angular velocity of SCCs. There are vortexes in the utricle and tube flow in the narrow duct by observing the relative flow of the endolymph (see Fig. 3), which is consistent with the study of Boselli et al.7. The cupula deformation is induced by the transcupular endolymph pressure19. 3 cupulae deformation at the temperature of 37℃ is shown in Fig. 4, 5A, 5B, and 5C. The maximal cupula deformation increases with the angular velocity of SCCs during 0–1 s. After 1 s, the cupula deformation becomes stable because the force on cupula induced by endolymphatic pressure is equal to elastic force of cupula. The maximal cupula deformation in HC is greater than that in AC and PC because AC and PC are less stimulated in perceiving angular motion of the head. Besides, sensory cells are distributed on the crista surface of which the strain triggers the activation of sensory hair bundles. The SPV depends on the cupula strain. Then, we explored the HC cupula strain (see Fig. 5D). The result reveals that the strain at the center of crista surface is the largest.
When the SCCs are stimulated by the same angular acceleration at 36℃ and 37℃ respectively, the transcupular pressure, the maximal cupula deformation and the strain in HC are almost identical (see Fig. 6), indicating that the changes of the endolymphatic density and viscosity caused by temperature reduction have no effects on the cupula deformation and strain. The cupula is a structure of collagen gel whose elastic modulus is sensitive to the change of temperature. When the temperature of SCCs reduced by 1℃, the variation of cupula’s elastic modulus might be the dominant factor affecting the cupula deformation. Furthermore, the cupula was considered as a linearly elastic material due to its small deformation8. Thus, the maximal cupula deformation reduced by 0.3–0.8 µm considering the individual differences of the volunteers when the temperature of SCCs decreased by 1℃. The elastic modulus of cupula increased by 0.3-1 Pa, and the percentage increased by 6%-20%, which could induce the corresponding reduction of cupula deformation.
Hair cell bundles are located on the crista surface in the cupula. The cupula strain at the crista surface induces the activation of the hair cell bundles that transmit the neural signals to the brain and trigger the corresponding nystagmus. Although the SPV depends on the cupula strain, the SPV is also dependent on the cupula deformation because the cupula deformation and strain are correlated15. Thus, the SPV can be the external characteristic of cupula deformation. Besides, the maximal cupula strain is at the center of crista surface (see Fig. 5D), which is consistent with the descriptions by Goyens et al.10. There are the most sensitive nerves in the center of crista surface20, increasing the sensitivity of detecting angular motion. When the SCCs experience horizontal rotation, the HC is more stimulated than other SCCs. Then, the HC cupula deformation is the greatest. The hair cell bundles embedded in the HC cupula become active, transmitting the neural signals to the brain and triggering horizontal nystagmus, which is consistent with the results of VOR experiment. The numerical model of 3 SCCs we constructed is reliable.
When the volunteers experienced head rotation with the constant angular acceleration of 30°/s2, each volunteer’s SPV at the SCCs’ temperature of 37℃ was greater than that at 36℃, indicating that the decrease of cupula deformation resulted in the decrease of SPV due to the temperature reduction of SCCs. The elastic modulus of cupula, the density and the viscosity of endolymph changed after the temperature of SCCs dropped by 1℃. However, the variation of endolymphatic density and viscosity caused by the temperature reduction had not the effects on the transcupular pressure that generated cupula deformation. Thus, the decrease of cupula deformation was not induced by the variation of endolymphatic density and viscosity, but probably by the increase of the cupula’s elastic modulus. Then, the stiffness of cupula might be strengthened with the temperature reducing by 1℃, which led to the decrease of cupula deformation. The neural signals generated by the hair cell bundles embedded in cupula were reduced, resulting in the decrease of volunteers’ SPV. When the SPV of volunteers decreased by approximately 3°/s, the cupula deformation in the numerical model of SCCs reduced by approximately 0.3–0.8 µm. Some studies suggested that the SPV of nystagmus is proportional to the cupula deformation2,21. Although the SPV increases nonlinear with the cupula deformation rising in our previous work15, the SPV can be considered to be proportional to the cupula deformation when the variation of SPV is very small. Besides, the real elastic modulus of human cupula remains ambiguous because the present technology can not directly measure it. In order to investigate the change of percentage for the cupula’s elastic modulus, we set the elastic modulus of cupula as 5 Pa in numerical simulation, which was estimated by Selva et al.6. The variation of percentage for the cupula’s elastic modulus increased by 6%-20%, which was not affected by the specific value of elastic modulus. The reason was that the cupula was considered to be a linearly elastic material, and then the elastic modulus of cupula varied linearly with the SPV when the temperature of SCCs changed with the same rotational stimulus. Additionally, the cupula is considered as a soft structure of collagen gel. The elastic modulus of 0 to 3 mg/ml collagen gel is the same order of magnitude as that of cupula6,13,14. When the temperature reduces from 37℃ to 36℃, the elastic modulus of 2 mg/ml collagen gel increases by approximately 6.8%14, which is in the percentage increment range of the cupula’s elastic modulus. The cupula has the similar properties to 2 mg/ml collagen gel, suggesting that the collagen concentration in the cupula may be around 2 mg/ml. In fact, due to individual differences, the properties of cupula for different people may be diverse, including elastic modulus, collagen concentration, and temperature characteristics. The differences in mechanical properties of cupula for people with normal vestibular function may be within reasonable ranges.
The cupula, the vestibular receptor, is associated with vestibular diseases, such as motion sickness, benign paroxysmal positional vertigo and Meniere’s disease. When the temperature of SCCs decreases from 37℃ to 36℃, the elastic modulus of the cupula may increase as well as the stiffness, which improves the ability of cupular anti-deformation. This study provides a novel method to explore the mechanical property of human cupula by VOR experiment and numerical simulation of SCCs model, which promotes the research of relieving vestibular diseases. Moreover, as the cupula is considered to be the material of collagen gel, it is helpful to realize the mechanical properties of the cupula by determining the collagen concentration in the cupula. In this paper, we found that the collagen concentration in the human cupula might be around 2 mg/ml. The elastic modulus and temperature characteristics of 2 mg/ml collagen gel are close to those of the cupula, indicating that they are likely to have more similar mechanical properties. In future studies, the exploration and verification of similar properties for both cupula and collagen gel will act a crucial role in the research of mechanical properties of human cupula, which benefit the study of clinical treatment for vestibular diseases.
Volunteers and equipment in the VOR experiment.
3 volunteers volunteered to participate in the VOR experiment, and gave written informed consent to their participation prior to the experiment. Besides, the statement that all volunteers agreed to publish their identifiable information or images in an online open-access publication was included in the written informed consent. They were informed of the experimental procedures, and were allowed to stop the experiment at any time. All of them had normal vestibular function without history of vestibular or ocular diseases. The experiments were approved by the ethics committee of Dalian university of technology, and were in line with the principles of the Declaration of Helsinki; the registration number was 2020-077.
As shown in Fig. 7A, a volunteer was sitting on the rotatable chair. The experimental equipment including eyepatch, gyroscope, wireless transmission module and battery was used for recording eye movement (see Fig. 7B). A small infrared camera was embedded in the eyepatch to monitor the left eye. Besides, the gyroscope fixed on the right side of the eyepatch could measure instantaneous angular velocity; the sampling frequency of the gyroscope was up to 200 Hz, and we used 50 Hz. There was no relative movement between the eyepatch and the head when the volunteers wearing the eyepatch sensed the angular motion of head. Then, the gyroscope could real-time monitor the angular velocity of volunteers’ head. The wireless transmission module included transmitter and receiver module; the transmitter module connecting to the infrared camera transmitted the video signals of recording eye movement to the receiver module which was connected to the computer; the eye movement videos were stored on the computer. Additionally, the headrest was fixed to the chair to reduce the relative movement of the volunteers' heads during rotation, and the volunteers were secured to the chair with safety belts. Figure 7C shows the working process of the refrigerating device which has the function of producing cold air and transmitting it to the external auditory canal of the volunteers. The refrigerating device consisted of five modules including power source, wind-supply department, refrigerating equipment, heat dissipation device, and air transmission apparatus. Lithium battery was selected as power source, which was portable and supplied appropriate electricity for other modules. The flowing air was generated by the wind-supply department, which was transmitted to the refrigerating equipment through the air transmission apparatus. The refrigerating equipment produced cold air depending on three semiconductor chilling plates in it. There were several air ducts constituting the air transmission apparatus. The terminal of the air ducts were inserted into the soft earplugs which were probed into the ear to transmit the cold air (below 26℃) to the external auditory canal. The earplugs had multiple air holes, which ensured that the air flowed out of the earplug air holes normally without harm for the volunteers. The heat generated by the electronic equipment was dissipated by the cyclic water in the module of the heat dissipation device.
Experimental procedure.
The volunteers sitting on the rotatable chair fastened the safety belts, and wore the eyepatch which was adjusted to record the pupil movement. As shown in Fig. 8A and 8C, the volunteers maintained a head position with the axis of rotation perpendicular to the ground and passing through the middle point P1 of both ears. We marked the auxiliary signs on the chair railings to keep volunteer’s head in the correct position. The chair rotated anticlockwise for 7 s with a constant angular acceleration of 30°/s2, and the initial rotational velocity of 0. The gyroscope recorded the real-time angular velocity of the volunteer’s head when the chair was rotating. All the volunteers firstly participated in the experiment without using the refrigerating device (normal experiment), and then participated in the experiment using the refrigerating device (control experiment). To eliminate the interference between the same volunteers in normal and control experiments, the volunteers rested for at least 30 minutes before the next experiment. All experiments were performed in a dark room to exclude the disturbance of light to eye movement. Besides, the volunteers used refrigerating device for 5 minutes in advance before the control experiment to ensure that the temperature of the SCCs decreased and achieved the balance between heat and cold.
Nystagmus processing.
The recorded videos of eye movement were processed by MATLAB R2017b to track and locate the center of the pupil (see Fig. 7D). We discarded the first and last data of nystagmus slow phases, and removed the phases shorter than 50 ms to reduce the statistical error22. Besides, the first 3 s of the nystagmus was discarded because the SPV generally increases within 3 s after the head senses the change in acceleration18. We also removed the nystagmus data when the volunteers blinked. The SPV was calculated based on the method provided by Wu et al.15.
Numerical simulation.
We reconstructed the 3D geometric model of SCCs in human left ear according to the parameters provided by Ifediba et al.23. The geometric model of SCCs included utricle, three semicircular canals, and three cupulae (see Fig. 8B). The geometric model of SCCs was meshed by Hypermesh (version 12.0). The endolymph domain consisted of 214 k tetrahedral elements and 48 k nodes. The three cupulae domains consisted of 41 k tetrahedral elements and 9 k nodes. Additionally, the results of numerical simulation were almost the same when we refined mesh in endolymph domain with 492 k tetrahedral elements and 103 k nodes, indicating that the selected number of elements for the original model was appropriate.
The cupula is a gelatinous structure10,24,25. Its density is 1000 kg/m3 26,27; the Young’s modulus is 5 Pa6; the Poisson ratio is 0.486,8,28. The endolymph is similar to water25,29−31, with the density of 1000 kg/m3 and viscosity of 0.0085 Pa8,26. These parameters commonly represented the physical properties of the SCCs at the temperature of 37℃. When the temperature of the SCCs dropped to 36℃, the density and viscosity of endolymph increased by 0.036% and 1.83% respectively based on the properties of water32.
We constructed the computer models including endolymph and three cupulae in ANSYS Workbench (version 16.0) based on the method provided by Goyens et al.10. When the SCCs are rotating, the endolymphatic movement in the earthbound reference frame can be expressed by Navier-Stokes equations33:
where ρ is the fluid density, u is the flow velocity vector, P is the static pressure, and µ is the dynamic viscosity. Besides, considering the relative reference frame moving with the walls of SCCs, the walls of endolymph are stationary. Hence, the movement of endolymph in the relative reference frame can be described by the following equations33:
where v is the fluid velocity vector relative to the velocity of the moving reference frame, Ω = (0, 0, ω) is the angular velocity vector of the moving reference frame, and r is the radial coordinate of the fluid element.
We used the Ansys Workbench (version 16.0) to simulate the rotation of SCCs at temperatures of 36℃ and 37℃, respectively. The numerical model was loaded with the anticlockwise rotational acceleration of 30°/s2, and the initial rotational velocity of 0. Furthermore, we observed that the results of cupula response were stable after 2 s in the early calculations. Thus, the computational time was set as 2 s in order to reduce unnecessary calculating time. The time step was set as 0.001 s.
Acknowledgments
This study was funded by the Fundamental Research Funds for the Central Universities [No. DUT20YG112], National Key Research and Development Program of China [No. 2018YFE0197700] and the National Natural Science Foundation of China [Nos. 11572079, 11772087, 31500765]. Furthermore, the authors gratefully thanked M. Jiang, and K. Tao for their help in VOR experiments.
Authors' contributions
The contributions of Xiang Wu were as follows. First, he designed and performed the VOR experiments, and numerical simulation. Second, he analyzed the data, and wrote the manuscript. Shen Yu provided the suggestions of designing the VOR experiment and numerical simulation. Besides, he also revised the manuscript. Shuang Shen contributed to revise and improve the manuscript. Wenlong Liu contributed to modify the design of VOR experiment, and revise manuscript.
Competing interests
The authors declare no competing interests.
- Rabbitt, R. D. et al. Dynamic displacement of normal and detached semicircular canal cupula. J Assoc Res Otolaryngol 10, 497-509 (2009).
- Squires, T. M., Weidman, M. S., Hain, T. C. & Stone, H. A. A mathematical model for top-shelf vertigo: the role of sedimenting otoconia in bppv. J Biomech 37, 1137-1146 (2004).
- Boselli, F., Kleiser, L., Bockisch, C. J., Hegemann, S. C. A. & Obrist, D. Quantitative analysis of benign paroxysmal positional vertigo fatigue under canalithiasis conditions. J Biomech 47, 1853-1860 (2014).
- Silver, R. B., Reeves, A. P., Steinacker, A. & Highstein, S. M. Examination of the cupula and stereocilia of the horizontal semicircular canal in the toadfish Opsanus tau. J Comp Neurol 402, 48-61 (1998).
- Mchenry, M. J. & Van, Netten, S. M. The flexural stiffness of superficial neuromasts in the zebrafish(Danio rerio) lateral line. J Exp Biol 210, 4244-4253 (2007).
- Selva, P., Oman, C. M. & Stone, H. A. Mechanical properties and motion of the cupula of the human semicircular canal. J Vestib Res 19, 95-110 (2009).
- Boselli, F., Obrist, D. & Kleiser, L. Vortical flow in the utricle and the ampulla: a computational study on the fluid dynamics of the vestibular system. Biomech Model Mechanobiol 12, 335-348 (2013).
- Shen, S. et al. Numerical simulation of the role of the utriculo-endolymphatic valve in the rotation-sensing capabilities of semicircular canals. J Biomech 49, 1532-1539 (2016).
- Djukic, T. & Filipovic, N. Numerical modeling of the cupular displacement and motion of otoconia particles in a semicircular canal. Biomech Model Mechanobiol 16, 1-12 (2017).
- Goyens, J., Pourquie, M. J. B. M., Poelma, C. & Westerweel, J. Asymmetric cupula displacement due to endolymph vortex in the human semicircular canal. Biomech Model Mechanobiol 18, 1577-1590 (2019).
- Shen, S. et al. Biomechanical Analysis of Angular Motion in Association with Bilateral Semicircular Canal Function. Biophys J 118, 729-741 (2019).
- Rabbitt, R. D., Damiano, E. R. & Grant, W. J. Biomechanics of the semicircular canals and otolith organs. In The vestibular system (eds. Highstein, S. M. & Fay, R. R.) 153–201 (Springer, 2004).
- Willits, R. K. & Skornia, S. L. Effect of collagen gel stiffness on neurite extension. J Biomater Sci Polym Ed 15, 1521-1531 (2004).
- Yang, Y., Leone, L. M. & Kaufman, L. J. Elastic Moduli of Collagen Gels Can Be Predicted from Two-Dimensional Confocal Microscopy. Biophys J 97, 2051-2060 (2009).
- Wu, X., Yu, S., Liu, W. & Shen. S. Numerical modeling and verifcation by nystagmus slow-phase velocity of the function of semicircular canals. Biomech Model Mechanobiol 19, 2343-2356 (2020).
- Cawthorne, T. & Cobb, W. A. Temperature Changes in the Perilymph Space in Response to Caloric Stimulation in Man. Acta Otolaryngol 44, 580-588 (1954).
- Keck, W. & Thoma, J. Conduction of thermal stimuli in the human temporal bone. Arch Otorhinolaryngol 245, 335-339 (1988).
- Bockisch, C. J., Khojasteh, E., Straumann, D. & Hegemann, S. C. A. Eye position dependency of nystagmus during constant vestibular stimulation. Exp Brain Res 226, 175-182 (2013).
- Yamauchi, A., Rabbitt, R. D., Boyle, R. & Highstein, S. M. Relationship between inner-ear fluid pressure and semicircular canal afferent nerve discharge. J Assoc Res Otolaryngol 3, 26–44 (2002).
- Eatock, R. A. & Songer, J. E. Vestibular hair cells and afferents: two channels for head motion signals. Annu Rev Neurosci 34, 501-534 (2011).
- Obrist, D. & Hegemann, S. C. Fluid-particle dynamics in canalithiasis. J R Soc Interface 5, 1215-1229 (2008).
- Bockisch, C. J., Khojasteh, E., Straumann, D. & Hegemann, S. C. A. Development of Eye Position Dependency of Slow Phase Velocity during Caloric Stimulation. PLoS One 7, e51409 (2012).
- Ifediba, M. A., Rajguru, S. M., Hullar, T. E. & Rabbitt, R. D. The role of 3-canal biomechanics in angular motion transduction by the human vestibular labyrinth. Ann Biomed Eng 35, 1247-1263 (2007).
- Van, Buskirk, W. C., Watts, R. G. & Liu, Y. K. The fluid mechanics of the semicircular canals. J Fluid Mech 78, 87-98 (1976).
- Ekdale, E. G. Form and function of the mammalian inner ear. J Anat 228, 324-337 (2016).
- Rajguru, S. M., Ifediba, M. A. & Rabbitt, R. D. Three-dimensional biomechanical model of benign paroxysmal positional vertigo. Ann Biomed Eng 32, 831-846 (2004).
- Kassemi, M., Deserranno, D. & Oas, J. G. Fluid–structural interactions in the inner ear. Comput Struc 83, 181-189 (2005).
- Wu, C. et al. Dynamic analysis of fuid–structure interaction of endolymph and cupula in the lateral semicircular canal of inner ear. J Hydrodyn Ser B 23, 777–783 (2011).
- Kandel, E. R., Schwartz, J. H. & Jessell, T. M. Principles of neural science, 3rd edn. International edition 1135 (Prentice-Hall International Inc, 1991).
- Latash, M. L. Neurophysiological basis of movement, 2nd edn. Human Kinetics, Champaign (2008).
- Squire, L. R., Berg, D. & Bloom, F. E. Fundamental neuroscience (Elsevier/Academic Press, 2013).
- Kestin, J., Sokolov, M. & Wakeham, W. A. Viscosity of liquid water in the range 8°c to 150°c. J Phys Chem Ref Data 7, 941-948 (1978).
- Batchelor, G. K. An introduction to fuid dynamics (Cambridge University Press, 2007).