#### Abstract

We propose a noninvasive method to estimate the time constant. The calculation of this factor permits us to understand the pressure variations of the inner ear and also predict the behavior of the flow resistance of the cochlear aqueduct. A set of mathematical relationships incorporating the intralabyrinthine pressure, the intracranial pressure, and the time constant was applied. The modeling process describes the hydrodynamic effects of the cerebrospinal fluid in the intralabyrinthine fluid space, where the input and output of the created model are, respectively, the sinusoidal variation of the respiration signal and the distortion product of otoacoustic emissions. The obtained results were compared with those obtained by different invasive techniques. A long time constant was detected each time when the intracranial pressure increased; this phenomenon is related to the role of the cochlear aqueduct described elsewhere. The interpretation of this model has revealed the ability of these predictions to provide a greater precision for hydrodynamic variation of the inner ear, consequently the variation of the dynamic process of the cerebrospinal fluid.

#### 1. Introduction

The intracranial pressure (ICP) has been shown to influence the perilymphatic pressure [1–3], because the cranial cerebrospinal fluid (CSF) and inner ear intralabyrinthine fluid communicate through the cochlear aqueduct [4–8]. The ICP changes on a second-to-second basis [9]. Superimposed on the baseline ICP are periodic changes due to cardiovascular activity and respiration [4, 10, 11]. In addition, the ICP varies with the posture change [12, 13]. To achieve a new equilibrium state between the hydrodynamical interactions, a variation of the flow resistance () of the cochlear aqueduct has been linked to the perpetual state of flux [2, 14, 15]. The cochlear aqueduct would be a lowpass filter that should be able to transmit infrasonic waves (i.e., below of 20?Hz) from CSF to the cochlea [6].

Densert et al. have investigated the inner ear pressure by measuring a time constant for pressure release [16–18]. This time constant reflects the pressure change of the inner ear. The experiment was carried out by perforating the tympanic membrane of the cat, applying a square wave pressure at a low frequency of 6.25?mHz, and measuring the time constant by introducing a sound in the perilymphatic fluid. Each pressure applied to the middle ear resulted in a double response in the perilymphatic fluid. The first time constant was referred to as the initial response “”, and it appeared instantaneously upon the application of the stimulus. The second response “” was a direct reaction to the cessation of the pressure stimulus. Later, Wit et al. repeated a similar experiment on the guinea pig and calculated a new time constant called () [19–21]. was mainly determined by the flow resistance of the aqueduct, combined with the compliance of the cochlear windows. The characteristic function of was deduced from a relationship between the intracranial pressure ICP, intralabyrinthine pressure ILP, and time constant . However, because both of these techniques are invasive and use a frequency of 6.25?mHz that belongs to the frequencies band diseases, they cannot be used to measure the time constant and so to characterize the resistance of cochlear aqueduct in normal or pathological conditions in humans.

In this context, the aim of our work was to estimate under hydrodynamic conditions, the variation of the time constant in the human inner ear with a noninvasive technique. In this mathematical model, the estimated time constant reflects the necessary duration of the cochlear aqueduct to reach an equilibrium state between the intracranial and the intralabyrinthine spaces. This idea is based on recent experiments of Büki et al. in 2000 [22]. These authors have shown that distortion product of the otoacoustic emissions (DPOAE) generated around 1?kHz respond to pressure-related stapes impedance changes with change in phase relative to the generator tones, and provide a noninvasive means of assessing intralabyrinthine pressure changes [23]. They also demonstrated from their protocol the absence of any significant confounding middle-ear [24] effect to intracranial pressure ICP. They were described [25] as the relationship between intralabyrinthine and cerebrospinal fluid pressure from the otoacoustic emissions (OAEs) techniques.

Our technique is based on the transmission of infrasonic pressure waves from cerebrospinal to intralabyrinthine fluids through the human cochlear aqueduct [26], by utilizing the recording of the DPOAE stimulus. The estimated values were characterized as a function of the dynamic modulation between intralabyrinthine and intracranial pressures, and then their mean values as a function of the posture body were deduced.

