Modeling the impact of thoracic pressure on intracranial pressure

A potential contribution to the progression of Spaceflight Associated Neuro-ocular Syndrome is the thoracic-to-spinal dural sac transmural pressure relationship. In this study, we utilize a lumped-parameter computational model of human cerebrospinal fluid (CSF) systems to investigate mechanisms of CSF redistribution. We present two analyses to illustrate potential mechanisms for CSF pressure alterations similar to those observed in microgravity conditions. Our numerical evidence suggests that the compliant relationship between thoracic and CSF compartments is insufficient to solely explain the observed decrease in CSF pressure with respect to the supine position. Our analyses suggest that the interaction between thoracic pressure and the cardiovascular system, particularly the central veins, has greater influence on CSF pressure. These results indicate that future studies should focus on the holistic system, with the impact of cardiovascular changes to the CSF pressure emphasized over the sequestration of fluid in the spine.

Unique to spaceflight, Spaceflight Associated Neuro-ocular Syndrome (SANS) manifests as a reduction in an astronaut's visual acuity and other anatomical changes suggestive of early-stage ocular disorders, such as optic disc edema, cotton wool spots, and choroidal folds 1,2 .The etiology associated with the manifestation of SANS in some astronauts and not others remains unclear but is likely a multifactorial response to the microgravity and vehicle environment 3,4 .Astronauts experience several well-documented physiological changes when exposed to the microgravity environment, one of the most prominent being the cephalad, or headward, redistribution of fluid 5 .Vascular deconditioning appears following these cephalad fluid redistribution changes, resulting in an incomplete orthostatic response on return to a terrestrial environment with standard gravity (1-G) [6][7][8][9] .NASA and Russian Space agencies utilize a number of inflight countermeasures or prereturn-to-earth procedures, such as fluid loading and exercise 10 , compression garments 10,11 , and lower body negative pressure (LBNP) 12 to mitigate the likelihood of orthostatic intolerance.
The prevalent hypothesis relates the fluid shift contribution to SANS symptoms to either elevated cerebrospinal fluid (CSF) pressure or fluid compartmentalization in the region of the eye during microgravity 1 .Recent studies illustrate the potential for moderately elevated intracranial pressure (ICP) in microgravity to pressure levels between 1-G supine and standing orientations.Unlike the terrestrial environment in which CSF pressure varies due to postural changes, the acute change fluid redistribution in microgravity produces mildly and chronically elevated CSF pressure 13 .From this, it is has been hypothesized that SANS symptoms manifest less from vascular changes and more from the relation of CSF redistribution due to physiological pressure changes, such as reduction in thoracic pressure or changes in CSF and cranial blood volume 14,15 .Specifically, thoracic pressure change in microgravity is postulated to mediate the mechanisms of fluid redistribution in the spinal dural space in a manner similar to mediation seen during standing 16 .Although the volume of the CSF in the spinal cord subarachnoid space, encapsulated by the spinal dural sac is less than 20% of the total CSF volume 17 , the contribution of these CSF redistributions resulting from microgravityassociated physiological changes are currently unresolved.
In the context of this manuscript, we will refer to thoracic pressure as associated with the average pressure in the entire thoracic cavity, including the pleural, pulmonary, blood, interstitial, and air domains.Defined in this manner, direct measurement of the averaged thoracic pressure historically relies on inferential measures, such as esophageal pressure 18 or thoracic impedance 19 .Parabolic flight testing shows that the reduction in esophageal pressure during microgravity portions of the flight exceeds 5 mHg on average, compared to <2 mmHg average reduction in central venous pressure (CVP) 20 .Under these conditions, thoracic pressure reduction is related to both the acute cephalad fluid shift 20 and the removal of so-called tissue compressive forces 21 , typically described as the pressure generated by the weight of the tissue surrounding fluid holding vessels.However, mean arterial pressure (MAP) may also exhibit a steady drop during most of the microgravity portion of parabolic flight, reaching as much as 30 mmHg below its maximum value at the initiation of the microgravity period 22 .This is likely because of hydrostatically affected regulation factors due to cephalad fluid shift.Adding to these observations, parabolic flight studies demonstrate that the average pressure of the intra-jugular vein increases during the microgravity portion of the parabolic flight 23 and the flow may stagnate or reverse during the cardiac cycle in more prolonged microgravity settings 24 .This illustrates the complex interaction of fluid redistribution and other gravity-influenced mechanisms derived from parabolic flight studies and that the time-course measurements have not always reached homeostatic pressure balance in the <20 s microgravity period, further complicating interpretation of these findings for long-duration missions.
Lumped-parameter computational models of the cardiovascular system represent an established tool for investigating the impact of cephalad fluid shift and thoracic pressure in altered gravity conditions [25][26][27][28] .Early models, premised on the work of Guyton 29 , provide an understanding of the importance of the decrease in hydrostatically induced thoracic pressure, and thus a change in intrapleural (extracardiac) pressure, in explaining the paradoxical relation of stroke volume and cardiac output through microgravity-induced compounding of physiological responses 25 .Models that include increasingly complex representations of the cardiovascular blood compartment interactions and regulatory mechanisms relate vascular deconditioning to the onset of syncope and orthostatic intolerance 26,28,30 and the relative contribution of altered (artificial) gravity to preventing orthostatic deconditioning 31 .Fewer lumped-parameter models include intracranial blood and CSF compartments [32][33][34][35] .Those who implement such tools often simplify assumptions, in particular, Monro-Kellie conditions, where a constraint restricts the net fluid volume in the head to a fixed quantity due to the presence of the rigid cranium. 32,33,35or CSF and blood compartment interactions 21 .Regardless of these limitations, reduced gravity and tilt simulations with lumped parameter models that include the intracranial blood, with or without CSF interactions, predict a moderate reduction in ICP between 1 mmHg 35 and 4 mmHg 21 under microgravity conditions with respect to the supine position.Although these findings qualitatively match invasive CSF pressure measurements under acute microgravity conditions (~2-4 mmHg) 13 , they do not include the potential regulation of the CSF pressure via repositioning of fluid to the spinal dural sac via changes in simulated thoracic pressure which may be important at longer durations of microgravity exposure.
In this study, we investigate the thoracic-to-spinal dural sac transmural pressure contribution to CSF regulation proposed by Laurie, et al. 15 through a set of numerical investigations simulating microgravity thoracic pressure change and the role such changes play on CSF pressure regulation.The numerical study extends a published lumped parameter model of CSF and blood interactions in the cranial vault 33 to include Monro-Kellie volume constraints, the existence of multiple cranial drainage pathways, and dynamic venous compliance contributions.We demonstrate the model's credibility in CSF pressure prediction by comparative validation to analogous tilt table studies.To elicit an understanding of the magnitude and pathways by which microgravity thoracic pressure changes contribute to changes in CSF pressure, we present the results of two parametric simulations that seek to isolate specific influences on CSF spaces.Each study complements the other, and together they inform a more comprehensive picture of the role of thoracic pressure in CSF pressure regulation.

