 Research
 Open Access
 Published:
Interactions of brain, blood, and CSF: a novel mathematical model of cerebral edema
Fluids and Barriers of the CNS volumeÂ 18, ArticleÂ number:Â 42 (2021)
Abstract
Background
Previous models of intracranial pressure (ICP) dynamics have not included flow of cerebral interstitial fluid (ISF) and changes in resistance to its flow when brain swelling occurs. We sought to develop a mathematical model that incorporates resistance to the bulk flow of cerebral ISF to better simulate the physiological changes that occur in pathologies in which brain swelling predominates and to assess the modelâ€™s ability to depict changes in cerebral physiology associated with cerebral edema.
Methods
We developed a lumped parameter model which includes a representation of cerebral ISF flow within brain tissue and its interactions with CSF flow and cerebral blood flow (CBF). The model is based on an electrical analog circuit with four intracranial compartments: the (1) subarachnoid space, (2) brain, (3) ventricles, (4) cerebral vasculature and the extracranial spinal thecal sac. We determined changes in pressure and volume within cerebral compartments at steadystate and simulated physiological perturbations including rapid injection of fluid into the intracranial space, hyperventilation, and hypoventilation. We simulated changes in resistance to flow or absorption of CSF and cerebral ISF to model hydrocephalus, cerebral edema, and to simulate disruption of the bloodâ€“brain barrier (BBB).
Results
The model accurately replicates wellaccepted features of intracranial physiology including the exponentiallike pressureâ€“volume curve with rapid fluid injection, increased ICP pulse pressure with rising ICP, hydrocephalus resulting from increased resistance to CSF outflow, and changes associated with hyperventilation and hypoventilation. Importantly, modeling cerebral edema with increased resistance to cerebral ISF flow mimics key features of brain swelling including elevated ICP, increased brain volume, markedly reduced ventricular volume, and a contracted subarachnoid space. Similarly, a decreased resistance to flow of fluid across the BBB leads to an exponentiallike rise in ICP and ventricular collapse.
Conclusions
The model accurately depicts the complex interactions that occur between pressure, volume, and resistances to flow in the different intracranial compartments under specific pathophysiological conditions. In modelling resistance to bulk flow of cerebral ISF, it may serve as a platform for improved modelling of cerebral edema and bloodâ€“brain barrier disruption that occur following brain injury.
Introduction
Intracranial pressure (ICP) is the cornerstone of most treatment decisions in the neurointensive care unit. Monro and Kellie presented the hypothesis linking volume and pressure in the closed cranium that served as the basis for our understanding of raised ICP [1, 2]. The physiologist Hugh Davson described steadystate ICP in terms of the flow of cerebrospinal fluid, the resistance to its outflow from the cranium, and the pressure in the sagittal sinus where CSF is absorbed [3, 4]. This equation, which defines steadystate ICP under certain conditions laid the cornerstone upon which all further models of intracranial physiology were built. Davsonâ€™s simple, oneelement model of CSF hydrodynamics accounts for the value of steadystate ICP when the primary pathology is a disturbance in the ability to reabsorb CSF but does not account for other pathologies that may lead to increased ICP. The pioneering work of Marmarou added the element of intracranial compliance in order to better describe the dynamic aspects of ICP and the pressure changes that occur with injection of different volumes into the closed cranium at varying rates [5]. Further investigators used these models to describe the importance of the rate of volume changes in the cranium on intracranial pressure dynamics [6, 7]. Others developed models to depict obstruction to the flow of CSF within the ventricular system and replicate various clinical aspects of obstructive hydrocephalus [8, 9]. Subsequently, Ursino and others added the important element of the cerebrovascular circulation to models of intracranial fluid dynamics to account for cerebral blood flow (CBF) and its key aspects such as cerebrovascular autoregulation [10,11,12,13,14,15,16,17]. Recent studies have also modeled the ICP waveform and translated an electrical analogue model into a working physical model of the cerebral vascular circulation [18, 19]. Together these models of intracranial hydrodynamics built our current understanding of ICP, CBF, and their interrelationships and served as the foundation for a vast body of literature describing important interactions between cerebral and systemic physiological parameters [20, 21].
While models of intracranial physiology to date have provided the basis of our understanding of ICP, it is important to point out that fundamentally they were based on two primary elements: the circulation of CSF and the cerebrovascular circulation. As such, they did not take into account a crucial aspect of intracranial physiology, swelling of the brain itself and the effect of raised ICP on the circulation of cerebral interstitial fluid (ISF). Cerebral swelling is a key aspect of many pathologies commonly treated in the neurointensive care unit, including traumatic brain injury (TBI) and stroke. The ability to model the effects of cerebral swelling on intracranial physiology is an important goal. In this work, we aimed to build on the previous classical models of intracranial fluid dynamics adding the additional component of â€śbrainâ€ť to those models that have described the flow of â€śCSFâ€ť and â€śblood.â€ť As such, we sought to develop a model of intracranial hydrodynamics that incorporates the flow of cerebral interstitial fluid. We sought to assess the ability of this model to predict clinically relevant changes in cerebral physiology that are known to occur with cerebral edema and raised ICP.
Methods
We developed a lumped parameter model to represent the interactions between brain tissue, CSF, and blood enclosed within the rigid cranium under both normal physiological conditions and under conditions of brain edema. As in the classical formulations of intracranial physiology, we describe the model using a mathematical formulation that is based on the physical relations in the system and included a representation of these relations using an electrical analog circuit. FigureÂ 1 shows the intracranial anatomy and the corresponding electrical analog model. The model consists of four major compartments and flow conduits. Each compartment contains fluid within an elastic barrier that is represented as an elastic chamber with nonlinear pressureâ€“volume relations. These nonlinear relations may take a form of either an exponent or a polynomial describing the chamberâ€™s transmural pressure as a function of total volume enclosed within that chamber. The selection of exponential or polynomial functions were chosen to accord with accepted intracranial physiological relationships and are detailed below. Importantly, the enclosed volume in some chambers may include the volumes of other elastic chambers that are contained within it. From the outsidein these chambers are: (1) The subarachnoid space that is limited outwardly by the semielastic dura and the rigid skull and inwardly by the pia mater but encloses all intracranial volumes including the brain tissue, ISF, CSF, and the cerebrovascular circulation. (2) The brain compartment through which the cerebral interstitial fluid circulates and is limited outwardly by the pia. (3) The ventricles where CSF is formed and through which it flows, and (4) the cerebral vascular tree. The superior sagittal sinus, owing to the thick dural walls of the sinus, is a separate compartment within the cranium. It is the outflow for venous blood and the site of reabsorption for CSF. It maintains a pressure lower than the ICP under most physiological conditions [22,23,24,25,26,27]. In addition, in the model the spinal thecal sac compartment is interconnected to the intracranial space through the link between the cranial and spinal subarachnoid spaces.
The outermost compartment in the model, the subarachnoid space, is represented as nonlinear elastic chamber in the electrical analog circuit by a capacitor with a nonlinear pressureâ€“volume relationship in the form of an exponential function (Eq.Â 1):
\({ICP}_{SAS}\) is the subarachnoid pressure, \({V}_{tot}\) is the total cranial volume, which is the sum of the brain compartment volume, ventricular volume, subarachnoid space volume, and blood volume, and \({V}_{0}\) is the cranial unstressed volume. Coefficient values \({a}_{s}\), \({b}_{s}\), and \({c}_{s}\) as well as all the coefficients for the equations listed below are detailed in Appendix. As the capacity of the intracranial space to accommodate additional volume is limited, this equation accounts for the exponentiallike rise in ICP when the total volume increases. Previous models have also included a similar equation to account for the exponentiallike rise in pressure when intracranial volume reserve capacity reaches its limit [28, 29].
The brain compartment in turn encloses two additional compartments that are, for the most part, contained within it: the ventricles and the vascular bed. This compartmentâ€™s transmural pressure (\({ICP}_{{BR}_{TM}}\)) as a function of its volume is described by a polynomial (Eq.Â 2):
where \({V}_{BR}\) is the total volume enclosed by the pia, including the volume in the ventricles and in the blood vessels. Coefficient values \({a}_{b}\), \({b}_{b}\), \({c}_{b}\) and \({d}_{b}\) are detailed in Appendix.
Ventricular pressure is defined as follows: external pressure is the \({ICP}_{BR}\) and the ventricular transmural pressure \(IC{P}_{{VEN}_{TM}}\) is defined as a function of its own volume as described by a polynomial (Eq.Â 3):
where \({V}_{VEN}\) is the ventricular volume. Coefficient values \({a}_{v}\), \({b}_{v}\), \({c}_{v}\) and \({d}_{v}\) are detailed in Appendix.
The cerebrovascular circulation within the cranium is subdivided into three volume components: arterial and arteriolar, C_{A}, capillary, C_{C} and venous, C_{V}. The pressures in the three subcompartments, arteriolar \({P}_{ARL}\), capillary \({P}_{C}\), and venous \({P}_{V}\) were described as a function of their blood volumes \({V}_{A}\), \({V}_{C}\), and \({V}_{V}\) respectively, by linear expressions (Eq.Â 4):
As noted above, the superior sagittal sinus is the only intracranial compartment that maintains a pressure lower than ICP and is the site for both venous outflow and CSF reabsorption.
The spinal thecal sac is connected to the intracranial space through fluid conduit between the intracranial and spinal subarachnoid spaces. Pressureâ€“volume relations in the sac were also described by a polynomial (Eq.Â 5):
where \({V}_{SPI}\) is the spinal thecal sac volume. Coefficient values \({a}_{t}\), \({b}_{t}\), and \({C}_{t}\) are detailed in Appendix.
Importantly, in the model transmural pressures of the ventricular and brain compartments are very small; however, they are required to maintain compartment volumes. For example, the ventricle is enclosed by the brain on its outer margin. The pressure in the ventricle (P_{VEN}) is slightly greater than the pressure in the brain at baseline physiological conditions, thereby maintaining a positive transmural pressure and ventricular volume. When this normal physiological condition is reversed by pathophysiological processes, collapse of the compartment occurs.
Production rates of CSF (Q_{CSF}) and ISF (Q_{ISF}) in the model are constant (Table 1), in line with experimental evidence that supports a constant rate of production independent of changes in ICP [30]. Initial baseline CBF (Q_{CBF}) is 910Â ml/min, but varies in relation to changing volume and pressure in the other model compartments or changes in resistance within the cerebrovascular tree.
Fluid conduits connect between some of the different compartments. In the model, resistors represent viscous resistance to fluid flows within the brain and ventricular system and through the cerebrovascular tree. In the cerebrovascular circulation, R_{A}, R_{C}, and R_{V} and represent viscous resistance to flow in the circulation. The resistors in the vascular bed are affected by the transmural pressure of the vascular wall that is exposed to ICP_{BR}. Changes in ICP_{BR} affect the volume inside the vessel that inturn alters the resistance. Radius changes affects resistance, according to HagenPoiseuille's law for a cylinder (Eq.Â 6):
\(where V=\pi {r}^{2}L\)
\(\mu\) is viscosity, L segment length and r is the radius. \({R}_{0}\) is initial resistance of the compartment with compartment volume \(V={V}_{0}\) at baseline conditions.
The normally very high resistance to the passage of fluid from the cerebrovascular circulation into the brain parenchyma is depicted by R_{BBB} which has a very high value under normal physiological conditions. The model depicts pathophysiological disturbance that impair the normally nearly impenetrable bloodâ€“brain barrier as a decreasing R_{BBB.}
The resistances to flow of CSF and ISF are depicted as follow: R_{VEN} represents resistance to flow within the ventricles and out of the ventricular system into the subarachnoid space. R_{BF} represents resistance to the bulk flow of cerebral ISF within the brainâ€™s extracellular space and out to the subarachnoid space where it merges with the circulating CSF. R_{OUT} represents the resistance to CSF absorption into the superior sagittal sinus. Lastly, R_{FM} represents the resistance to flow of CSF from the intracranial subarachnoid compartment to the spinal subarachnoid compartment. The resistors to the flow of CSF and ISF are affected by the volume of CSF or ISF in a compartment through which they flow. For example, an extracellular compartment with large fluid content will result in larger space through which cerebral ISF can flow leading to a lower resistance (low R_{BF}). In contrast, when the volume of the extracellular space is low as occurs in brain swelling, the resulting denser cellular space leads to greater resistance to cerebral ISF flow through the constricted ECF channels (high R_{BF}).
The solution of the model involves solving seven differential equations. Due to the nonlinear behavior of the parametersâ€™ values and their dependency on pressure and volume values, the Euler method was employed to solve the model and reach a steady state solution for each pathological condition. The Euler method was used as follows:
Volume update for every capacitor:
where n is capacitor index and \(i\) is discrete time.
Pressure update:
where the function is one of those described above (exponential or polynomial)
Parameter values update:
Flow calculation:
Table 1 details initial baseline values and the model componentsâ€™ characteristic pressure and volume relations. Most baseline values were designated based on normal values obtained from the literature, while others whose values may be uncertain due to the paucity of data in the literature were derived and estimated on the basis of normal physiological pressures and flows.
Simulation of physiological perturbations
Selected physiological perturbations were assessed to test model validity and to assess whether model simulations accorded with expected cerebral physiological responses. We simulated rapid intracranial injection of increasing volumes of fluid at normal ICP in order to establish the pressureâ€“volume response to a rapid change in intracranial volume. In addition, we simulated hyperventilation and hypoventilation by altering the arterialarteriolar resistance (\({R}_{A}\)) in the cerebrovascular component of the model, since this is the primary mechanism by which ventilation changes affect CBF. We assessed the modelâ€™s ability to appropriately simulate the changes in cerebral physiological parameters that are induced by different pathophysiological phenomena that lead to raised ICP, including hydrocephalus (obstructive and absorptive) and cerebral edema induced by either cellular swelling or bloodâ€“brain barrier breakdown. Ventricular obstruction was modelled by an increasing resistance to flow in the ventricles (\({R}_{VEN}\)), absorptive hydrocephalus was modelled by increasing an increasing resistance to CSF reabsorption (\({R}_{OUT}\)), cellular swelling with an increasing resistance to the flow of cerebral interstitial fluid was modelled by increasing the resistance to the bulk flow of cerebral interstitial fluid (\({R}_{BF}\)), while a breakdown in the bloodâ€“brain barrier was simulated by a decrease in the resistance to flow across the BBB (\({R}_{BBB}\)). Changes in volume distributions across the different intracranial compartments (brain tissue, ventricles, subarachnoid space, and the vascular compartments) induced by these alterations once the model reached steadystate were determined. The changes in the ICP waveform with rising ICP were also investigated. MATLAB (The MathWorks Inc., Natick, Massachusetts) was used for intracranial fluid dynamics model simulations.
Results
Initial assessment of the model demonstrated stable values for pressure, volume, and flow in all three model elements of brain, CSF, and blood. The hierarchy of pressures between the ventricles, cerebral parenchyma, and subarachnoid space on the one hand, and the cerebral venous compartment on the other, was maintained in the expected manner (Fig.Â 2). Importantly, simulations of rapid intracranial volume injections, demonstrated the exponentiallike intracranial pressureâ€“volume curve (Fig.Â 3A). Rising ICP_{BR} also resulted in changes to the ICP waveform that are typically seen when ICP is elevated, including a progressively increasing ICP pulse amplitude typical to states of high intracranial elastance (Fig.Â 3B). Changes in arteriolar resistance led to the expected changes in CBF and ICP_{BR}, simulating the decrease in ICP induced by hyperventilation and the increase in ICP induced by hypoventilation (Fig.Â 4). Hyperventilation resulted in a 20% decrease in CBF from baseline and an 18% decrease in ICP. Conversely, hypoventilation resulted in a 30% increase in CBF from baseline and a 15% increase in ICP.
As expected, the lumped parameter model demonstrated rising intracranial pressure under conditions of obstructive hydrocephalus, absorptive hydrocephalus, and cerebral edema (Fig.Â 5Aâ€“D). While ICP rises in each of these cases, the effect on the volumes in the different intracranial compartments varies between simulations of different pathophysiological conditions. When nonabsorptive hydrocephalus (high R_{OUT}) or obstructive hydrocephalus (high R_{VEN}) is simulated, ventricular volume increases substantially to 140% and 152% of baseline values, respectively, as ICP reaches 20Â mmHg. In contrast, when cerebral edema is simulated by an increased R_{BF} or decreased R_{BBB}, ventricular volume decreases markedly to 8% and 12% of baseline values, respectively, as ICP reaches 20Â mmHg.
Table 2 details the changes in model output parameters that are induced by a change in the input values of R_{OUT}, R_{VEN}, R_{BF}, and R_{BBB} which lead to an ICP rise to a level of 30Â mmHg. The changes in the volumes of the brain, ventricles, and subarachnoid space compartments are detailed in the Table. Importantly, in the model CBF decreases substantially with the rise in ICP. The results demonstrate that model output corresponds to accepted patterns associated with each of the distinct pathophysiological conditions. The model clearly indicates that simulating cerebral edema with an increase in the resistance to the flow of cerebral interstitial fluid (\({R}_{BF}\)) leads not only to raised ICP_{BR} (Fig.Â 5C), but also to an increase in total brain volume, and a marked decrease in ventricular and subarachnoid space volume (Table 2), accurately replicating the physiological changes known to occur with brain edema. We also modeled the effect of changes in the integrity of the bloodâ€“brain barrier on cerebral physiological parameters (Table 2). Not surprisingly, a breakdown of bloodâ€“brain barrier integrity represented by a decrease in \({R}_{BBB}\), resulted in an increase in the volume of the brain component in the model (Table 2) leading to a rise in ICP_{BR} (Fig.Â 5D) and a marked decrease in ventricular volume (Fig.Â 6D) that contrasts with the increasing ventricular volume that occurs with nonabsorptive or obstructive hydrocephalus (Fig.Â 6A, B). ICP_{BR} increased dramatically as R_{BBB} decreased to low levels indicating an exponentiallike effect of a rising permeability to fluid shifts across the BBB (Fig.Â 5D).
Discussion
The fourcompartment model we describe simulates the important interactions of brain, blood, and CSF within the closed cranium and accurately replicates important and wellaccepted aspects of intracranial physiology. The model replicates the pressureâ€“volume curve in response to rapid volume injection into the intracranial space (Fig.Â 3) and the responses to hyper and hypoventilation (Fig.Â 4). Importantly, the incorporation of a novel compartment to the model, the brain, through which flows cerebral interstitial fluid, adds an essential component that has not been previously included in models of intracranial hydrodynamics. This is in effect a modification of Davsonâ€™s original formulation:
where the superscript, ^{D}, denotes the parameters referred to by Davson: ICP^{D} a lumped ICP for a single intracranial compartment, Q^{D}_{CSF} representing flow of all CSF within the cranium, which in the context of the current model corresponds to flow of CSF in the subarachnoid space or the sum of ISF and CSF flow, and R^{D}_{OUT} representing the resistance to reabsorption of CSF into the superior sagittal sinus. Rather than treating the intracranial space as a single lumped compartment with a unitary pressure, the model we describe stresses the small but important differences between pressures in each of the different cerebral compartments, ventricles, brain, and subarachnoid space. In our model, steadystate pressure in the brain parenchyma and subarachnoid space, respectively, may be described as follows:
Substituting Eq.Â 9 in Eq.Â 8 yields Eq.Â 10.
This modification of the Davson equation considers both CSF and ISF production and flow as essential elements in determining steadystate cerebral intraparenchymal pressure. This additional element is vital because it allows modelling of pathologies in which cerebral edema is the primary component leading to raised ICP and an increased resistance to bulk flow of ISF. Our results clearly demonstrate the rise in ICP caused by an increase in R_{BF} (Fig.Â 5) and the consequent changes in volumes of the other intracranial compartments that occur with brain swelling (Table 2). Pathologies in which cerebral edema predominates, either as the primary insult or the consequence of secondary brain injury occur frequently in neurocritical care and continue to present difficult therapeutic challenges. The current model offers the ability to separate two distinct pathways leading to ICP elevation, those caused by obstruction to CSF flow (\({R}_{VEN}\), \({R}_{OUT}\)) and those caused by swelling of the brain itself which leads to an increased resistance to cerebral interstitial fluid flow (\({R}_{BF}\)). The model also depicts the changes caused by a breakdown of the normally nearly impenetrable bloodâ€“brain barrier (\({R}_{BBB}\)) (Figs.Â 5D, 6D). The ability to model these important pathophysiological conditions may provide potential benefit by affording an improved categorization of the underlying pathophysiology that leads to raised ICP. Crucially, in addition to demonstrating raised ICP, the model replicates well known clinical manifestations that occur with brain swelling such as a marked decrease in volume of the ventricles and the subarachnoid space (analogous to cisternal and sulcal effacement) which most previous models of intracranial hydrodynamics have not addressed (Table 2). Opening of the bloodâ€“brain barrier is another wellrecognized pathological process known to lead to severe cerebral edema and raised ICP [31,32,33]. The current model is successful in emulating the multifaceted effects of increased BBB permeability, indicating its potential utility in modelling the bloodâ€“brain barrier disruption that occurs with brain injury. Importantly, the model simulates the oftencomplex multivariate interactions between key cerebral physiological parameters. By simulating a primary disturbance in one of four key resistances (\({R}_{VEN}\), \({R}_{OUT}\), \({R}_{BF}\), and \({R}_{BBB}\)) the model is able to demonstrate the ensuing changes in the volumes, pressures, and resistances in each of the intracranial compartments. The model also depicts the changes in intracranial compliance that occur with rising ICP and with a rapid volume load to the intracranial space (Fig.Â 3A). Similar to previous models of ICP pressure dynamics it replicates the changes in the ICP waveform that result from raised ICP and decreased intracranial compliance (Fig.Â 3B) [11, 34]. The ability to study interrelationships is an important advantage of lumped parameter modelling. Although this is a preliminary study to establish the theoretical underpinnings of modelling intracranial pressureâ€“volume dynamics in relation to the resistance to cerebral ISF flow, our hope is that future translational studies will be able to build on this work and use large animal experimental models or data from brain injured patients in order to assess the modelâ€™s ability to help predict thresholds of intracranial volumereserve capacity in cerebral edema.
Previous models of intracranial hydrodynamics and cerebrovascular circulation established the underlying principles needed to depict the interactions between the pressure, volume, and flow of CSF and CBF [35]. However, a substantial gap in these models remained in that none described the important changes that occur when the brain itself swells and the flow of cerebral interstitial fluid is impaired by that swelling. This lacunae in the models of intracranial hydrodynamics is of crucial importance because many of the pathologies dealt with by neurosurgeons and neurocritical care specialists involve severe cerebral edema. Without a means to accurately represent the changes that occur when the brain swells a substantial gap exists in the ability of models of intracranial physiology to describe the interactions between cerebral physiological parameters in a clinically meaningful way. The physiological importance of the flow of cerebral interstitial fluid has received renewed attention in recent years [36,37,38]. When cerebral intracellular swelling occurs the normal flow of brain ISF may be impaired [39, 40] with important consequences for ICP. Intracellular swelling leads to a substantial reduction in the volume of the cerebral extracellular space [41,42,43], leaving smaller constricted channels through which the cerebral ISF must flow. Under these conditions, the flow of cerebral ISF is impaired by cellular swelling. The model depicts this pathophysiological condition as an increased resistance to bulk flow of ISF through the brain parenchyma (increased R_{BF}). An increasing R_{BF} value represents the heightened resistance to the efflux of the water that is a byproduct of cerebral metabolism out of the brain parenchyma by way of its normal passage from the extracellular space into the subarachnoid space. The effects of the resulting increased resistance to cerebral ISF flow are apparent in the model where a marked rise in ICP results from an increase in \({R}_{BF}\). The resulting rise in ICP may, in turn, further affect other cerebral physiological parameters leading to a cascade effect with deleterious consequences. An advantage of the lumped parameter model is that it allows the simulation of these interactions.
The model we describe also needs to be examined in light of recent controversies regarding cerebral ISF production and flow. In recent years, the development of the glymphatic and intramural periarterial drainage (IPAD) theories has challenged traditional conceptions of ISF flow within the brain [44,45,46,47,48]. Strongly held views are held by both sides of the ongoing debate regarding the experimental evidence and whether it supports the glymphatic theory or more traditional theories of ISF production and circulation. We sought to describe a model that describes the importance of ISF circulation as a key component in determining ICP, especially under conditions of brain edema. We aimed to ensure that the model is compatible with both traditional and glymphatic theories of cerebral ISF circulation. The model only assumes that ISF is produced at a constant rate (0.00083Â ml/sec) as noted in Table 1. This comprises just under 20% of the summed total of CSF plus ISF production, which is in line with ratio of production of metabolic water produced by cellular metabolism in relation to the rate of CSF production. Whether this fluid is produced by mechanisms proposed by the glymphatic theory, the IPAD theory, or by traditional theories of ISF production and flow does not bear directly on the modelâ€™s results. The model postulates that cerebral ISF flows from the extracellular space of the brain parenchyma into the subarachnoid space. An important aspect of the model is that it obeys the principles of hydrodynamics in that fluid flows from a higher pressure into a compartment of lower pressure. The subarachnoid space which envelopes the brain is postulated to have a slightly lower pressure than both the brain parenchyma which it surrounds and the ventricles which empty directly into it. The slightly lower pressure of the intracranial subarachnoid space results from the fact that it itself is the site of CSF exit out of the cranium by way of reabsorption into the superior sagittal sinus. Importantly, when cerebral swelling or hydrocephalus lead to a smaller subarachnoid space, resistance to flow within this compartment increases and CSF reabsorption is impaired. While it is difficult to definitively establish whether small pressure differences between the intracranial compartments exist, it is not unreasonable to postulate that they may occur. The small number of clinical studies that have measured ICP concomitantly in the ventricles and brain parenchyma have found differences in these simultaneous measurements [49,50,51,52,53]. A metaanalysis of these studies identified that 70% of concomitant measurements in the same patient could vary by as much asâ€‰Â±â€‰6Â mmHg and an overall mean pressure difference of 1.6Â mmHg exists [54]. Admittedly, it remains difficult to ascertain to what degree these differences are due to actual differences in pressure between compartments and to what degree they may be due to measurement error. Further meticulous experimental and clinical studies will be required to establish whether the small differences in transmural pressures proposed in the model can be demonstrated and quantified by experimental data.
Limitations
Our study has several limitations. This is a preliminary work that seeks to add the additional element of cerebral interstitial fluid flow to the previous models that described ICP in terms of CSF hydrodynamics and the cerebral vascular circulation alone. In order to study the important firstorder effects of adding the novel element of cerebral interstitial fluid flow, we did not incorporate a dynamic model of CBF with autoregulation into the current version of the model. This will be an important goal of future studies that seek to investigate the important interplay between cerebral edema, raised ICP, and the compensatory ability to maintain adequate cerebral perfusion. Second, the current model does not inherently separate the cerebral intracellular and the cerebral extracellular compartments. Further modifications to the model may achieve an improved description of different pathophysiological conditions if this separation is incorporated. Another potential limitation of this study is that some physiological parameters incorporated into the model, such as resistance to flow of CSF in the ventricles and resistance to cerebral ISF flow, are difficult to measure in clinical practice. The estimation of these parameters, which relies on derivations from more easily measured pressures, volumes, and flows in the intracranial system, introduces uncertainty into the model. Though unlikely to alter model behavior, the exact degree to which these uncertainties may modify model output will require further investigation. Lastly, the clinical utility of the model can only be truly assessed with translational studies that incorporate the collection of cerebral physiological parameters from either a large animal model or brain injured patients that are â€śplugged into the modelâ€ť in order to assess its ability to predict thresholds of intracranial volumereserve capacity. Although we did not have the resources to perform these resourceintensive studies, our hope is that future investigations will aim to determine whether the model may be a useful platform to achieve this goal.
Conclusions
The fourcompartment model of intracranial physiology incorporates a depiction of cerebral interstitial fluid flow and resistance that allows modelling of brain edema. The model accounts for the complex interactions that occur between different cerebral physiological parameters and accurately depicts the changes in pressure, volume, and resistances to flow in the different intracranial compartments under specific pathophysiological conditions. It may serve as a platform for improved modelling of the cerebral edema and bloodâ€“brain barrier disruption that occur following brain injury.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Abbreviations
 \({P}_{A}\) :