#### 2. Materials and Method

##### 2.1. Subjects

All experiments were performed on young, healthy volunteers with normal-hearing of both genders (4 males, 4 females, age ranging from 22 to 32). Volunteers were instructed not to swallow, to keep reasonably quiet and breathe naturally. They were placed on a tilting table enabling three postures: up-right, supine flat on their back on a horizontal plane, and finally head down with respect to horizontal plane.

##### 2.2. Time Constant Estimation Method

###### 2.2.1. Data Acquisition

To estimate the time constant of the inner ear, DPOAE, and thoracic signals, 250 to 450 points were recorded over 2–3 minutes with a sampling frequency equal to 6.25?Hz, as already previously described [26]. Briefly, DPOAEs are recorded at about 1?kHz after stimulation of two tones of frequency and where ratio is equal to 1.2. The corresponding primary levels and of DPOAE were set at 70?dB SPL in the sealed ear canal. The DPOAE and thoracic signals data (Figure 1) are then filtered to obtain nonlinear smoothing signals and to select the infrasonic waves of the ICP and ILP, corresponding to the respiratory frequency band between 0.17 and 0.5?Hz. These methods have been shown to be effective in removing noise and artifacts from DPOAE without blurring eventual physiological modulations. Thereafter, the normalized ICP and ILP infrasonic waves were used as the input and output signals, respectively, for the model system.

###### 2.2.2. Model System

Based on the work of Gopen et al. [6] and Wit et al. [21] on the inner ear, a simplified mechanical model was created (Figure 2(a)). This model describes the connections and interconnections between the three different compartments of the human body: the respiratory (), the head (), and the auditory systems (). Thereafter, the equivalent electrical parameters to this model (Figure 2(b)): CSF characteristics are represented by the lumpedelement Rcsf, Lcsf, and Cc [26], the cochlear aqueduct connecting the cochlea to these spaces acts as a resistance . The cochlear pressure signal transmitted from the CSF spaces is applied both at the round window and stapes but their large impedance difference implies that only the round window and its compliance can play a significant part.

**(a)**

**(b)**

This model has been created for two objectives. The first one was to calculate the input signal which defines the relationship between the thoracic signal and the cranial fluid pressure. The second was to analytically prove from the output signal, the relationship between the variations of the cochlear response with the movement of the cerebrospinal fluid (CSF).

###### 2.2.3. Input Signal

The input signal is represented by which is the dynamic intracranial pressure ICP related to the rapid fluctuation of breathing. To calculate , we started by modeling the thoracic impedance with the blood variation in the chest. If the imposed changes in pressure and resulting flow are sufficiently small, is represented in terms of a single resistance and single compliance [27] as where is the maximum chest impedance, represents the respiratory variation and is equal to [1/s] with varying between ?Hz and is a constant.

is also defined by the total impedance of the blood and tissue longitudinal impedances given [28] by where and are the blood and tissue longitudinal impedance, respectively, and are expressed in .

At the same time, it is known that where “” is the variation blood volume in the chest, “” expressed in square meters is the cylindrical base area, “” expressed in meters is the chest length, and “” is the resistivity of the blood [28].

Because the respiratory and cardiac movements are transported by the blood to the cerebrospinal fluid [11] in (3), can be modeled by the sinusoidal variation of .

The blood pressure variation [Pa] in the chest, is related to the compliance factor [/Pa] and to the variation of the respiration volume [27] by the following relation: If we consider that, the connection between the respiratory system and the brain has a laminar flow, that is, like a short-circuit (or like a resistance with constant value), then values can be defined as the input signal of the electrical circuit in figure (Figure 2(b)). From (3) and (4) the variation of blood pressure with the impedance was obtained.

This differential equation can be considered as a definite integral with varying between and between []. Therefore, the equation of the input signal is as follows: where is the dynamic intracranial pressure ICP related to the rapid fluctuation of breathing in the cerebrospinal fluid CSF. is the maximal amplitude of .