Methods
The head model The present mathematical model (henceforth, head model) is an extension of the lumped-parameter model introduced by Stevens, et al. 33 , to study steady-state solutions to CSF Infusion, Head-Down Tilt (HDT), and Microgravity challenges.The purpose of the model described herein is to investigate changes in fluid drainage from the head and changes in the pressures and volumes of cranial compartments in response to prescribed changes in thoracic pressure, MAP, and CVP.Except where otherwise noted, we shall use the pressure of Ventricular CSF compartment and ICP interchangably.Since no human studies were performed and all parameter data is available via published sources, ethical approval was not sought for this work.
Briefly, the head model, illustrated in Fig. 1, abstracts the fluid holding anatomy of the cranial space as a series of discrete, compliant compartments with fixed, resistive flow paths specifying the movement of blood and CSF through the cranial volume.Imposing a conservation of mass constraint at each compartment allows for the simultaneous solution of the flow and volume distribution within each fluid space.
Notable changes from the original Stevens model include partitioning the "Venous Sinus/Jugular Vein" compartment into a series of three "necklevel" compartments and a cranial compartment.Additionally, the "Extraventricular CSF" compartment is partitioned into two new compartments.The new neck and CSF compartments are as follows: 1. Cranial CSF: cranial compartment capturing the volume of CSF present in cerebral cisterns and cranial subarachnoid space; 2. Spinal CSF: "body-level" compartment capturing the volume of CSF present in spinal subarachnoid spaces; 3. Venous Sinus: cranial compartment capturing the blood flow and volume of the venous sinus which communicates with the Cranial CSF and the "neck-level" compartments; 4. Jugular Veins: "neck-level" compartment capturing flow through the left and right internal jugular veins; 5. Secondary Veins: "neck-level" compartment capturing the flow through the vertebral veins and other (smaller) spinal venous structures; 6. Extrajugular: "neck-level" compartment capturing the flow through interstitial spaces, deep neck veins, and any other extra-spinal drainage.
This division of "neck-level" compartments is inspired by the ultrasound analysis completed in refs.36 and 37, MRI reconstructions and assessments given in ref. 38, and descriptions given in ref. 39.The division of cerebral drainage between the three "neck-level" compartments (specified as a percentage of total drainage) was taken from the mean values for healthy control subjects in ref. 38 which are compatible with those given in refs.36,37.In addition to the conservation of mass, we enforce two additional constraints: 1. Pressure in the Thoracic Compartment: the pressure for the Thoracic compartment is specified as a function of time.2. Monro-Kellie constraint: Originally considered by Monro and, later, Kellie 40 , the volume of the cranial compartment is fixed.The Monro-Kellie constraint encompasses all compartments labeled as at the "Head Level" in Fig. 1.
While the first of these simply replaces the relevant dynamic equation, the second represents an additional algebraic constraint imposed upon the system.Therefore, the system of ordinary differential equations is overdetermined.To resolve this, we use techniques from geometric numerical integration 41,42 to project the equations onto a lower-dimensional manifold where the constraints are satisfied.This is equivalent to eliminating equations/variables by solving the algebraic constraints for a particular variable, but in a more general setting.This will be discussed in more detail in "Model construction".

