 Research
 Open Access
 Published:
Mathematical modelling of the CSF system: effects of microstructures and posture on optic nerve subarachnoid space dynamics
Fluids and Barriers of the CNS volume 19, Article number: 67 (2022)
Abstract
Background
The pressure difference between the eye and brain in upright postures may be affected by compartmentalization of the optic nerve subarachnoid space (ONSAS). Both pressure and deformation will depend on the microstructures of the ONSAS, and most likely also on ocular glymphatic clearance. Studying these factors could yield important knowledge regarding the translaminar pressure difference, which is suspected to play a role in normaltension glaucoma.
Methods
A compartment model coupling the ONSAS with the craniospinal CSF system was used to investigate the effects of microstructures on the pressure transfer through the ONSAS during a posture change from supine to upright body postures. ONSAS distensibility was based on MRI measurements. We also included ocular glymphatic flow to investigate how local pressure gradients alter this flow with changes in posture.
Results
A compartmentalization of the ONSAS occurred in the upright posture, with ONSAS porosity (degree of microstructural content) affecting the ONSAS pressure (varying the supine/baseline porosity from 1.0 to 0.75 yielded pressures between − 5.3 and − 2 mmHg). Restricting the minimum computed porosity (occurring in upright postures) to 0.3 prevented compartmentalization, and the ONSAS pressure could equilibrate with subarachnoid space pressure (− 6.5 mmHg) in \(\le\) 1 h. The ocular glymphatics analysis predicted that substantial intraocularCSF flows could occur without substantial changes in the ONSAS pressure. The flow entering the ONSAS in supine position (both from the intraocular system and from the cranial subarachnoid space) exited the ONSAS through the optic nerve sheath, while in upright postures the flow through the ONSAS was redirected towards the cranial subarachnoid space.
Conclusions
Microstructures affect pressure transmission along the ONSAS, potentially contributing to ONSAS compartmentalization in upright postures. Different pathways for ocular glymphatic flow were predicted for different postures.
Introduction
Severe ocular conditions such as glaucoma and papilledema have been related to abnormal flows and pressures in the cerebrospinal fluid (CSF) within the optic nerve subarachnoid space (ONSAS) [1,2,3,4,5]. The CSF dynamics of the ONSAS are not fully understood, but the CSF communicates between the ONSAS and the intracranial subarachnoid space (SAS), and evidence suggest that CSF is absorbed along the ONSAS through the optic nerve sheath (ONS) to reach the lymphatic system [6,7,8], providing a CSF absorption route and a possible outlet for waste products [8]. Thus, disturbances in the local flow and pressure dynamics of the ONSAS may contribute to disorders related to both the eye and the brain. Due to the small size and sensitive anatomical location of the ONSAS, it is difficult to perform in vivo pressure and flow measurements in this compartment. Therefore, mathematical models provide an attractive tool for investigating the ONSAS dynamics by predicting the flow and pressure numerically, based on available in vivo data.
Most mathematical models and clinical fluid dynamic investigations of the ONSAS mainly reflect horizontal body positions and in special cases headdown tilt [1, 2, 9,10,11,12,13,14], not including upright postures. As craniospinal CSF pressures and volumes change between different postures due to gravity [15,16,17], the pressure and flow in the ONSAS are altered in upright postures compared to horizontal postures. Based on MRI data of the ONSAS, we previously predicted a potential compartmentalization of the ONSAS in the upright posture, hindering flow and pressure transmittance between the eye and brain [18], possibly protecting the eye from low CSF pressures. However, the predictive model used in this previous study did not include the effects of the complex system of microstructures occupying the ONSAS [19]. The ONSAS is traversed by trabeculae, pillars, and septae, which connect the pia around the optic nerve to the arachnoid membrane [20, 21] and may therefore affect flow and pressure in the ONSAS. Kaskar et al. modelled the CSF circulation in the cranial SAS and ONSAS in the supine position and found that the resistance to flow through the ONSAS was heavily dependent on the density of these microstructures [9]. Moreover, they concluded that the ONSAS pressure was highly sensitive to the ONSAS resistance. To fully investigate the ONSAS dynamics, these effects should be included in the predictive modelling.
Furthermore, the findings of a possible glymphatic system for the eye, with CSF entering the ONSAS from the intraocular side and then passing through the optic nerve sheath [8, 22], provide an additional fluid dynamic component that may affect the ONSAS pressure and flows. Investigating the effects of such a pathway in both horizontal and upright postures could contribute valuable insight into ocular glymphatic function and its driving pressure gradients.
Thus, the aim of the current study was to apply mathematical modelling to investigate how ONSAS microstructures affect ONSAS pressure and flow during a change from supine to an upright body posture, thereby further evaluating the previously suggested postural compartmentalization of the ONSAS. A secondary aim was to investigate how these flows and pressures would align with an ocular glymphatic flow pathway. This manuscript will describe the model implementation and showcase results of predicted flows and pressures for the ONSAS, including both temporal changes and steady state levels.
Methods
The postural model for the craniospinal CSF space was first presented by Gehlen et al. [23] building on work by Qvarlander et al. [24] and Magnaes [17]. The ONSAS compartment was subsequently added in Holmlund et al. [18] and in the current study ONSAS microstructures are added through porous media modelling in addition to a CSF pathway from the intraocular system to the ONSAS, reflecting ocular glymphatics [8, 22]. The calculated model variables are listed and described in Table 1. The basal input data (input parameters) to the model can be found in Table 2. After the mathematical model description, a section describing the different tests/simulations is provided.
Mathematical model
Differential equations for the CSF compartments
The model consists of a system of three ordinary differential equations (ODEs) describing the CSF volume change cranially (\(\Delta {V}_{c}\)), spinally (\(\Delta {V}_{s}\)), and in the ONSAS (\(\Delta {V}_{ONSAS}\)), where \(\Delta V\) indicates the deviation in volume from supine equilibrium (\(\Delta V=V{V}_{su{p}_{eq}}\)). The ODE system is solved for t > 0 and reads:
\({Q}_{ou{t}_{c}}\) and \({Q}_{ou{t}_{s}}\) represent cranial and spinal absorption/outflow to venous blood, for example via the arachnoid granulations [25, 26], and \({Q}_{ou{t}_{ONS}}\) the outflow from the ONSAS across the ONS [6, 7]. \({Q}_{cs}\) is the exchange of fluid between the cranial and spinal SAS (cs indicates that the positive flow direction is from the cranial towards the spinal compartment) and \({Q}_{cONSAS}\) is the fluid exchange between the cranial SAS and the ONSAS, and \({Q}_{f}\) corresponds to the CSF formation (applied cranially) and was assumed to be constant [27]. We also included a source of CSF formation across the lamina cribrosa (\(LC\)) membrane at the interface between the ONSAS and the eye, \({Q}_{LC}\) [8, 22]. For an overview of the model and its compartments, see Fig. 1.
Pressure and flow
The CSF outflow rates are determined by the pressure differences and outflow resistances across the absorption pathways [28,29,30], and take the form
The compartmental CSF pressures (\({p}_{comp}\)) are denoted \({p}_{c}\), \({p}_{s}\) and \({p}_{ONSAS}\), corresponding to the cranial CSF pressure (or ICP), spinal CSF pressure and ONSAS pressure, respectively, while the backpressures for CSF absorption (\({p}_{back}\)) are denoted \({p}_{vc}\), \({p}_{vs}\), and \({p}_{orb}\) and represent dural sinus pressure, spinal venous pressure, and intraorbital pressure (the pressure surrounding the ONS) Fig. 1. The hydrostatic reference levels for these pressures are the auditory meatus, the venous hydrostatic indifference point (HIP) and the \(LC\), for the cranial, spinal and ONSAS compartment, respectively. The distribution of cranial and spinal absorption were calculated in the same manner as in Gehlen et al., utilizing the total CSF outflow resistance of the craniospinal system (\({R}_{out}\)) [23]. The outflow across the ONS was defined as:
where \({R}_{ou{t}_{ONS}}\) represents the outflow resistance across the optic nerve sheath (see next subsection for the definition). All CSF absorption routes are modelled as diodes, i.e., no backflow is allowed [31, 32].
For flows in between different CSF compartments, we used the following equations in the model:
where \({R}_{CS}\) is the flow resistance between the cranial and spinal CSF spaces, and \({R}_{ONSAS}\) is the (axial) flow resistance between the cranial SAS and the ONSAS posterior to the eye. The hydrostatic terms on the righthand side of Eqs. (6–7) are included to account for hydrostatic differences between the reference levels for \({p}_{c}\), \({p}_{s}\) and \({p}_{ONSAS}\). \({\rho }_{CSF}\) represents CSF density, \(g\) the gravitational acceleration, and \({h}_{sc}\) and \({h}_{clc}\) are the vertical distances between the different reference points (see Fig. 2).
The inflow across the lamina cribrosa was, in a similar way, set to
where \(IOP\) is intraocular pressure and \({R}_{LC}\) is the flow resistance across the \(LC\).
Specific flow resistances
The resistances related to CSF absorption include \({R}_{out}\) and \({R}_{ou{t}_{ONS}}\). \({R}_{out}\) has previously been assessed by in vivo infusion measurements in healthy [33]. \({R}_{ou{t}_{ONS}}\) was estimated from permeability measurements in porcine ONS [34] (Table 2).
For flows in between compartments, we have \({R}_{CS}\) and \({R}_{ONSAS}\). \({R}_{CS}\) was calculated based on the size of the spinal CSF space at the level of foramen magnum [35] (Table 2). Similarly, the \({R}_{ONSAS}\) was calculated based on the size of the ONSAS. We separated the ONSAS into 14 sections, each being 2 mm in length. These 14 sections together cover 4 regions; bulbar, midorbital 1, midorbital 2, and intracanicular parts of the optic nerve (see Table 3). Subdivision of the ONSAS into several sections is important as the ONSAS decreases in size in upright postures [18, 36], introducing substantial flow resistance in this position. Thus, allowing for variations in size, distensibility and microstructural content along the ONSAS will be crucial for the estimation of this resistance. Assuming that the ONSAS is a porous medium, we may calculate the regional resistances based on Darcy’s law:
where the subscript \(i\) denotes each section of the ONSAS, \(\mu\) is the dynamic CSF viscosity, \({L}_{i}\) is the length of the section, \(\varphi\) is the porosity, and \(A\) is the annular cross section of the ONSAS, i.e., \({A}_{i} =\uppi ({r}_{{ONS}_{i}}^{2}{r}_{{ON}_{i}}^{2})\) where \({r}_{ON}\) is the radius of the optic nerve and \({r}_{ONS}\) is the radius of the optic nerve sheath. \({\kappa }_{i}(\varphi )\) is permeability along the ONSAS calculated as in Kaskar et al. [9, 37, 38] assuming flow perpendicular to the ONSAS microstructures. The baseline (supine) porosity was altered between 0.95, 0.9, 0.85 and 0.75 in separate simulations to determine its effects on ONSAS dynamics. The first two values are similar to that estimated in the SAS in general [39,40,41], and introduce only minor effects on the ONSAS pressure in supine. The last two are included to account for the possibility of a denser ONSAS [20]. The baseline radii and permeabilities used for the ONSAS are presented in Table 3. For comparison, we also calculated the resistance based on Poiseuille flow in an annulus, reflecting the case with a negligible resistance contribution from the trabeculae, pillars and septae:
where \({k}_{i}= \frac{{r}_{{ONS}_{i}}}{{r}_{{ON}_{i}}}\).
\({R}_{LC}\) was calculated using the formula in Eq. (9) and the relationship
where values for the hydraulic conductivity \(\mathrm{K}\) of the \(LC\) could be acquired from Ayyalasomayajula et al. [42] and the thickness/length and area from work by Jonas et al. [43, 44]. There are indications, however, that the flow may mainly pass through specific routes through the \(LC\) [8, 22], suggesting a lower resistance than predicted by Eq. (11). For this reason, we performed separate simulations where \({R}_{LC}\) was altered in a sensitivity analysis.
Relationship between CSF pressure and volume
As is presented by Gehlen et al. [23], the CSF pressure and volume relationship can be separated into a cranial and spinal part written as:
where the constant \({k}_{s}\) is a measure of the spinal contribution to the total compensatory reserve of the CSF system in the supine position [23], since the reference pressures \({p}_{re{f}_{s}}\) and \({p}_{re{f}_{c}}\) are the same in this position, as are the cranial and spinal CSF pressures, i.e. \({p}_{c}\approx {p}_{s}\) in supine [45].
An important assumption of the original Gehlen model was that the reference pressure was dependent on venous pressure:
where \({p}_{v}\) denotes the venous pressure and \({p}_{0}\) is a constant. In the original model, cranial and spinal venous pressures were assumed to be the same in the supine position \({(p}_{vs}={p}_{vc}=CVP)\). However, due to viscous resistance along the veins from the heart to the cranium [46], the cranial and spinal venous pressures are somewhat different even in this position, i.e. \({p}_{vs}\ne {p}_{vc}\) [47]. For Eqs. (12–13) to hold, the constant component of the reference pressure \(({p}_{0})\) was therefore slightly adjusted cranially to accommodate for this difference, such that \({{p}_{re{f}_{s}}=p}_{0s}+{p}_{vs}={p}_{0c}+{p}_{vc}={p}_{re{f}_{c}}\) in the supine position. The viscous resistances were acquired from the literature [46,47,48].
The relationship between ONSAS volume \(\Delta {V}_{ONSAS}\) and pressure \({p}_{ONSAS}\) does not have any established formulation (as opposed to Eqs. (12–13)). From the definition of compliance we can derive an expressionfor this relation:
where we assume that the ONSAS is shaped like an annulus. The last factor on the RHS of Eq. (15) reflects the distensibility of the optic nerve sheath:
The distensibility was acquired from MRIs of healthy volunteers [18]. Since each of the 13 ONSAS segments each had a unique distensibility (\({D}_{i}\)), the \(D\) in Eq. (15) corresponds to the average distensibility of all ONSAS sections (the 14^{th}, or intracanicular section, was assumed as rigid due to its connection to bone). The distensibility was assumed to be constant for each section (for the investigated pressure range) as the change in volume of the ONSAS is relatively small. The effects of this assumption were investigated by limiting the collapsibility of the ONSAS in a separate analysis (see description below).
Postural changes
To model changes in body posture, gravitational effects must be incorporated into the model, which affect both the CSF and the venous system. Within the CSF, hydrostatic effects are already included through Eqs. (6–7). On the venous side, spinal venous pressure in the HIP remains at CVP regardless of body posture, while the cranial venous pressure is affected by hydrostatic effects and the collapse of the jugular veins [47, 49, 50]. The jugular vein collapse is introduced when jugular pressure reaches the surrounding atmospheric zero pressure:
The \({h}_{sijv}\) is simply the vertical distance from the HIP (where we have CVP, or \({p}_{vs}\)) to the jugular veins, and \({p}_{v{i}_{s}}\) is the viscous resistance from the HIP to the jugulars. This creates a zero pressure reference point on the venous side, and cranial venous pressure can then be calculated by the hydrostatic column from the top of the collapse [47, 49, 50]. It is important to note that the ICP decreases with the cranial venous pressure when moving towards the upright posture [53], as reflected in Eqs. (4 and 12–14).
To allow for compartmentalization of the ONSAS during the posture change, the ONS was allowed to change its size (radius) with changes in pressure
where \(D\) is given by Eq. (16) for each individual segment \(i\) (again, except for the intracanicular segment) and baseline corresponds to the supine position. For the porous media models, the porosity will also change with the pressure since the CSF volume decreases with the pressure. By setting a starting porosity in supine (the model baseline), the crosssectional area taken up by the microstructures can be determined as
and as the ONSAS radius is altered with posture (i.e., pressure), the porosity will change as
where \({{A}_{ONSAS}}_{i}\) is the baseline ONSAS crosssectional area. This implementation is equivalent to the area of the microstructures remaining unaltered while any changes in ONS/ONSAS size is due to a reduction in CSF space. Since we measured a unique distensibility for each of the 13 segments of the ONSAS (minus the rigid canicular section), the porosity was allowed to vary along the ONSAS (for pressures above or below the baseline).
Finally, the \(IOP\) also decreases (slightly) when moving towards an upright posture [51]. This decrease was simply implemented as a linear function of the tiltangle (\(\alpha\)). It is unknown if \({p}_{orb}\) changes with posture since in vivo measurements of \({p}_{orb}\) have only been performed in horizontal positions [52]. In the current analysis, we assume a constant \({p}_{orb}\).
Simulations and setups
When running the simulations, the posture change was achieved by linearly increasing the body tiltangle (\(\alpha\)) from 0° to 90° over the course of 10 s. The posture change was set to occur after 5 min in supine equilibrium, and the simulations continued roughly an hour after the posture change. All CSF flows, volumes and pressures were calculated continuously over time. The simulations correspond to a healthy subject in supine equilibrium with a corresponding ICP (\({p}_{{c}_{sup}}\)) of 11.2 mmHg, which is within the normal range measured in vivo [33, 53,54,55]. This corresponds to an ICP of 7.35 mmHg at the level of the LC (\({p}_{{c}_{LC}}\)). The \(IOP\) (at the \(LC\)) started from 18.9 mmHg and was changed to 15.1 mmHg in the upright posture [51]. The simulations were run both with and without the porous media modelling. The baseline porosities tested were 0.95, 0.9, 0.85 and 0.75, where for example 0.75 means that 75% of the ONSAS is occupied by CSF and 25% of microstructures. We also repeated the simulations while imposing limits on the allowed minimum porosity, which is equivalent to limiting the collapsibility (and thus the compartmentalization) of the ONSAS (i.e., a noncollapsible state). The porosity limits tested were 0, 15, 30, 45 and 60%. We also performed a sensitivity analysis regarding the intraocular glymphatic flow by varying the resistance to flow across the \(LC\) (\({R}_{LC})\). The resistance was adjusted to allow for an intraocularONSAS flow from 0.24% (default) to as much as 24% of the aqueous humour production (the fluid that circulates within the eye), which is about 2.5 \(\mathrm{\mu L}\) in total [56]. ONSAS pressure, volume, flow, and (minimum) porosity were the main outcome parameters of interest in the analyses.
Results
Craniospinal CSF pressures and volumes
When moving from supine towards an upright body posture, ICP (\({p}_{c}\)) decreased and the spinal CSF pressure (\({p}_{s}\)) increased. After fully standing/sitting up, both pressures increased slightly, due to the CSF formation being slightly larger than the absorption, eventually reaching equilibrium (where CSF formation and absorption are again equal) after about an hour (Fig. 3). At this point, the \({p}_{c}\) had reached − 3.7 mmHg (from 11.2 mmHg at baseline), corresponding to a postural pressure drop that agrees well with reported measurements [24, 53, 55, 57, 58], while the ICP at the level of the LC (\({p}_{{c}_{LC}}\)) was − 6.5 mmHg (dropping from 7.35 mmHg at baseline) (Fig. 3).
The simulated volumetric changes as a function of time are plotted in Fig. 4. Immediately after the posture change roughly 3.3 mL CSF had moved from the cranial to the spinal compartment, which is also in agreement with reported in vivo measurements [17]. From instant upright to steady state, the spinal volume changed only slightly, whereas the cranial volume increased from 3.3 ml below the supine value to 1.2 mL below the supine value, reaching a new equilibrium.
In supine position, absorption/outflow cranially and spinally was 0.31 mL/min and 0.04 mL/min, respectively. In the upright posture, the corresponding numbers were 0.23 mL/min and 0.12 mL/min. From supine to the upright position, \({Q}_{cs}\) thus increased from 0.04 mL/min to 0.12 mL/min and was always directed from the cranial to the spinal compartment.
ONSAS pressure and flow
For all simulations, the ONSAS pressure (\({p}_{ONSAS}\)) decreased when going from the supine to the upright posture due to the postural decrease in ICP (Fig. 5). We start by presenting the results predicted by Poiseuille flow. Until a certain point, \({p}_{ONSAS}\) decreased in a similar fashion as the \({p}_{{c}_{LC}},\) however, the two pressures were eventually decoupled, and the upright pressures differed (− 5.3 mmHg vs. − 6.5 mmHg) (Fig. 3). This was due to ONSAS compartmentalization, as the ONS almost occluded in the midorbital section of the optic nerve (in a posterior slice of the second midorbital region, see Table 3 for the four regions). The corresponding decrease in ONSAS volume was 0.068 mL (Fig. 4). The results for this configuration reflect those in our previous study [18].
Figure 5 displays \({p}_{ONSAS}\) when varying the baseline porosity of the ONSAS, with \({p}_{{c}_{LC}}\) included as a reference (\({p}_{ONSAS}\) equals this pressure for low ONSAS flow resistances). All the simulations predicted a compartmentalization of the ONSAS and a break in the communication between \({p}_{ONSAS}\) and \({p}_{{c}_{LC}}\) in the upright posture. The ONSAS pressures reached nearequilibrium within 15 min after the posture change had occurred (Fig. 5). The minimum porosity occurred at the point of ONSAS collapse and varied between 0.14 to 0.06 depending on the assumed porosity in supine (0.95 to 0.75). The baseline porosity clearly affected the equilibrium \({p}_{ONSAS}\). For Poiseuille flow the equilibrium \({p}_{ONSAS}\) was − 5.3 mmHg, for a porosity of 0.9 the \({p}_{ONSAS}\) almost equaled the (midbrain) ICP (− 3.9 vs. − 3.7 mmHg). For the lowest (0.75) porosity setting, \({p}_{ONSAS}\) was upheld to − 2 mmHg. The supine \({p}_{ONSAS}\) slightly decreased with decreasing baseline porosity (Fig. 5). The abovementioned findings were also reflected in the ONSAS volume (see Fig. 6) where the total change in volume decreased with decreasing baseline porosity (range 67 to 50 \(\mathrm{\mu L}\) for the porosity values tested).
Imposing a lower limit for the calculated porosity introduced a substantial effect on the dynamics of the different models, as well as a lowering of the equilibrium \({p}_{ONSAS}\) since complete compartmentalization was thus not allowed (see Fig. 7). From lower to higher porosity limits there was a transition from a compartmentalized state (with an elevated equilibrium pressure) to a state where the pressure could not be upheld indefinitely and eventually reached that of \({p}_{{c}_{LC}}\). However, for the lowest porosity limit tested (0.15) the decrease in pressure took a very long time (> 5 h) before reaching a state where the two pressures (\({p}_{ONSAS}\) and \({p}_{{c}_{LC}}\)) were the same. A minimum porosity limit of roughly 0.3 allowed the pressures to equilibrate within an hour upon standing/sitting up. Furthermore, a prohibited compartmentalization allowed for lower pressures to be transmitted along the entire ONSAS, thus, the change in radius close to the \(LC\) and the total volume change of the ONSAS was larger for higher porosity limits. Between the minimum (zero) and maximum (60%) porosity limits tested, this increase was roughly 10 \(\mathrm{\mu L}\) for baseline porosities 0.95 and 0.9, 15 \(\mathrm{\mu L}\) for baseline porosity 0.85 and about 25 \(\mathrm{\mu L}\) for baseline porosity 0.75.
The flow rate through the ONSAS in the supine position was close to 0.002 mL/min (0.6% of \({Q}_{f}\)) for most simulations and was directed from the cranial SAS to the ONSAS compartment. In the upright position, the equilibrium \({Q}_{ONSAS}\) was reversed but nearly stagnant (− 1e−5 mL/min). The outflow via the optic nerve sheet (\({{Q}_{out}}_{ONS}\)) was 0.002 mL/min in supine and zero in upright postures due to the ONSAS pressure dropping below the intraorbital pressure (implemented as a diode).
Intraocular to ONSAS glymphatic flow
The ocular glymphatic flow across the lamina cribrosa (LC) was very low for the calculated \({R}_{LC}\), roughly 6e−6 mL/min in supine and 10e6 mL/min in upright (0.24% and 0.43% of the aqueous humour production, respectively). The difference between the two values was due to the decreasing \({p}_{ONSAS}\) during the change in posture, as the pressure difference \({IOP p}_{ONSAS}\) and the resulting glymphatic flow then increased, which is in line with the literature [8]. For all simulations, the flow through the \(LC\) passed directly across the ONS in the supine position while it instead flowed along the ONSAS to reach the cranial SAS in the upright posture. Thus, the porosity (or resistance) of the ONSAS mainly affects the glymphatic flow in upright postures. The sensitivity analysis for the intraocularONSAS flow revealed that for the lowest tested resistance \({R}_{LC}/100,\) the ONSAS pressure barely changed (from that of the default simulation) in supine, and by less than 1 mmHg in upright. For example, in upright, the abovementioned setting increased \({p}_{ONSAS}\) by 0.3 mmHg for Poiseuille flow (compared to the default simulation), and by 0.73 mmHg for the 0.75 baseline porosity simulation. The corresponding flow rate was 24% of aqueous humour production in supine and 42% in upright for Poiseuille flow while the change was slightly smaller for the 0.75 baseline porosity simulation, 24% to 34%, due to a higher ONSAS pressure in upright for the porous media model. The full sensitivity analysis is displayed in Table 4.
Discussion
To date, the ONSAS pressure and flow alterations between horizontal and upright body postures remain essentially unexplored. Based on existing mathematical models, we developed a new model that allows for prediction of the ONSAS dynamics as a function of body posture including the effects of the microstructures occupying the ONSAS. As the ONSAS is part of the CSF system, a crucial aspect of the model is the interaction between the ONSAS and the craniospinal SAS. This interaction was affected by the content of the ONSAS, possibly allowing for a compartmentalization of the ONSAS in upright postures that introduces significant alterations in steady state ONSAS pressures. The intraocular glymphatic flow increased in upright compared to the supine position, but the flow was redirected from the ONS towards the cranial SAS. Sustaining the glymphatic flow did not require any larger changes in ONSAS pressure even for relatively large intraocularCSF flows in any body posture.
In clinical studies, the ICP is often used as a surrogate for the ONSAS pressure (\({p}_{ONSAS}\)) when estimating potential pressure disturbances between the eye and brain [59,60,61]. Furthermore, the size of the ONSAS posterior to the globe has for some years been investigated as a potential noninvasive indicator of ICP made accessible by ultrasound or MRI [62, 63]. These investigations are based on the assumption of a fully communicating ONSAS and a negligible flow resistance along this pathway. While these assumptions may hold under certain circumstances, e.g., in horizontal body postures (in healthy subjects) or for elevated pressures, they may not hold for all situations. Furthermore, the postural dependency of ICP is often completely overlooked, as the ICP is seldom assessed in upright body postures [51, 58]. Our results open for the possibility that \({p}_{{c}_{LC}}\) and \({p}_{ONSAS}\) are not always the same in upright postures, as even the prediction using Poiseuille flow estimated a small (1 mmHg) difference between the two pressures (Fig. 3). The dependency on ONSAS porosity could mean that differences in porosity, e.g., between different persons or groups of people, could lead to different pressures and flows in the ONSAS (Fig. 5). The model results indicate that a baseline porosity of 0.9 is needed for the \({p}_{ONSAS}\) to mimic midbrain ICP (the common reference for CSF pressure), although the temporal behavior of these pressures differ somewhat (Fig. 5). Thus, if the baseline porosity lies somewhere close to 0.9, midbrain ICP may be a reasonable surrogate for the ONSAS pressure even in upright postures.
One main assumption in our default model [18] is that we extrapolate the distensibility results from 13 to 90 degrees. Imposing a lower limit on the calculated porosity is a way to remedy this, and to investigate the effects of this assumption (Fig. 7). The results show that the limit must be low (roughly 0.1) to uphold a pressure above \({p}_{{c}_{LC}}\) indefinitely, while a limit in between 0.3 and 0.1 may keep the pressure elevated in between one to several hours. Both these options would require a high flow resistance and a low flow in the ONSAS. CT contrast infusions [64] in porcine models has indicated a limited and posture (or gravity) dependent filling of the ONSAS, and could support a limited flow rate in upright postures. If the minimum porosity is higher than 0.3, we can expect relatively quick changes in ONSAS pressure (roughly less than 30 min–1 h to reach \({p}_{{c}_{LC}}\)). Another result of imposing a lower limit on the porosity is that the size of the ONSAS keeps decreasing closer to the bulbar region all the way to the 90degree upright posture, something that is prevented if “complete” compartmentalization occurs. Thus, a larger change in the bulbar ONSAS size may indicate a lack of compartmentalization and a transfer of lower pressures to the back of the eye, or, alternatively, an increased distensibility of the optic nerve sheath.
In the supine position, our simulations predicted a \({p}_{ONSAS}\) equal to the ICP at the level of the \(LC\) (that is, \({p}_{{c}_{LC}}\)) or slightly lower (Fig. 5). While data are sparse, measurements in animals have shown horizontal \({p}_{ONSAS}\) values similar to \({p}_{{c}_{LC}}\) [65] or values a few mmHg lower [66], and measurements in cadavers yielded a \({p}_{ONSAS}\) in the range of 0–6 mmHg [67]. While the latter may not represent the physiological situation, the animal studies support that nonfluid content within the ONSAS may contribute to a nonnegligible resistance already in horizontal positions. Thus, a decreased ONSAS porosity leads to a decreased ONSAS pressure in supine position, but also to an increased ONSAS pressure in upright postures, diminishing the postural effect on the ONSAS pressure in two different ways (Fig. 5).
Our model can be compared to the work by Kaskar et al. who modeled the CSF system in the supine position. Based on the pressure measurements in the cadavers [67], their model predicted an ONSAS resistance of 200–241 mmHg/(mL/min) and an ONS absorption of 5–10% of the total CSF outflow [9]. Their outflow rates are high compared to ours (roughly 5–10% vs. 0.5–0.6%), which may indicate that our \({R}_{ou{t}_{ONS}}\) was overestimated. Decreasing our \({R}_{ou{t}_{ONS}}\) by a factor of 10 would yield a lowered supine \({p}_{ONSAS}\) since the flow rate would then increase, but only for a porosity of 0.75 would the change be larger than 1 mmHg (1.2 mmHg for 0.75, 0.36 mmHg for 0.85, and < 0.15 mmHg for the rest). Since \({p}_{ONSAS}\) quickly dropped below the intraorbital pressure during the posture change, stopping any ONS absorption, any error in \({R}_{ou{t}_{ONS}}\) will have a limited effect on the dynamic changes during the alteration in posture (Fig. 5). A limitation of the Kaskar model is that the ONS absorption is set as a constant and is not pressuredriven, i.e., it is independent of factors such as the intraorbital pressure and does not explicitly include any resistance to absorption across the ONS. Because \({p}_{ONSAS}\) is affected by both \({R}_{ONSAS}\) and \({R}_{ou{t}_{ONS}}\), our model does allow for analyses of these contributions separately (in addition to their postural dependency).
Clinical applications
The ONSAS dynamics are believed to relate to ocular disorders through the translamina cribrosa pressure difference, i.e., the difference between intraocular pressure and ONSAS pressure. A relevant example is normal tension glaucoma where a low ICP has been suspected as a contributing factor to the development of the disorder [5, 59, 68, 69]. Understanding the link between ICP and the pressure within the ONSAS is crucial to accurately interpret findings related to this pressure difference, especially since the transLC pressure difference is often just calculated as IOPICP, which likely does not hold for all postures (Figs. 3, 5). For example, a lack of compartmentalization could allow for lower pressures to reach the posterior eye in upright postures, thus increasing the pressure difference despite normal IOP and ICP [70]. The ONSAS dynamics may also be disturbed in other ways. For example, there is evidence of a hindered ONSAS flow both in patients with normaltension glaucoma [1, 71] and those with papilledema [2] in horizontal positions, which could mean that the ONSAS resistance is instead increased in these patients.
Posture, and thus gravitational effects, is also of importance for understanding the spaceflight associated neuroocular syndrome, a syndrome where astronauts lose their visual acuity after extended visits in microgravity [72]. An abnormal ICP (and by extension an abnormal ONSAS pressure) is one suggested hypothesis [51, 73, 74] and compartmentalization of the ONSAS another [72, 73]. The pressure behaviour as studied with our model may contribute valuable reference data for comparison to the microgravity state. While ICP measurements during longduration spaceflight are nonexistent, measurements in acute microgravity suggest a decrease in both ICP and CVP by roughly 3–4 mmHg [74]. By removing gravitational effects and lowering the CVP accordingly (3.5 mmHg), our model predicts an ONSAS pressure of 7.7 mmHg in microgravity. This is slightly larger than the supine baseline value of 7.35 mmHg. This happens because the hydrostatic difference between the LC and the auditory meatus disappears in microgravity, and this hydrostatic effect is slightly larger than the microgravityinduced change in CVP/ICP. Furthermore, microgravity removes the normal lowering of ICP that occurs in upright postures on earth, resulting in an, on average, higher ICP in microgravity (assuming 16 h per day are spent upright). The difference in this average will depend on the porosity of the microstructures, where a denser (low porosity) ONSAS will lead to a smaller postural effect on the ONSAS pressure on earth.
The extended model could potentially be used for understanding the clearance of eye metabolites, through the ocular glymphatic system [75] of the ONSAS pathway [8]. Assessing local pressure gradients that drive the fluid flow will be essential for increasing our understanding of this clearance system. The sensitivity analysis of the LC flow resistance (Table 4) revealed that flow from the intraocular system could be quite substantial (even up to 24%) without requiring any larger (< 1 mmHg) changes in the ONSAS pressure, despite the high resistance pathway of the ONSAS in the upright posture. As shown by Wang et al. [8], the flow through the posterior of the eye to the ONSAS depends on the \(IOP{p}_{ONSAS}\) pressure difference (by design in our model). This suggests that flow from the intraocular system increases in upright postures. This is in contrast to glymphatic function in general, which has been shown to be increased during sleep (when we are horizontal) [76]. Our model does predict that clearance across the ONS mainly occurs in supine (horizontal) positions however, while in upright positions the flow is directed to the rest of the SAS, suggesting different pathways for the two postures. Thus, intraocular fluid is absorbed in supine while it contributes to CSF production in upright. If absorption across the ONS still occurs in upright postures, it would likely indicate that the counter pressure to this absorption, i.e., the intraorbital pressure, must change in a similar fashion as \({p}_{ONSAS}\). However, we consider this unlikely as the intraorbital pressure can be suspected to remain fairly constant. An alternative is that some other pressure constitutes the counter pressure to this absorption.
Finally, additional CSF compartments can subsequently be added to the current model, moving towards a more complete description of the CSF dynamics and glymphatic functions. As CSF dynamics do vary with body position, including the effects of posture in such models are likely essential. Additional outflow/absorption routes related to the glymphatic system could be included. However, this should not affect the dynamics investigated in this study as absorption in our model is based on the total outflow resistance \({R}_{out}\).
Limitations
We extrapolated the distensibility values from 13 degrees headup tilt (maximum achievable in our scanner) to 90 degrees headup tilt. Imaging in sitting MRIs could contribute important data here, however, in this study, this limitation was instead addressed by the porosity limitation analysis. We did not include any pulsations in this study, including general CSF pulsations as well as local arterial pulsations in the optic nerve and ONS, which may contribute to the mixing and distribution of CSF between the ONSAS and cranial SAS. The relative contribution from the spinal and cranial compartments to CSF compliance has been debated, where some estimations indicate a larger spinal contribution compared to the cranial contribution [17, 77]. Changing this parameter in this way will lower ICP slightly, exacerbating the effect of the compartmentalization. However, this will increase the postural shift in volume to values much higher than that measured in vivo, as opposed to the current setup, suggesting that our choice recreates the in vivo situation satisfactorily. Only CSF absorption to venous blood was included for the craniospinal CSF compartments despite lymphatic outflow existing for these compartments as well [78, 79]. Since accurate postural craniospinal CSF pressures and volume changes were achieved with the current model, they were deemed sufficient for the specific application of the current study (where we focus on the ONSAS), but these other flow pathways need to be added to study for example general glymphatic flow throughout the brain. There is a lack of studies of ocular glymphatics in the upright posture to validate our glymphatics flow predictions. While our model can predict the outcome of different scenarios, more work in this area is needed to validate these predictions and for further improving the model.
Conclusions
The current study presents a model for analyzing the ONSAS fluid dynamics and its interactions with the craniospinal CSF spaces with varying body posture. The porosity related to microstructures within the ONSAS was predicted to affect the ONSAS dynamics during the change in posture from supine to upright as well as the equilibrium ONSAS pressure in upright postures. The results thus suggest that measured ICP may not be an appropriate surrogate for the ONSAS pressure in upright postures. While the current study focused on the ONSAS and the analysis of ocular glymphatics, the model can be expanded to include other CSF compartments, offering a tool for gaining a more complete description of the CSF system and its posture dependency.
Availability of data and materials
No original datasets were used in the present study. All parameters used in the model are described in the methods section and relevant figures.
Abbreviations
 CSF:

Cerebrospinal fluid
 HDT:

Headdown tilt
 HUT:

Headup tilt
 HIP:

Hydrostatic indifference point
 ICP:

Intracranial pressure
 IOP:

Intraocular pressure
 LC:

Lamina cribrosa
 MRI:

Magnetic resonance imaging
 ODE:

Ordinary differential equation
 ON:

Optic nerve
 ONS:

Optic nerve sheath
 ONSAS:

Optic nerve subarachnoid space
References
Killer HE, Miller NR, Flammer J, Meyer P, Weinreb RN, Remonda L, et al. Cerebrospinal fluid exchange in the optic nerve in normaltension glaucoma. Br J Ophthalmol. 2012;96:544–8.
Killer HE, Jaggi GP, Miller NR, Huber AR, Landolt H, Mironov A, et al. Cerebrospinal fluid dynamics between the basal cisterns and the subarachnoid space of the optic nerve in patients with papilloedema. Br J Ophthalmol. 2011;95:822–7.
Berdahl JP, Yu DY, Morgan WH. The translaminar pressure gradient in sustained zero gravity, idiopathic intracranial hypertension, and glaucoma. Med Hypotheses. 2012;79:719–24.
Wall M. Idiopathic intracranial hypertension. Neurol Clin. 2010;28:593–617.
Ren R, Jonas JB, Tian G, Zhen Y, Ma K, Li S, et al. Cerebrospinal fluid pressure in glaucoma: a prospective study. Ophthalmology. 2010;117:259–66.
Killer HE, Laeng HR, Groscurth P. Lymphatic capillaries in the meninges of the human optic nerve. J Neuroophthalmol. 1999;19:222–8.
Lüdemann W, von Rautenfeld DB, Samii M, Brinker T. Ultrastructure of the cerebrospinal fluid outflow along the optic nerve into the lymphatic system. Child’s Nerv Syst. 2005;21:96–103.
Wang X, Lou N, Eberhardt A, Yang Y, Kusk P, Xu Q, et al. An ocular glymphatic clearance system removes βamyloid from the rodent eye. Sci Transl Med. 2020. https://doi.org/10.1126/scitranslmed.aaw3210.
Kaskar OG, Fleischman D, Lee YZ, Thorp BD, Kuznetsov AV, Grace L. Identifying the critical factors governing translaminar pressure differential through a compartmental model. Investig Opthalmol Vis Sci. 2019;60:3204.
Pircher A, Montali M, Berberat J, Remonda L, Killer HE. The optic canal: a bottleneck for cerebrospinal fluid dynamics in normaltension glaucoma? Front Neurol. 2017;8:47.
Salerni F, Repetto R, Harris A, Pinsky P, Prudhomme C, Szopos M, et al. Biofluid modeling of the coupled eyebrain system and insights into simulated microgravity conditions. PLoS ONE. 2019;14: e0216012.
Laurie SS, Vizzeri G, Taibbi G, Ferguson CR, Hu X, Lee SMC, et al. Effects of shortterm mild hypercapnia during headdown tilt on intracranial pressure and ocular structures in healthy human subjects. Physiol Rep. 2017. https://doi.org/10.14814/phy2.13302.
Lee C, Rohr J, Sass A, Sater S, Zahid A, Macias B, et al. In vivo estimation of optic nerve sheath stiffness using noninvasive MRI measurements and finite element modeling. J Mech Behav Biomed Mater. 2020;110: 103924.
Gerlach DA, MarshallGoebel K, Hasan KM, Kramer LA, Alperin N, Rittweger J. MRIderived diffusion parameters in the human optic nerve and its surrounding sheath during headdown tilt. NPJ Microgravity. 2017;3:1–6.
Magnæs B. Body position and cerebrospinal fluid pressure. Part 1: Clinical studies on the effect of rapid postural changes. J Neurosurg. 1976;44:687–97.
Magnæs B. Body position and cerebrospinal fluid pressure. Part 2: clinical studies on orthostatic pressure and the hydrostatic indifferent point. J Neurosurg. 1976;44:698–705.
Magnaes B. Clinical studies of cranial and spinal compliance and the craniospinal flow of cerebrospinal fluid. Br J Neurosurg. 1989;3:659–68.
Holmlund P, Støverud KH, Wåhlin A, Wiklund U, Malm J, Jóhannesson G, et al. Posturedependent collapse of the optic nerve subarachnoid space: a combined MRI and modeling study. Invest Ophthalmol Vis Sci. 2021;62:26.
Killer HE, Berberat J, Wostyn P, Pircher A, Rossinelli D. Posturedependent collapse of the optic nerve subarachnoid space: a combined mri and modeling study. Investig Ophthalmol Vis Sci. 2021;62:5783.
Killer HE, Laeng HR, Flammer J, Groscurth P. Architecture of arachnoid trabeculae, pillars, and septa in the subarachnoid space of the human optic nerve: anatomy and clinical considerations. Br J Ophthalmol. 2003;87:777–81.
Saboori P. Subarachnoid space trabeculae architecture. Clin Anat. 2021;34:40–50.
DeikeHofmann K, Reuter J, Haase R, Paech D, Gnirs R, Bickelhaupt S, et al. Glymphatic pathway of gadoliniumbased contrast agents through the brain: overlooked and misinterpreted. Invest Radiol. 2019;54.
Gehlen M, Kurtcuoglu V, Schmid DM. Is posturerelated craniospinal compliance shift caused by jugular vein collapse? A theoretical analysis. Fluids Barriers CNS. 2017;14:5.
Qvarlander S, Sundström N, Malm J, Eklund A. Postural effects on intracranial pressure: modeling and clinical evaluation. J Appl Physiol. 2013;115:1474–80.
Weed LH. Studies on cerebrospinal fluid. No. III : the pathways of escape from the subarachnoid spaces with particular reference to the Arachnoid Villi. J Med Res. 1914;31:51–91.
Kido DK, Gomez DG, Pavese AM, Potts DG. Human spinal arachnoid villi and granulations. Neuroradiology. 1976;11:221–8.
Eklund A, Smielewski P, Chambers I, Alperin N, Malm J, Czosnyka M, et al. Assessment of cerebrospinal fluid outflow resistance. Med Biol Eng Comput. 2007;45:719–35.
Davson H. Formation and drainage of the cerebrospinal fluid. Sci Basis Med Annu Rev. 1966;238–59.
Czosnyka M, Czosnyka Z, Momjian S, Pickard JD. Cerebrospinal fluid dynamics. Physiol Meas. 2004;25:R5176.
Andersson N, Malm J, Eklund A. Dependency of cerebrospinal fluid outflow resistance on intracranial pressure. J Neurosurg. 2008;109:918–22.
Glimcher SA, Holman DW, Lubow M, Grzybowski DM. Ex vivo model of cerebrospinal fluid outflow across human arachnoid granulations. Investig Ophthalmol Vis Sci. 2008;49:4721–8.
Alitalo K. The lymphatic vasculature in disease. Nat Med. 2011;17:1371–80.
Malm J, Jacobsson J, Birgander R, Eklund A. Reference values for CSF outflow resistance and intracranial pressure in healthy elderly. Neurology. 2011;76:903–9.
Raykin J, Best L, Gleason R, Mulugeta L, Myers J, Nelson E, et al. Optic nerve sheath mechanics and permeability in VIIP syndrome. Invest Ophthalmol Vis Sci. 2014;55:4591.
Støverud KH, Langtangen HP, Ringstad GA, Eide PK, Mardal KA. Computational investigation of cerebrospinal fluid dynamics in the posterior cranial fossa and cervical subarachnoid space in patients with Chiari I malformation. PLoS ONE. 2016;11:1–16.
Wang N, Xie X, Yang D, Xian J, Li Y, Ren R, et al. Orbital cerebrospinal fluid space in glaucoma: the Beijing intracranial and intraocular pressure (iCOP) study. Ophthalmology. 2012;119:2065–74.
Van der Westhuizen J, Prieur Du Plessis J. An attempt to quantify fibre bed permeability utilizing the phase average Navier–Stokes equation. Compos Part A Appl Sci Manuf. 1996;27:263–9.
Westhuizen J, Plessis JPD. Quantification of unidirectional fiber bed permeability. J Compos Mater. 1994;28:619–37.
Tada Y, Naaoshimo T. Modeling and simulation of brain lesions by the finiteelement method. IEEE Eng Med Biol Mag. 1994;13:497–503.
Gupta S, Soellinger M, Boesiger P, Poulikakos D, Kurtcuoglu V. Threedimensional computational modeling of subjectspecific cerebrospinal fluid flow in the subarachnoid space. J Biomech Eng. 2009;131:21010.
Saboori P, Sadegh A. Material modeling of the head’s subarachnoid space. Sci Iran. 2011;18:1492–9.
Ayyalasomayajula A, Park RI, Simon BR, Vande Geest JP. A porohyperelastic finite element model of the eye: the influence of stiffness and permeability on intraocular pressure and optic nerve head biomechanics. Comput Methods Biomech Biomed Engin. 2016;19:591–602.
Jonas JB, Mardin CY, SchlotzerSchrehardt U, Naumann GOH. Morphometry of the human lamina cribrosa surface. Investig Ophthalmol Vis Sci. 1991;32:401–5.
Jonas JB, Holbach L. Central corneal thickness and thickness of the lamina cribrosa in human eyes. Investig Ophthalmol Vis Sci. 2005;46:1275–9.
Lenfeldt N, Koskinen LOD, Bergenheim AT, Malm J, Eklund A. CSF pressure assessed by lumbar puncture agrees with intracranial pressure. Neurology. 2007;68:155–8.
Vinje V, Eklund A, Mardal KA, Rognes ME, Støverud KH. Intracranial pressure elevation alters CSF clearance pathways. Fluids Barriers CNS. 2020;17:29.
Avasthey P. Venous pressure changes during orthostasis. Cardiovasc Res. 1972;6:657–63.
Bateman GA, Bateman AR. Differences in the calculated transvenous pressure drop between chronic hydrocephalus and idiopathic intracranial hypertension. AJNR Am J Neuroradiol. 2019;40:68–73.
Dawson EA, Secher NH, Dalsgaard MK, Ogoh S, Yoshiga CC, Gonza J, et al. Standing up to the challenge of standing : a siphon does not support cerebral blood flow in humans. Am J Physiol Regul Integr Comp Physiol. 2004;287:911–4.
Valdueza JM, von Münster T, Hoffman O, Schreiber S, Einhäupl KM. Postural dependency of the cerebral venous outflow. Lancet. 2000;355:200–1. https://doi.org/10.1016/S01406736(99)048047.
Eklund A, Johannesson G, Johansson E, Holmlund P, Qvarlander S, Ambarki K, et al. The pressure difference between eye and brain changes with posture. Ann Neurol. 2016;80:269–76.
Moller PM. The pressure in the orbit. Acta Ophthalmol Suppl. 1955;1–100.
Holmlund P, Eklund A, Koskinen LOD, Johansson E, Sundstrom N, Malm J, et al. Venous collapse regulates intracranial pressure in upright body positions. Am J Physiol Regul Integr Comp Physiol. 2018;314:R377–85.
Andresen M, Hadi A, Petersen LG, Juhler M. Effect of postural changes on ICP in healthy and ill subjects. Acta Neurochir (Wien). 2014;157:109–13.
Petersen LG, Petersen JCG, Andresen M, Secher NH, Juhler M. Postural influence on intracranial and cerebral perfusion pressure in ambulatory neurosurgical patients. Am J Physiol Regul Integr Comp Physiol. 2016;310:R100–4.
Sunderland DK, Sapra A. Physiology, aqueous humor circulation. Treasure Island (FL); 2022.
Petersen LG, Lawley JS, LiljaCyron A, Petersen JCG, Howden EJ, Sarma S, et al. Lower body negative pressure to safely reduce intracranial pressure. J Physiol. 2019;597:237–48.
Jasien JV, Samuels BC, Johnston JM, Downs JC. Effect of body position on intraocular pressure (IOP), intracranial pressure (ICP), and translaminar pressure (TLP) via continuous wireless telemetry in nonhuman primates (NHPs). Invest Ophthalmol Vis Sci. 2020;61:18.
Berdahl JP, Fautsch MP, Stinnett SS, Allingham RR. Intracranial pressure in primary open angle glaucoma, normal tension glaucoma, and ocular hypertension: a casecontrol study. Invest Ophthalmol Vis Sci. 2008;49:5412–8.
Jonas JB, Wang NL, Wang YX, You QS, Bin XX, Yang DY, et al. Estimated translamina cribrosa pressure difference versus intraocular pressure as biomarker for openangle glaucoma. The Beijing Eye Study 2011. Acta Ophthalmol. 2015;93:7–13.
Pircher A, Remonda L, Weinreb RN, Killer HE. Translaminar pressure in Caucasian normal tension glaucoma patients. Acta Ophthalmol England. 2017;95:e524–31.
Xie X, Zhang X, Fu J, Wang H, Jonas JB, Peng X, et al. Noninvasive intracranial pressure estimation by orbital subarachnoid space measurement: the Beijing Intracranial and Intraocular Pressure (iCOP) study. Crit Care. 2013;17:R162.
Padayachy L, Brekken R, Fieggen G, Selbekk T. Pulsatile dynamics of the optic nerve sheath and intracranial pressure: an exploratory in vivo investigation. Neurosurgery. 2016;79:100–7.
Fleischman D, Kaskar O, Shams R, Zhang X, Olson D, Zdanski C, et al. A novel porcine model for the study of cerebrospinal fluid dynamics: development and preliminary results. Front Neurol. 2019;10:1–7.
Morgan WH, Yu DY, Alder VA, Cringle SJ, Cooper RL, House PH, et al. The correlation between cerebrospinal fluid pressure and retrolaminar tissue pressure. Investig Ophthalmol Vis Sci. 1998;39:1419–28.
Hou R, Zhang Z, Yang D, Wang H, Chen W, Li Z, et al. Intracranial pressure (ICP) and optic nerve subarachnoid space pressure (ONSP) correlation in the optic nerve chamber: the Beijing Intracranial and Intraocular Pressure (iCOP) study. Brain Res. 2016;1635:201–8.
Liu D, Kahn M. Measurement and relationship of subarachnoid pressure of the optic nerve to intracranial pressures in fresh cadavers. Am J Ophthalmol. 1993;116:548–56.
Lindén C, Qvarlander S, Jóhannesson G, Johansson E, Östlund F, Malm J, et al. Normaltension glaucoma has normal intracranial pressure. Ophthalmology. 2017;1–8.
Guy AH, Wiggs JL, Turalba A, Pasquale LR. Translating the low translaminar cribrosa pressure gradient hypothesis into the clinical care of glaucoma. Semin Ophthalmol. 2016;31:131–9.
Johannesson G, Eklund A, Linden C. Intracranial and intraocular pressure at the lamina cribrosa: gradient effects. Curr Neurol Neurosci Rep. 2018;18:25.
Boye D, Montali M, Miller NR, Pircher A, Gruber P, Killer HE, et al. Flow dynamics of cerebrospinal fluid between the intracranial cavity and the subarachnoid space of the optic nerve measured with a diffusion magnetic resonance imaging sequence in patients with normal tension glaucoma. Clin Exp Ophthalmol. 2018;46:511–8.
Lee AG, Mader TH, Gibson CR, Brunstetter TJ, Tarver WJ. Space flightassociated neuroocular syndrome (SANS). Eye (Lond). 2018;32:1164–7.
Mader TH, Gibson CR, Pass AF, Kramer LA, Lee AG, Fogarty J, et al. Optic disc edema, globe flattening, choroidal folds, and hyperopic shifts observed in astronauts after longduration space flight. Ophthalmology. 2011;118:2058–69.
Lawley JS, Petersen LG, Howden EJ, Sarma S, Cornwell WK, Zhang R, et al. Effect of gravity and microgravity on intracranial pressure. J Physiol. 2017;595:2115–27.
Iliff JJ, Wang M, Liao Y, Plogg BA, Peng W, Gundersen GA, et al. A paravascular pathway facilitates CSF flow through the brain parenchyma and the clearance of interstitial solutes, including amyloid beta. Sci Transl Med. 2012;4: 147ra111.
Xie L, Kang H, Xu Q, Chen MJ, Liao Y, Thiyagarajan M, et al. Sleep drives metabolite clearance from the adult brain. Science. 2013;342:373–7.
Tain RW, Bagci AM, Lam BL, Sklar EM, ErtlWagner B, Alperin N. Determination of craniospinal canal compliance distribution by MRI: methodology and early application in idiopathic intracranial hypertension. J Magn Reson Imaging. 2011;34:1397–404.
Bradbury MW, Cole DF. The role of the lymphatic system in drainage of cerebrospinal fluid and aqueous humour. J Physiol. 1980;299:353–65.
Ma Q, Ineichen BV, Detmar M, Proulx ST. Outflow of cerebrospinal fluid is predominantly through lymphatic vessels and is reduced in aged mice. Nat Commun. 2017. https://doi.org/10.1038/s41467017014846.
Jacobsson J, Qvarlander S, Eklund A, Malm J. Comparison of the CSF dynamics between patients with idiopathic normal pressure hydrocephalus and healthy volunteers. J Neurosurg. 2018;1–6.
Ekstedt J. CSF hydrodynamic studies in man. 2. Normal hydrodynamic variables related to CSF pressure and flow. J Neurol Neurosurg Psychiatry. 1978;41:345–53.
Eisenträger A, Sobey I, Czosnyka M. Parameter estimations for the cerebrospinal fluid infusion test. Math Med Biol. 2013;30:157–74.
Wahlin A, Ambarki K, Birgander R, Alperin N, Malm J, Eklund A. Assessment of craniospinal pressurevolume indices. AJNR Am J Neuroradiol. 2010;31:1645–50.
Marmarou A, Shulman K, LaMorgese J. Compartmental analysis of compliance and outflow resistance of the cerebrospinal fluid system. J Neurosurg. 1975;43:523–34.
Sater SH, Sass AM, Seiner A, Natividad GC, Shrestha D, Fu AQ, et al. MRIbased quantification of ophthalmic changes in healthy volunteers during acute 15° headdown tilt as an analogue to microgravity. J R Soc Interface. 2021. https://doi.org/10.1098/rsif.2020.0920.
Funding
Open access funding provided by Umea University. This research project is financially supported by the Swedish National Space Agency (Grant no. 193/17) and the Swedish foundation for Strategic Research.
Author information
Authors and Affiliations
Contributions
PH, KHS, AE conceived the study. All authors contributed to the model development and analysis. All authors participated in the first draft, critically edited and revised the manuscript, read and approved the final version of the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Holmlund, P., Støverud, KH. & Eklund, A. Mathematical modelling of the CSF system: effects of microstructures and posture on optic nerve subarachnoid space dynamics. Fluids Barriers CNS 19, 67 (2022). https://doi.org/10.1186/s12987022003664
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12987022003664
Keywords
 CSF dynamics
 Optic nerve subarachnoid space
 Translaminar pressure
 Glaucoma
 Posture
 Compartmentalization
 Ocular glymphatics
 Numerical modelling