Equation (6) is a simplified model of ; this equation is valid and may represent the movement of the rapid fluctuation linked to the respiratory system.

*Validation of the Input Signal Equation*

The resulting amplitude demonstrates that the influence of the respiration on CSF flow is decreased in cranial direction [11], and if we calculate the unit of ,

###### 2.2.4. Output Signal

The intralabyrinthine pressure ILP in the inner ear is considered as the output signal (Figures 2(a) and 2(b)), this output signal is represented by . Our objective was to find the relationship between ILP, ICP, and . Based on the physical model described by Feijen et al. [20], the differential equation was resolved, but rather than sending a square signal, a sinusoidal signal of was injected (6), therefore we obtained: where is the maximal amplitude pressure in [Pa] of intralabyrinthine fluid in the cochlea at the time constant in [s] and , and vary with the variation of the respiration frequency [1/s].

We can deduce the relationship between the intracranial and the intralabyrinthine pressure by substituting (6) in (9): As we can see in (11), the relationship between the two factors is strongly related to the frequency of respiration and to the time constant of the inner ear.

By replacing by , a new relationship between ICP and ILP was deduced. The following relationship was obtained: where is , is the maximum intracranial pressure, is the intralabyrinthine pressure (ILP) or the cochlear response at instant “”, and is the dynamic intracranial pressure (ICP) represented by the rapid fluctuation of the breathing at instant “”. Considering as a constant at a given frequency of respiration. The time constant has a role of independent variable, where his positivity is necessary for a valid estimate. The is the constant parameter at “” second. For each selected changes.

#### 3. Results

To estimate the variation of the time constant , the infrasonic waves (between 0.17 and 0.5?Hz) of ICP and ILP from DPOAE and respiration databases were extracted. ICP and ILP are represented by and, respectively (see also paragraph .2). A database for each subject for 3 different postures was recorded separately: up-right (), supine flat on their back on a horizontal plane (), and finally head down () with respect to horizontal plane. Figure 3 gives a representative example of a database for one posture. The recordings of dynamic pressure modulation () and are sinusoidal signals. Each recording can be divided in two different parts depending on the physiological condition. Indeed, the two sinusoidal curves () and of Figure 3 are:

Thereafter, using (12), we estimated for each sequence of each recording obtained for each subject and each posture. A sequence (or wavelength) is composed of two peaks representing the consecutive positive and negative pressures (Figure 4).

Each peak is divided into an up-going edge and a down-going edge. The resulting solution was found to contain two time constants variables ( and ). Because was significantly smaller than , we only represent . The edges varying between 6 and 7 points, we obtained 6 to values. We observed that was giving a uniform curve variation for the sequences contained in the same phase signal part; therefore for all the results presented thereafter, we discarded the perturbation signal part.

##### 3.1. as a Function of the Dynamic Pressure Modulation

was estimated for each sequence of the same phase signal part. Then, was plotted as a function of for the positive and for the negative pressures (Figures 5(a) and 5(b), resp.). followed a uniform and regular pattern for all sequences located without the turbulent transition. For the positive pressure, progressively increased on the up-going edge, until it reached a maximum corresponding to the maximum of the and , and then it progressively decreased on the down-going edge. For the negative pressure, it is the opposite: on the down-going edge, increased until it reached its maximum which corresponds to the minimum of () and , while it decreased on the up-going edge. For a given magnitude, is not the same in the negative or positive edge. For all subjects, does not vary symmetrically in the positive and negative pressures (see, e.g., at 0.4 and -0.4 of Figures 5(a) and 5(b)). Therefore, these results show that to increase the modulation of dynamic pressure leads to an increase in the cochlear aqueduct resistance.

**(a)**

**(b)**

##### 3.2. as a Function of Body Posture

was estimated for each sequence in the same phase signal part. The 12-14 values obtained for the positive or negative pressures of one sequence were averaged to obtain a mean value. Then, the mean values were plotted as a function of the body posture: up-right, supine flat and head-down at with respect to horizontal plane (Figures 6(a) and 6(b)).