Model construction
To derive the head model, much like in similar lumped-parameter models 32,33,43 , we make the following assumptions: • All fluids (i.e., blood and CSF) are incompressible and isothermal.
• Flow across the blood-brain barrier (i.e. between the Capillaries and Brain compartments) is given by the Starling-Landis Equation: where Q Capillaries,Brain is the flow (mL min −1 ) from the Capillaries compartment to the Brain compartment, K CB is the filtration coefficient (mL min −1 mmHg −1 ), P Capillaries is the pressure (mmHg) in the Capillaries compartment, P Brain is the pressure (mmHg) in the Brain compartment, σ CB is the reflection coefficient, π C is the blood colloid osmotic pressure (mmHg), and π B is the brain interstitial fluid osmotic pressure (mmHg).The values for these parameters can be found in Table 1.The only osmotic forces considered are those due to differences in protein concentration.• When not otherwise specified, all other flows are proportional to the pressure differential between compartments.That is, where Q ij denotes the flow from compartment i to compartment j in mL min −1 , Z ij is the fluidity (or inverse of flow resistance) in mL min −1 mmHg −1 , and P i , P j refer to the pressures in compartments i and j, respectively, in mmHg.• The change in volume between two compartments that share a compliant boundary is linear in the change in the pressure differential between them.That is, where V ij denotes the volume of the cup formed by the interface between compartments i, j in mL and C ij = C ji is the local compliance between the compartments in mL mmHg −1 .
• CSF production is constant as a result of sufficiently robust regulatory mechanisms at the pressure levels of interest 33 .(See Table 2.) By imposing conservation of mass (equiv.volume, since the fluids are assumed incompressible) in each compartment, we form the following system of differential equations: where PðtÞ ¼ P 1 ðtÞ . . .where 0 and I are suitably sized zero and identity matrices, respectively.As in 33 , we assign a 0.2 mmHg difference in pressure between the Spinal CSF compartment and the Ventricular CSF compartment to represent the transmantle pressure.In order for CSF to flow in the stead-state condition from the Ventricular CSF compartment, through the Cranial CSF compartment, to the Spinal CSF compartment, it is necessary for the Cranial CSF compartment to take a baseline pressure between 11.0 and 11.2 mmHg.Our numerical experimentation indicates that the model is not sensitive to the value within this range, so a mean value of 11.1 mmHg is assigned to the baseline Cranial CSF compartment pressure.
Also as in 33 , we use the approximation that the fluidity between the Brain and Ventricular CSF compartments is 1000 times the fluidity across the blood-brain barrier, i.e., between the Capillaries and Brain compartments.We therefore set Z Ventricular CSF,Brain = 1000K CB , where K CB is the filtration coefficient from the Starling-Landis equation.
We compute the remainder of the fluidities using measured mean pressures and flows between compartments.That is, we compute where Q ij is the mean flow rate from compartment i to compartment j, and P i ; P j are the mean pressures for those compartments.This ensures that the given mean pressures are a steady-state solution of the differential equation.The values and their sources used for this paper are given in Tables 2 and 3.
The baseline compliance values at mean pressures are given in Table 4.We use the pressure-dependent cranial compliance functions from 32,43,44 for the compliances between the following compartments: • Brain and Ventricular CSF • Brain and Venous Sinus • Brain and Cranial CSF • Venous Sinus and Cranial CSF.Details and derivation may be found in the original sources, we present the approach in brief below.The general form is given by where we use the parameter values C 0 = 6.5333, r = 0.633431, and γ = 0.604229 32,33,43 .As in 32 , we portion the compliance values by relative volumes and assign the bulk (95%) of compliance values to interfaces with the venous compartments.This yields the following pressure-dependent compliance formulae: For the neck-level (Jugular, Secondary Venous, Extrajugular) and Central Veins compartments, we use the formulation for venous compliance as given by 26 to compute a total compliance value for the compartment.That is, for compartment i, the total compliance is determined by  where C 0 i is the peak compliance of the compartment, N represents the ratio of asymptotic to peak compliance, α i is the compliance steepness factor, and P i,trans is the transmural pressure of the compartment.For the neck-level compartments, the compartment pressure is used for the transmural pressure.For the Central Veins compartment, the transmural pressure is computed as P Central Veins (t) − P Thoracic (t).The peak compliance values, C 0 i , are computed by matching the baseline compliance values from Table 4 with Equation (9) using the mean pressures given in Table 3 and parameter values from Table 1.
To enforce linear constraints such as the Monro-Kellie constraint, we make use of tangent space parameterization techniques from ref. 42.This is equivalent to algebraically eliminating equations using the constraint equations to solve for given variables.By using this parameterization, the constraints will be satisfied to machine precision at every time step.We consider the following general form for a system of ordinary differential equations: subject to the linear constraints Let Q be an orthonormal basis for the null space of A, that is, Q T Q = I and Ay = 0 implies y = Qz for some z.Assuming that AX(0) − B = 0, we can write X(t) = QZ(t) + X(0), where Z(0) = 0.The system in Eq. ( 10) then can be written as the over-determined system: Premultiplication of Eq. ( 12) with Q T yields our reduced set of equations: which may be solved using standard techniques.