Arterial pressure
 \({P}_{ARL}\) :

Arteriolar pressure
 \({P}_{C}\) :

Capillary pressure
 \({P}_{V}\) :

Venous pressure
 \({ICP}_{VEN}\) :

Ventricular pressure
 \({ICP}_{BR}\) :

Brain intraparenchymal pressure
 \({ICP}_{SAS}\) :

Subarachnoid space pressure
 \({P}_{SPI}\) :

Spinal thecal pressure
 \({P}_{SSS}\) :

Superior sagittal sinus pressure
 \({C}_{A}\) :

Arterialarteriolar compliance
 \({C}_{C}\) :

Capillary compliance
 \({C}_{V}\) :

Venous compliance
 \({C}_{VEN}\) :

Ventricular compliance
 \({C}_{BR}\) :

Brain compliance
 \({C}_{SAS}\) :

Subarachnoid space compliance
 \({C}_{SPI}\) :

Spinal thecal sac compliance
 \({R}_{A}\) :

Arterialarteriolar resistance to flow
 \({R}_{C}\) :

Capillary resistance to flow
 \({R}_{V}\) :

Venous resistance to flow
 \({R}_{BV}\) :

Bridging veins resistance to flow
 \({R}_{BBB}\) :

Brain Blood Barrier
 \({R}_{VEN}\) :