**(a)**

**(b)**

The mean value is not always going ascendancy during the transition from up-right to supine flat, or from supine flat to head-down and does not vary linearly. As one began to move down subject's head from up-right position to horizontal plane , there was always a variation of value, the mean can continue to increase or decrease according to each person. With head-down tilting at , mean has continued to increase except for two persons and of Figures 6(a) and 6(b), but its value has always remained with incremental changes being much greater than from head-up tilting. For a given posture, the mean varies from one subject to the other. However, the variation range is similar in the positive and negative pressure. Indeed, on the up-right posture, mean varies from 1.2 to 17.1 in the positive pressure and from 1.6 to 22.8 in the negative pressure; on the horizontal posture, mean varies from 1.3 to 24.0 in the positive pressure and from 1.6 to 34.9 in the negative pressure; on the head-down tilting at , mean varies from 6.0 to 42.4 in the positive pressure and from 7.5 to 56.6 in the negative pressure. We can deduce from these mean values, the maximum variation range was provided in the head-down posture at , being 36.4 for the positive pressure and 49.1 for the negative pressure.

#### 4. Discussions

Intralabyrinthine pressure monitoring in humans is potentially interesting in two situations [29]: hydrocephalus (with abnormal ICP) and Menière's disease (with an alleged pathology of endolymphatic hydrostatics). Using invasive measurements is difficult to perform because of the vulnerability of the human inner ear. The changes in OAE levels in ears with Menière's disease have been proposed as being potentially useful indicators, for clinical monitoring of labyrinth function [29–31]. The DPOAE signals provide an original way to estimate the time constant by noninvasive technique.

Very long time constants were observed by Densert et al. [16] after the blockage of the cochlear aqueduct. In our results, the reaches its maximum value relating mainly to overpressure and under-pressure. The positive peak values of the modulation indicates an over pressure, while the negative peak is the under-pressure. The modulation also provides pressure at the cochlear aqueduct. The monitoring and processing of pressure states ILP and ICP, indicates the flow parameters in a fluid assembly; the information of the cochlear aqueduct [20] is provided by this equation: , where is the fluid flow between intralabyrinthine and cerebral space in [/s]. The phenomenon of long-time constants , detected in our results at peak , was explained from this mathematical relationship by the increasing value of the resistance of the cochlear aqueduct, or decreasing of the fluid flow, where the cochlear aqueduct plays a prominent role and inhibits a large amount of fluid between the two spaces [1, 6, 14].

Knowing that is the inner ear pressure, it presents the global variation pressures including the endolymphatic system, perineural and perivascular spaces. The magnitude on the area without turbulent transition in our databases (Figure 3), was almost equal to the modulation magnitude of phase was proportionally varied with the of the cochlear aqueduct pressures. Therefore the change in the time constant can be explained by the change of the flow resistance [2] during inner ear pressure variation provided by the permeability change of the cochlear aqueduct, caused by a change of structures filling the aqueduct and its entrance in scala tympani. Nevertheless, the regulation of the physiological variations of the inner ear pressures appears to be well balanced even in individuals with poor patency of the cochlear aqueduct, probably due to the close hydrodynamic relationship between the endolymphatic and the cerebrospinal fluid systems [32]. On the other hand, during the turbulent transition indicated on Figure 3, the was inferior and nonproportionally varied with the of the cochlear aqueduct pressures. Therefore the change of the time constant is characterized by fast nonlinear curves and generally by variability of one sequence to another. This can be explained by the role of the cochlear aqueduct [4–7]. The cochlear aqueduct quickly modifies its flow resistance; because it acts as a lowpass filter to ICP changes and attenuates frequency components above its cutoff frequency, where pressure equalization caused by the different factors takes place within seconds.