Model Validation
To validate the construction and calibration of the head model, we turn to clinical data for tilt-table studies, where participants are placed in a supine position and hemodynamic and CSF pressure data is collected while subjected to a series of tilt angles.While not in a supine position (i.e., at tilt angles other than 0), the body is subject to hydrostatic effects which can induce fluid shifts and pressure changes.We replicated the experimental conditions of three tilt-table studies that included data on ICP (mapped here to the Ventricular CSF pressure) [45][46][47] .The results of these validation are presented in Fig. 2. We note a fairly reasonable agreement with the referent studies with a slight underprediction (in terms of magnitude) for small tilt-angles and a slight overprediction (in terms of magnitude) for large tilt-angles.Considering the head model does not incorporate a lower body or arteriole-regulatory mechanisms that would be present in response to large-scale fluid shifts such as those imparted by hydrostatics, this level of agreement seems reasonable.

Experimental configurations
To help orient the reader, we summarize the changes made for each numerical experiment in Fig. 1.Note that the model is held in the supine position for each experiment, removing hydrostatic effects.Each experiment is run until the model has converged to a steady-state.
For the Fixed MAP and CVP study and to isolate the effects of the compliant boundary between the Thoracic and the Spinal CSF compartments, we make the following changes to the head model: • the pressure in the Thoracic compartment linearly changes from -6 --10 mmHg over the first minute of the experiment, that is, > : ; • the pressure in the Central Arteries compartment is fixed to a specified value, and • the pressure in the Central Veins compartment is fixed to a specified value.
In this study, we consider a total of 91 combinations of (fixed) pressure in the Central Arteries and Central Veins compartments.In particular, we consider a Central Arteries pressure from 86; 87; . . .; 98 f g mmHg (representing ± 6 mmHg from baseline) and a Central Veins pressure from 2; 3; . . .; 8 f gmmHg (representing ± 3 mmHg from baseline).This range of  The model is first initialized with the given pressure for the Central Arteries and Central Veins compartments and allowed to run until steadystate.These steady-state pressures are then used as the initial condition for the experiment in which the pressure in the Thoracic compartment is varied.While this linear pressure profile for the Thoracic compartment does not correspond to any particular experimental or observed profile, the analysis was repeated for various non-linear profiles and over a range of times.Changes to the length of time during which the changes take place as well as non-linear profiles do impact the magnitude of the response: shorter time intervals and "sharper" profiles (i.e., those with larger rates of change) will increase the magnitude of the initial "ring" but do not alter the steadystate solution nor the qualitative behavior of the solution.Therefore, a linear profile was used for ease of implementation as a dynamic equation.
For the Dynamic MAP and CVP study and to determine the impact of changes in pressure in the Thoracic compartment on cardiac output, we make the following change to the head model: • the pressure in the Thoracic compartment linearly changes from -6 to -10 mmHg over the first minute of the experiment, that is, • the cardiac output is as originally described in ref. 25 and used in ref. 26.
That is, the cardiac output is modeled as a mass-balance between pressure-driven flow from the Central Veins to the heart and a Starling-like flow from the heart to the Central Arteries.At each heartbeat, the pressure for a virtual heart compartment, P Heart , is determined so that the following equation is satisfied: F in represents the pressure-driven flow into the heart and is given by where Z Heart represents the fluidity between the Central Veins compartment and the virtual heart compartment and P Collapse represents the pressure of partial collapse 26 .In our application, the partial collapse pressure is given by P Collapse = P Thoracic + 2. F out represents the flow from the heart into the Central Arteries compartment and is given by where S represents the effectiveness of the heart pumping (values S < 1 can be used to model damaged/atrophied tissue), HR(t) is the heart rate (beats min −1 ), SV is the stroke volume (mL beat −1 ), and C is a subjectspecific tuning factor used to fit the cardiac output to the specified baseline flow rate.In our application, the external pressure, P External is equal to the pressure in the Thoracic compartment, P Thoracic .This reflects the impact that the thoracic space pressure has on diastolic filling 48 .Once a virtual heart pressure that achieves this balance is found, the resulting flow rate is used for the duration of the heartbeat ((HR(t)) −1 minutes).
To model the baroreceptor regulatory mechanism, we use the formulation developed in ref. 49.Since the activation amounts are a time-average of the previous cardiac cycle, we may solve 49 Eq. 5 to yield where the notation HR [n] denotes the heart rate over the nth cardiac cycle, τ is the characteristic time for heartrate, and σ [n] is given by where P 0 is the baseline Central Arteries pressure, α σ is the coefficient for sympathetic activity, β σ is the coefficient for parasympathetic activity, γ σ is the baseline activation rate, ν is the slope parameter, and HR 0 is the baseline heart rate.

Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Fixed MAP and CVP
The first computational experiment attempts to isolate the direct effect of changes in the pressure of the Thoracic compartment on the Ventricular CSF compartment pressure by holding the MAP and CVP constant while the pressure in the Thoracic compartment varies.Here we investigate 91 combinations of pressures in the Central Arteries and Central Veins compartments consisting of ± 6 mmHg from baseline mean arterial pressure and ±3 mmHg from baseline central venous pressure.As shown in Fig. 3, these 91 combinations exhibit an initial transient in response to the change in pressure in the Thoracic compartment but the Ventricular CSF compartment pressure then relaxes to a steady state with a maximum change of <0.05 mmHg from the initial condition.