Resistance to CSF flow in the ventricles
 \({R}_{BF}\) :

Resistance to bulk flow of cerebral interstitial flow
 \({R}_{OUT}\) :

Resistance to CSF reabsorption and to bulk flow in the subarachnoid space
 \({R}_{FM}\) :

Resistance to CSF flow through the foramen magnum
 \({V}_{A}\) :

Arteriolar volume
 \({V}_{C}\) :

Capillary volume
 \({V}_{V}\) :

Venous volume
 \({V}_{VEN}\) :

Ventricular volume
 \({V}_{BR}\) :

Brain volume
 \({V}_{SAS}\) :

Subarachnoid space volume
 \({V}_{TOT}\) :

Total intracranial volume
 \({V}_{SPI}\) :

Spinal thecal sac volume
 \({Q}_{CSF}\) :

CSF production rate
 \({Q}_{ISF}\) :

Cerebral interstitial fluid production rate
 \({Q}_{CBF}\) :

Cerebral blood flow
References
Kellie G. An account of the appearance observed in the dissection of two or three individuals presumed to have perished in the storm of the third and whose bodies were discovered in the vicinity of Leith on the morning of the fourth November 1821, with some reflections on the pathology of the brain. Trans Med Chir Soc (Edin). 1824;I:84â€“122 abd 3â€“69.
Monro A. Observations on the structure and function of the nervous system. Edinburgh: Creech and Johnston; 1783.
Davson H. Physiology of the ocular and cerebrospinal fluid. London: Churchill; 1956.
Davson H, Hollingsworth G, Segal MB. The mechanism of drainage of the cerebrospinal fluid. Brain. 1970;93(4):665â€“78.
Marmarou A, Shulman K, Rosende RM. A nonlinear analysis of the cerebrospinal fluid system and intracranial pressure dynamics. J Neurosurg. 1978;48(3):332â€“44.
Anile C, Portnoy HD, Branch C. Intracranial compliance is timedependent. Neurosurgery. 1987;20(3):389â€“95.
Sullivan HG, Miller JD, Griffith RL 3rd, Carter W Jr, Rucker S. Bolous versus steadystate infusion for determination of CSF outflow resistance. Ann Neurol. 1979;5(3):228â€“38.
Rekate HL. The usefulness of mathematical modeling in hydrocephalus research. Childs Nerv Syst. 1994;10(1):13â€“8.
Rekate HL, Brodkey JA, Chizeck HJ, el Sakka W, Ko WH. Ventricular volume regulation: a mathematical model and computer simulation. Pediatr Neurosci. 1988;14(2):77â€“84.
Ursino M. A mathematical study of human intracranial hydrodynamics. Part 2â€“Simulation of clinical tests. Ann Biomed Eng. 1988;16(4):403â€“16.
Ursino M. A mathematical study of human intracranial hydrodynamics. Part 1â€”The cerebrospinal fluid pulse pressure. Ann Biomed Eng. 1988;16(4):379â€“401.
Ursino M, Iezzi M, Stocchetti N. Intracranial pressure dynamics in patients with acute brain damage: a critical analysis with the aid of a mathematical model. IEEE Trans Biomed Eng. 1995;42(6):529â€“40.
Hoffmann O. Biomathematics of intracranial CSF and haemodynamics. Simulation and analysis with the aid of a mathematical model. Acta Neurochir Suppl (Wien). 1987;40:117â€“30.
Ursino M. Computer analysis of the main parameters extrapolated from the human intracranial basal artery blood flow. Comput Biomed Res. 1990;23(6):542â€“59.
Ursino M, Lodi CA. A simple mathematical model of the interaction between intracranial pressure and cerebral hemodynamics. J Appl Physiol. 1997;82(4):1256â€“69.
Lodi CA, Ter Minassian A, Beydon L, Ursino M. Modeling cerebral autoregulation and CO2 reactivity in patients with severe head injury. Am J Physiol. 1998;274(5):H1729â€“41.
Kadas ZM, Lakin WD, Yu J, Penar PL. A mathematical model of the intracranial system including autoregulation. Neurol Res. 1997;19(4):441â€“50.
Ficola A, Fravolini ML, Anile C. A Physical model of the intracranial system for the study of the mechanisms of the cerebral blood flow autoregulation. IEEE Access. 2018;6:67166â€“75.
Unnerback M, Ottesen JT, Reinstrup P. Validation of a mathematical model for understanding intracranial pressure curve morphology. J Clin Monit Comput. 2020;34(3):469â€“81.
Czosnyka M, Piechnik S, Richards HK, Kirkpatrick P, Smielewski P, Pickard JD. Contribution of mathematical modelling to the interpretation of bedside tests of cerebrovascular autoregulation. J Neurol Neurosurg Psychiatry. 1997;63(6):721â€“31.
Czosnyka M, Smielewski P, Timofeev I, Lavinio A, Guazzo E, Hutchinson P, et al. Intracranial pressure: more than a number. Neurosurg Focus. 2007;22(5):E10.
Sahar A, Hochwald GM, Ransohoff J. Cerebrospinal fluid and cranial sinus pressures. Relationship in normal and hydrocephalic cats. Arch Neurol. 1970;23(5):413â€“8.
Shulman K, Yarnell P, Ransohoff J. Dural Sinus Pressure. In normal and hydrocephalic dogs. Arch Neurol. 1964;10:575â€“80.
Ekstedt J. CSF hydrodynamic studies in man. 2. Normal hydrodynamic variables related to CSF pressure and flow. J Neurol Neurosurg Psychiatry. 1978;41(4):345â€“53.
Fargen KM, Garner RM, Kittel C, Wolfe SQ. A descriptive study of venous sinus pressures and gradients in patients with idiopathic intracranial hypertension. J Neurointerv Surg. 2020;12(3):320â€“5.
Lee K, Kittel C, Aldridge JB, Wolfe SQ, Brown P, Fargen KM. Correlation between intracranial pressure and venous sinus pressures in patients undergoing cerebral venography and manometry. J Neurointerv Surg. 2021. https://doi.org/10.1136/neurintsurg2020017161.
Nakagawa Y, Tsuru M, Yada K. Site and mechanism for compression of the venous system during experimental intracranial hypertension. J Neurosurg. 1974;41(4):427â€“34.
Takemae T, Kosugi Y, Ikebe J, Kumagai Y, Matsuyama K, Saito H. A simulation study of intracranial pressure increment using an electrical circuit model of cerebral circulation. IEEE Trans Biomed Eng. 1987;34(12):958â€“62.
Agarwal GC, Berman BM, Stark L. A lumped parameter model of the cerebrospinal fluid system. IEEE Trans Biomed Eng. 1969;16(1):45â€“53.
Heisey SR, Held D, Pappenheimer JR. Bulk flow and diffusion in the cerebrospinal fluid system of the goat. Am J Physiol. 1962;203:775â€“81.
Kimelberg HK. Current concepts of brain edema. Review of laboratory investigations. J Neurosurg. 1995;83(6):1051â€“9.
Unterberg AW, Stover J, Kress B, Kiening KL. Edema and brain trauma. Neuroscience. 2004;129(4):1021â€“9.
Kimelberg HK. Water homeostasis in the brain: basic concepts. Neuroscience. 2004;129(4):851â€“60.
Linninger AA, Xenos M, Sweetman B, Ponkshe S, Guo X, Penn R. A mathematical model of blood, cerebrospinal fluid and brain dynamics. J Math Biol. 2009;59(6):729â€“59.
Wakeland W, Goldstein B. A review of physiological simulation models of intracranial pressure dynamics. Comput Biol Med. 2008;38(9):1024â€“41.
Abbott NJ. Evidence for bulk flow of brain interstitial fluid: significance for physiology and pathology. Neurochem Int. 2004;45(4):545â€“52.
Sykova E, Nicholson C. Diffusion in brain extracellular space. Physiol Rev. 2008;88(4):1277â€“340.
Thrane AS, Rangroo Thrane V, Nedergaard M. Drowning stars: reassessing the role of astrocytes in brain edema. Trends Neurosci. 2014;37(11):620â€“8.
Iliff JJ, Chen MJ, Plog BA, Zeppenfeld DM, Soltero M, Yang L, et al. Impairment of glymphatic pathway function promotes tau pathology after traumatic brain injury. J Neurosci. 2014;34(49):16180â€“93.
Hrabetova S, Nicholson C. Contribution of deadspace microdomains to tortuosity of brain extracellular space. Neurochem Int. 2004;45(4):467â€“77.
Hansen AJ, Olsen CE. Brain extracellular space during spreading depression and ischemia. Acta Physiol Scand. 1980;108(4):355â€“65.
Hossmann KA. Cortical steady potential, impedance and excitability changes during and after total ischemia of cat brain. Exp Neurol. 1971;32(2):163â€“75.
Matsuoka Y, Hossmann KA. Brain tissue osmolality after middle cerebral artery occlusion in cats. Exp Neurol. 1982;77(3):599â€“611.
Diem AK, Carare RO, Weller RO, Bressloff NW. A control mechanism for intramural periarterial drainage via astrocytes: How neuronal activity could improve waste clearance from the brain. PLoS ONE. 2018;13(10):e0205276.
Nakada T, Kwee IL. Fluid dynamics inside the brain barrier: current concept of interstitial flow, glymphatic flow, and cerebrospinal fluid circulation in the brain. Neuroscientist. 2019;25(2):155â€“66.
Jessen NA, Munk AS, Lundgaard I, Nedergaard M. The glymphatic system: a beginnerâ€™s guide. Neurochem Res. 2015;40(12):2583â€“99.
Plog BA, Nedergaard M. The glymphatic system in central nervous system health and disease: past, present, and future. Annu Rev Pathol. 2018;13:379â€“94.
Brinker T, Stopa E, Morrison J, Klinge P. A new look at cerebrospinal fluid circulation. Fluids Barriers CNS. 2014;11:10.
Schickner DJ, Young RF. Intracranial pressure monitoring: fiberoptic monitor compared with the ventricular catheter. Surg Neurol. 1992;37(4):251â€“4.
Gopinath SP, Robertson CS, Contant CF, Narayan RK, Grossman RG. Clinical evaluation of a miniature straingauge transducer for monitoring intracranial pressure. Neurosurgery. 1995;36(6):1137â€“40 (discussion 401).
Chambers IR, Siddique MS, Banister K, Mendelow AD. Clinical comparison of the Spiegelberg parenchymal transducer and ventricular fluid pressure. J Neurol Neurosurg Psychiatry. 2001;71(3):383â€“5.
Koskinen LO, Olivecrona M. Clinical experience with the intraparenchymal intracranial pressure monitoring Codman MicroSensor system. Neurosurgery. 2005;56(4):693â€“8 (discussion8).
Lescot T, Reina V, Le Manach Y, Boroli F, Chauvet D, Boch AL, et al. In vivo accuracy of two intraparenchymal intracranial pressure monitors. Intensive Care Med. 2011;37(5):875â€“9.
Zacchetti L, Magnoni S, Di Corte F, Zanier ER, Stocchetti N. Accuracy of intracranial pressure monitoring: systematic review and metaanalysis. Crit Care. 2015;19:420.
Acknowledgements
None
Funding
None.
Author information
Authors and Affiliations
Contributions
Conception and design: OD, OB, YZ, GR. Acquisition of data: YZ, OD. Analysis and interpretation of data: OB, OD, YZ, GR. Drafting the article: GR, OD. Critically revising the article: OB, YZ, OD. Statistical analysis: YZ. Administrative /technical/ material support: OB. Approved the submitted version: GR, OD, YZ, OB. All authors read and approved the final 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.
Appendix
Appendix
Equation A1 presents the ventricular steady state pressure:
where \({Q}_{CSF}\) is CSF flow, \({R}_{VEN}\) is ventricular resistance to flow and \({ICP}_{SAS}\) is subarachnoid pressure.
Ventricular resistance to flow (\({R}_{VEN}\)) presented as an ICPdependent nonlinear function:
where \({R}_{{VEN}_{0}}\) is a constant set to 250 [mmHg sec \({\text{ml}}^{  1}\)] and \({ICP}_{BR}\) is brain intraparenchymal pressure.
Values for each of the coefficients in the polynomial equations appearing in the Methods section are as follows:
Coefficient values for Eq.Â 1, subarachnoid space pressure:
Coefficient values for Eq.Â 2, brain parenchymal pressure (transmural):
Coefficient values for Eq.Â 3, ventricular pressure (transmural):
Coefficient values for Eq.Â 5, spinal thecal sac pressure:
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
Doron, O., Zadka, Y., Barnea, O. et al. Interactions of brain, blood, and CSF: a novel mathematical model of cerebral edema. Fluids Barriers CNS 18, 42 (2021). https://doi.org/10.1186/s1298702100274z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1298702100274z
Keywords
 ICP
 Cerebral edema
 Cerebral interstitial fluid
 Bulk flow
 CSF
 Outflow resistance
 Lumped parameter model
 Bloodâ€“brain barrier