We can quantify the characteristic of , if we divide the obtained by the window compliance [21]. As we know is defined by the ratio of window pressure and volume for the human (round + oval) window. Ivarsson and Pedersen demonstrated the relationship between the window pressure and this compliance [33]. This relationship demonstrates that, only changes a few percentage points for the range of window pressure variation during an experiment. Consequently, the characteristic behavior of the flow resistance of the cochlear aqueduct can be similar to those obtained from .

Comparisons were made between our curves (Figures 6(a) and 6(b)) with those obtained by Chapman et al. [13], and are similar. The author has determined, by invasive technique, the quantitative relationship between changes in body position and ventricular fluid pressure (intracranial pressure) in normal subjects, using a chronically implanted telemetric pressure sensor. A comparison between this invasive technique and our noninvasive estimation, has been possible because on the normal subject the inner ear fluid pressure can be suddenly changed by the interaction between the intracranial and intralabyrinthine fluids in the normal labyrinth, consequently the inner ear can mirror the variation of the intracranial pressures [15]; adding also, without the normal condition, many other factors can disturb the regulation pressure of the inner ear, such as a change of ear canal [34], middle ear pressure [18] or injection of fluid [19].

Studies in guinea pigs have shown a relation between acute inner ear pressure changes and cochlear function at low-level DPOAE. The inner ear pressure was represented by the variation of the time constant [3]. The authors studied the amplitude and phase of the DPOAE and showed the effect of the DPOAE stimulus and inner ear pressure changes, during and after injection of 0.5?L of artificial perilymph in the inner ear. The two primary frequencies were set at 6?kHz () and 7.5?kHz () with ratio = 1.5, with intensities set, respectively, at 65?dB SPL () and 55?dB SPL (). These high frequencies are used to provide a DC special condition, where the direct current flow is caused to vary very slightly in a special way, with the tiny variations happening at a very high speed, so as to produce a square wave acoustical signal from the DPOAE stimulus when compared to a conventional condition of the input square signal applied by the previous work of Feijen et al. [20]. However, we cannot compare our prediction results to this experiment because we do not have the same frequencies or the same DPOAE stimulus. In our experiment, the DPOAE frequencies have been around 1?kHz and the phase of DPOAE only has been studied. The modeling process of these estimations from our experiment provides a good “apparent” behavior of , because the ICP variations were large (350?mm water or daPa), and OAEs seemed little affected above 2?kHz. Below 1?kHz, their levels tended to decrease by hardly more than 2?dB. The most compicuous change concerned the phase of low-frequency components of OAEs: below 1.8?kHz, it tended to lead the reference phase measured in low ICP conditions, and the size of the phase lead was found to be maximum at around 1?kHz. The maximum phase lead, averaged across ears, turned out to be proportional to ICP increase, according to the linear regression [23]: The main advantage of this estimation is to obtain the time constant of the pressure variation of the inner ear by a noninvasive technique. We believe that, if we take into account this simple predictor provided by the modeling process, it is possible to use the results of infrasonic waves extracted from the DPOAE stimulus around 1?kHz as an aide to predict the hydrodynamic behavior of the cochlear aqueduct below 20?Hz, therefore the relationship of the hydromechanical interactions between the intracranial and intralabyrinthine fluids.

#### 5. Conclusion

The modeling process is the only way to understand the hydrodynamical interactions between the intracranial and intralabyrinthine fluids in the inner ear, because it provides a noninvasive measurement. Using the fact that intralabyrinthine pressure changes induce characteristic phase shifts of DPOEAs around 1?kHz. The results are consistent with previous estimations of the time constant of the inner ear derived from invasive animal experiments. Noninvasive measurements of the time constant at low-frequency pressure waves may turn out to be applicable to monitoring the normal physiology and pathophysiology of the inner ear. This modeling was applied to healthy subjects, changes in the mathematical equations of this model will be necessary, if we want to estimate the time constant for different pathological conditions.

#### Acknowledgments

The authors would like to acknowledge the support of Eric Le Page, Director of OAEricle Laboratory (Australia). They would like to thank the professors of ESITPA for their encouragement for the finalization of this work.