Dynamic MAP and CVP
The second experiment allows the MAP and CVP to be driven by a dynamic cardiac output function as described in "Experimental configurations".This dynamic cardiac output function allows the changes in thoracic pressure to influence the Central Arteries and Central Veins compartments and its  direct interaction with the ventricular CSF compartment.As shown in Fig. 4, at the default MAP and CVP values (see Table 3), we observe an ~2 mmHg drop in ventricular CSF pressure in response to the drop in thoracic pressure from -6 mmHg to10 mmHg.Additionally, we observe an ~2 mmHg increase in MAP and a 3 mmHg decrease in CVP.These changes are accompanied by an increase in cardiac output (not shown in the Figure ) from 5300 mL min −1 to 5750 mL min −1 .

Discussion
The presentation of SANS is a complex system response to microgravity and other environmental factors likely linked to cephalad fluid shift in both the cardiovascular and cerebral spinal fluid domains 1 .Although research continues to elucidate the critical pathways associated with vision impairment, the relative contribution of many of these complex fluid shift-influenced features can be well evaluated by computational simulations 35,50 .One such potentially affecting feature is the change of the spinal dura transmural pressure due to the microgravity-induced alterations of thoracic pressure 15 .
Establishing the true influence of this relationship through direct measurements of CSF and thoracic pressure remains unfeasible in spacefaring humans leading to computational modeling as one of the most viable means of analysis and insight.Towards this end, in establishing plausible mechanistic pathways for sustained CSF drainage from cranial ventricular to subarachnoid spaces, this study has explored two computational simulations with a lumpedparameter model-fixed MAP and CVP as well as dynamic MAP and CVP.The Fixed MAP and CVP study seeks to isolate the influence of an acute change in spinal dura transmural pressure from its effects on other physiological parameters.As illustrated in "Fixed MAP and CVP" and Fig. 3, our model does not support the hypothesis that a decrease in thoracic pressure contributes to a substantial or sustained reduction in ventricular CSF pressure resulting from the isolated sequestration of CSF in the spinal CSF spaces (including the spinal dural sac).The changes in pressure attributable to this compliant interaction between the Thoracic compartment and the Spinal CSF compartment are vanishingly small (less than ± 0.5% of the mean value) and recover quickly over a wide range of cerebral perfusion pressures.These findings suggest that the impact of thoracic pressure on ventricular CSF pressure should not be viewed in isolation.Rather, a holistic view of the effects of changes in thoracic pressure and the accompanying CSF drainage and production pathways is necessary.
When we include the impacts of thoracic pressure-induced alterations in transmural pressure on the central arteries, veins, and heart similar to those seen in acute microgravity 21 , we observe a decrease in Ventricular CSF compartment pressure of ~2 mmHg that is consistent with the drop observed in acute microgravity 22,51 .It should be considered that the Dynamic MAP and CVP study (Dynamic MAP and CVP) does not exactly mirror a microgravity environment and the resulting MAP is relatively constant, increasing by 2.1% from the baseline value.This is accompanied by an increase in the cardiac blood flow from 5300 ml min −1 to 5750 ml min −1 .From these findings of the Dynamic MAP and CVP simulation and those from the Fixed MAP and CVP simulation, we infer that the acute microgravity-induced changes in ICP result predominantly from changes imposed on the venous system (Fig. 4) with little or no influence resulting from transmural pressure changes at the spinal dura.Not surprisingly, the Ventricular CSF compartment pressure is observed to mirror CVP changes during the transition phase, further supporting the link of these two pressures in these simulations.Such observations reinforce the importance of investigating the thoracic pressure influence on the cardiovascular system when examining changes in CSF pressure.
Although this study investigates a wide array of parameters and pathways that may influence the movement and sequestration of CSF in the spinal dura, these results should be interpreted in light of the physiological abstractions made in the computational formulation.Specifically, the formulation uses a lumped parameter abstraction of the fluid system, minimal active regulator mechanisms, and a fixed compliance and resistance values related to some compartments.The range of thoracic pressure changes is premised on the approximation of the thoracic pressure changes as being proportional to changes in esophageal pressure observed in microgravity 20 .This implies the range of thoracic pressure changes in our analysis may be wider than those experienced at the spinal dura.In-beat pulsatility of flow is also neglected which ignores potential small alterations in instantaneous transmural pressure at the spine due to phase differences of the compartment pressure waves.However, these limitations would likely have only a slight impact on the short-term model response and would have minimal to no impact on the observations of this analysis.
In conclusion, we have presented a lumped-parameter model of the blood and CSF flow in the head, neck, and thorax.Using this model, we investigated the role of thoracic pressure on CSF pressure and, in particular, the interactions with the Extraventricular CSF space as suggested by Laurie, et al. 15 .The current modeling effort does not support the hypothesis that this direct interaction with the spinal CSF space is sufficient to explain the experimentally observed drop in CSF pressure.However, we do find numerical evidence supporting that the impact of thoracic pressure on the cardiovascular system, in combination with the direct CSF space interaction, is capable of reproducing the observed changes.This work warrants additional investigation in how microgravity can influence the flow dynamics in those areas.In particular, investigating the unloading of tissue weight in microgravity, as suggested by Buckey, et al. 21, on the compliant Central Veins compartment is future work we intend to explore.

Fig. 1 |
Fig. 1 | Connectivity Diagram of the Head Model with the Changes for the Numerical Experiments.
values covers the baseline Central Arteries and Central Veins compartment pressures used by the models in refs.32,33 in a symmetric interval and is intended to demonstrate the independence of the results from any particular fixed values for the pressures.

Fig. 2 |
Fig. 2 | Comparison of model predictions of ventricular CSF compartment pressure against tilttable studies.Error bars indicate mean ± one standard deviation.

Fig. 3 |
Fig. 3 | Transient response of ventricular CSF compartment pressure to changes in thoracic compartment pressure when MAP and CVP are fixed.The simulations show <0.05 mmHg change from initial CSF pressure over 91 permutations of MAP and CVP.

Fig. 4 |
Fig. 4 | Change in Central Arteries, Central Veins, and Ventricular CSF pressure in response to a drop in Thoracic pressure for dynamic MAP and CVP.

Table 1 |
Parameter Values

Table 2 |
Mean Flows between Compartments at Steady-State

Table 3 |
Baseline Compartment Pressures at Steady-State 32

Table 4 |
Baseline compliance values Note: Since C ij = C ji , the order of the compartments is arbitrary.a Reassigned from the "Venous Sinus/Jugular Veins" compartment.