- Research article
- Open Access
- Open Peer Review
An anatomy-based lumped parameter model of cerebrospinal venous circulation: can an extracranial anatomical change impact intracranial hemodynamics?
BMC Neurology volume 15, Article number: 95 (2015)
The relationship between extracranial venous system abnormalities and central nervous system disorders has been recently theorized. In this paper we delve into this hypothesis by modeling the venous drainage in brain and spinal column areas and simulating the intracranial flow changes due to extracranial morphological stenoses.
A lumped parameter model of the cerebro-spinal venous drainage was created based on anatomical knowledge and vessels diameters and lengths taken from literature. Each vein was modeled as a hydraulic resistance, calculated through Poiseuille’s law. The inputs of the model were arterial flow rates of the intracranial, vertebral and lumbar districts. The effects of the obstruction of the main venous outflows were simulated. A database comprising 112 Multiple Sclerosis patients (Male/Female = 42/70; median age ± standard deviation = 43.7 ± 10.5 years) was retrospectively analyzed.
The flow rate of the main veins estimated with the model was similar to the measures of 21 healthy controls (Male/Female = 10/11; mean age ± standard deviation = 31 ± 11 years), obtained with a 1.5 T Magnetic Resonance scanner. The intracranial reflux topography predicted with the model in cases of internal jugular vein diameter reduction was similar to those observed in the patients with internal jugular vein obstacles.
The proposed model can predict physiological and pathological behaviors with good fidelity. Despite the simplifications introduced in cerebrospinal venous circulation modeling, the key anatomical feature of the lumped parameter model allowed for a detailed analysis of the consequences of extracranial venous impairments on intracranial pressure and hemodynamics.
The relationship between extracranial venous system abnormalities and central nervous system (CNS) disorders has been subjected to scrutiny in the latest years [1, 2]. In recent studies, researchers investigated the association of venous abnormalities with aging , leukoaraiosis  and neurological disorders like transient global amnesia [5–7], transient monocular blindness [8–10], exertional headache , multiple sclerosis  (MS) and Alzheimer's disease [13, 14]. Various mechanisms have been proposed to explain the role of venous abnormalities in CNS pathologies without a definite and final consensus of the scientific community. For instance, it was hypothesized that internal jugular veins (IJV) abnormalities (valve incompetence, reflux or stenosis) may cause intracranial hypertension. This fact would lead to modifications in cerebrospinal fluid dynamics, cerebral blood flow and cerebral blood perfusion; however the exact mechanism underlying these effects remains unknown .
Recently, a pathological vascular condition termed chronic cerebrospinal venous insufficiency (CCSVI) was introduced . CCSVI was defined as the impairment of the main extracranial outflow routes, i.e., IJV, vertebral veins (VV), azygos vein (AZY). Hemodynamic and morphological criteria were proposed for its radiological diagnosis [12, 15]. The hypothesis assumed that CCSVI could have a causal relationship with MS , which has sparked significant debate in the literature [1, 16], as well as the criteria defining CCSVI itself and its diagnosis [15, 17].
Nevertheless, while even the concept of CCSVI itself remains highly contentious, the debate triggered interest in gaining a better understanding of the relationship between extracranial venous system abnormalities and CNS disorders . Recent studies aimed at defining a protocol to examine the extra- and intracranial venous system in vivo, describing guidelines, limits, benefits, sensitivity and specificity of various non-invasive and invasive imaging methods [15, 17, 18].
In addition, biomechanical concepts were used to characterize the outflow from the brain through the main extracranial veins in physiological conditions and in subjects with occluded IJVs [19–21]. A hydrodynamic model of the main intracranial fluid flows , and a patient specific model  were proposed in order to explain how an extracranial alteration could produce intracranial impairment. However, a mathematical tool for modeling the principal extracranial and intracranial veins with anatomical fidelity and the analysis of their hemodynamic changes after morphological alterations are missing. Regarding the blood fluid-dynamics per se, different modeling techniques are available and have been previously adopted in physiological and pathological conditions with a particular focus on the arterial circulation due to the high mortality rate associated with arterial diseases. The venous models described in literature are often a part of more general models designed for whole body descriptions [22, 23], or for extracranial  or intracranial compartment only [24–26]. They were used for analyzing the global system instead of the local hemodynamics. The impact of IJV morphological alterations on the main intracranial vessels flow was recently described with a 3D model of the IJVs flow coupled with a 1D model of the main intracranial vessels ; however the smallest intracranial veins and the other extracranial routes are missing. For the same purpose, a model recently developed by Zamboni and colleagues  for studying the role of the collateral neck veins in the brain drainage was integrated  with the validated mathematical model for the simulation of the intracranial circulation of Ursino et al. .
In this work we propose a novel lumped parameter model of the venous cerebrospinal system based on anatomical knowledge of the venous network. This kind of modelization borrows concepts from electric circuit design: flows are modeled as currents, vessel resistance to the flow as impedance and pressure differences as voltage differences. We applied such model to quantitatively assess the relationship between local modifications in vessels lumen and the overall changes in drainage capability of the venous cerebrospinal system, with particular focus on the intracranial area. For this work, we considered subjects in the supine position, where the principal extracranial outflow routes are the IJVs, and the flow rate in the VVs is lower compared to the upright position [19, 28–32].
The flows estimated with the model were compared with those obtained with the phase contrast (PC) magnetic resonance images (MRI) from 21 healthy controls (HC).
The intracranial reflux location obtained simulating IJV obstruction was compared with EchoColor Doppler (ECD) examinations of 112 MS patients.
Cerebrospinal venous lumped parameter model
Hydraulic lumped-parameter models allow for quantifying pressures and flow rates at different locations within a vessel network in which different tracts can be described as combinations of i) hydraulic resistances relating the flow rate in a conduit and the corresponding pressure drop, ii) compliance elements quantifying the tendency of the conduit to dilate due to inner pressure, and iii) inertance elements accounting for the inertia of the fluid in the conduit.
In our case, compliance and inertance elements were neglected due to the fact that simulations were carried out in steady state. We modeled each vessel by means of a hydraulic resistance R, related to the vessel’s dimensions and blood rheological properties through Poiseuille's law:
where L and D are the length and the equivalent diameter of the vessel and μ is blood viscosity, which was assumed equal to 3 cPoise .
Consistently, we assumed the standard linear relation between pressure drop (∆P) and flow rate (Q):
Hydraulic resistances can be connected in series or parallel circuits: in the former case the total resistance is the sum of all single components, in the latter the reciprocal of the total resistance is the sum of the reciprocal of all single components. More complex networks can be built by combining these two elemental configurations.
In order to model the main cerebrospinal veins, we considered 164 vessels and therefore 164 resistive elements, connected as shown in Fig. 1. The following venous vessels were included in the model, distinguishing between the right (subscript r) and left (subscript l) vessel whenever appropriate. For the intracranial area, ophthalmic veins (OVl, OVr), basal veins of Rosenthal (RVl, RVr), internal cerebral veins (ICVl, ICVr), inferior sagittal sinus (ISS), superior sagittal sinus (SSS), great vein of Galen (GV), straight sinus (SRS), transverse sinuses (TSl, TSr), posterior occipital sinuses (POS, POSl, POSr), basilar plexus (subdivided in right and left subportions, BPl, BPr), cavernous sinuses (CSl, CSr), superior petrosal sinuses (SPSl, SPSr), inferior petrosal sinuses (IPSl, IPSr), sigmoid sinuses (SSl, SSr). For the extracranial area, internal jugular veins (IJVl, IJVr), vertebral veins (divided into six segments, VVl 1…6, VVr 1…6), cervical plexus (divided into seven segments CP 1…7 and anterior and posterior segments, CPa, CPp), connective vessels between cervical plexus and vertebral veins (CPVVl 1…6, CPVVr 1…6), thoracic plexus (divided into twelve segments, TP 1…12), azygos vein (divided into twelve segments, AZ 1…12), connective vessels between thoracic plexus and azygos vein (TPAZ 1…12), inferior vena cava (divided into three segments, CV, CV1, CV2), lumbar plexus (divided into two segments, LP1, LP2), lumbar vein (divided into two segments, LV1, LV2) and connective vessels between lumbar plexus and lumbar veins (LPLV1, LPLV2). Facial venous circulation was excluded from the model, as well as external jugular veins, which drain the extra-cranium area. With regards the intra-cranium compartment, minor vessels such as venules were neglected.
Inputs of the system were arterial flows of cerebral (Qc), femoral (Qf) and vertebral (Qv1…Qv21) regions. With regards to the vertebral area, a single Qvi was considered for the i-th vertebra. Pressure was set equal to zero at two points of the hydraulic system: i) at the level of subclavian vein, based on its notably higher cross sectional diameter, and hence minimal hydraulic resistance, compared to the other vessels included in the model (Equation 1), and ii) at the level of the right atrium, consistent with negligible atrial pressure in physiological conditions. Based on the latter assumption, the hydraulic system was divided into two decoupled parts; the first representing the cerebral and upper spinal areas, and the second representing the cardiac area and the remaining part of the spinal column. This subdivision implied independent inputs at cervical (Qc and Qvi, i=1,…, 7) and thoracic/lumbar (Qf and Qvi, i=8,…, 21) levels.
We computed pressures (Pa) at each node and flow rates (ml/s) for each vessel by solving the linear system of equations associated with the network through a Matlab® script (The MathWorks Inc., Natick, MA, United States).
MRI acquisition and processing
Twenty-one HC (Male/Female = 10/11; mean age ± standard deviation = 31 ± 11 years) were recruited for measuring normal flow rates in the main cerebrospinal venous vessels with PC MRI. Data were acquired using a 1.5 T Siemens Magnetom Avanto scanner equipped with a 12-channel head coil and a 4-channel neck coil at IRCCS S. Maria Nascente of Fondazione Don Carlo Gnocchi ONLUS, Milan, Italy. The study was approved by the Don Gnocchi Foundation Ethics Committee (Milan, Italy) and written informed consent was obtained from all subjects prior to study entry.
The following sequences were acquired: 1) T1 weighted localizer of the brain and of the neck (three perpendicular slices) for the slice positioning of subsequent sequences; 2) dual-echo turbo spin echo (TR=2650 ms, TE=28/113 ms; echo train length=5; flip angle=150°, 50 interleaved, 2.5-mm-thick axial slices, with a matrix size=256×256, interpolated to 512×512, and a FOV=250×250 mm2) for the exclusion of vascular pathology, or white matter hyperintensities; 3) 2D time-of-flight magnetic resonance venography (TOF MRV) of the head and neck (TR=21 ms, TE=4.78 ms, 128 axial slices, FoV= 320×255 mm2, voxel of 1×0.5×3 mm3, distance factor of −20 %), with a caudal saturation band, in order to saturate the arterial signal, for positioning the following sequences; 4) 2D PC sequences for the flow quantification of cerebrospinal venous vessels, as follows. For the IJVs and the VVs, a PC slice was positioned between the fifth and the sixth cervical vertebrae (C5-C6), and between C2-C3, respectively. The IJVs were scanned with a slice positioned at the cervical level between C5 and C6, the VV were scanned with a slice at the cervical level between C2 and C3. Other vessels acquired, when visible, were: ISS, SRS, left and right RV, ICV, TS, SPS, SS, anterior and posterior portions of SSS. Their anatomical positions were identified on the proton density weighted images, or on TOF sagittal or coronal maximum intensity projection. The PC slices were positioned perpendicular to the axis of considered vessel. The maximum encoding velocity was selected depending on the expected maximum velocity of the examined vessel: 60 cm/s for IJVs and VV; 50 cm/s for SRS, TS and SSS; 20 cm/s for the ISS, ICV and R, and SPS. All the PC sequences were acquired with a retrospective cardiac gating through a finger pulse oxymeter, allowing the reconstruction of 25 time points over the cardiac cycle.
The PC data were processed with Argus software (Siemens Medical Solutions): for every vein and every time point, the vessel contour was manually drawn on the magnitude image and then transferred onto the corresponding phase image. The phase offset, i.e., the phase value of stationary structures, was estimated as the mean phase value inside a region of interest drawn in the muscular tissue. For every pixel inside the segmented vessels, and for every time point, the phase value, corrected by the offset, was mapped to velocity. From the area of the segmented vessel and the velocity in each pixel, the average flow rate (in ml/s) over the cardiac cycle was computed.
Comparison between flow rates computed with physiological model and MRI
The flow direction in each modeled vessel was compared with known physiology (Fig. 2).
The one-sample Wilcoxon signed rank test was used to assess whether the mean flow rate measured with PC MRI in HC differed significantly from the model estimate. For the VVs, we compared the flow rates measured with MRI with rates computed with the model at the second VV segment. A p-value <0.05 was deemed significant while <0.1 was considered a trend.
Venous obstruction simulation
The effect of the main outflow routes obstruction was simulated with the model. The following cases were implemented: monolateral and bilateral obstruction of the internal jugular veins, the vertebral veins and the proximal azygos. In details, our analysis was performed as following.
Obstruction of the IJV
In order to simulate partial or total IJV occlusion, the diameter of each IJV was reduced from 0 % (complete patency) to 100 % (complete block), with increments of 10 %. For each of the 121 simulations, we computed the flow rates of each intracranial vein and the pressure at each intracranial node and at IJV upstream. We also checked which kind of bilateral IJV obstruction produced intracranial retrograde flows.
Vertebral vein obstruction
The diameter of all the VV segments was reduced from 0 to 100 %, with increments of 10 %, to simulate their partial or total obstruction. For each configuration, we assessed pressures and flow rates of each vessel and when and where the reflux was produced.
The diameters of the first four proximal azygos segments were reduced from 0 to 100 % with increments of 10 %. For each configuration, we assessed the pressure of each azygos portion and checked if any azygos or thoracic plexus segment was refluent.
The ECD of 112 MS patients were acquired with the ultrasound system MyLabVinco (Esaote S.p.A., Firenze, Italy) by an expert vascular sonographer. Fourty-nine patients were examined as part of a previously published multicenter study  whereas the others were acquired during clinical examinations. An overview of the demographic and the clinical characteristics of all the MS patients are detailed in Table 1.
The ultrasound system was equipped with a linear array transducer probe (LA332, Esaote SpA, Italy: operating bandwidth 3 – 11 MHz; imaging frequencies: 3.5 – 5.0 – 6.6 – 10.0 MHz; Doppler frequencies: 3.3 – 5.0 MHz) for the IJV and VV examination and with a phased array transducer probe (PA240, Esaote SpA, Italy: operating bandwidth: 1–4 MHz; B-Mode frequencies: 2.0 – 2.5 – 3.3 MHz) for intracranial examination. A thick layer of US gel (Aquasonic 100 – Parker Laboratories Inc, Fairfield, New Jersey USA) was used to ensure a complete coupling between the transducer and the examined subject’s skin, to avoid black cones and dark areas on the US image.
During the IJV examination, the sonographer paid particular attention to prevent excessive pressure on the examined subject’s neck, in order not to change IJV shape and dimension. B-mode was used to assess flow obstruction due to stenosis, annulus, or septums.
The intracranial examination was performed through the temporal and transcondylar windows in order to assess the flow direction of RV and petrosal sinuses, respectively . The Multigate Quality Doppler Profiles was also used to detect flows with slow velocity or different directions in the same vessel and to detect intracranial reflux .
While the full protocol requires examining the subject in supine and sitting positions, we considered only the data obtained in the supine examination for this work.
The occurrence and topography of venous reflux observed through ECD examinations with the contemporaneous presence of mono- or bilateral IJV obstruction was compared with estimates from the model.
Comparison between model estimates and MRI measures
The model correctly reproduced local flow direction in all its segments (Fig. 2).
Figure 3 shows the comparison between the flow rate measured with PC MRI and the corresponding model estimate. The number of vessels correctly imaged and the p-value of the signed rank test are also reported. We did not measure SSr or SSl flow rates as it was not possible to correctly image these vessels solely in the perpendicular plane. The flow rate measured with MRI was not significantly different from the model estimate for the following veins: IJVr, IJVl, VVr, VVl, ICVr, RVr, RVl, SPSl, SPSr, SRS, posterior part of SSS, TSl, TSr. The flow rates of ISS, ICVl and anterior SSS estimated with the model were statistically different from those measured with MRI. The median flow rate obtained with MRI for anterior SSS (not shown in Fig. 3) and its 95 % CI were 2.5 [1.1–3.1] ml/s, which is statistically different (p<0.0001) from that estimated by the model (7.1 ml/s).
Venous obstruction simulation
The structure of the lumped parameter set allowed for separate analysis of the effects of upper and lower circulation anatomical anomalies. As the upper and lower parts of the model were decoupled, changes in the portion above the cardiac level did not affect the one below, and vice versa.
The intracranial pressure estimated by the model had a specific trend related with the degree of IJV diameter reduction as shown in Fig. 4a. The amplitude depends on the specific node. Figure 4a provides, as an example, the pressure at the confluence of RVs to GV. At this level, intracranial pressure doubled its baseline value when both the IJVs were obstructed, in particular when one of them had a diameter reduction of at least 80 % and the other one of at least 90 %. Fig. 4b shows how the SSS posterior end pressure increases with different IJVs constrictions, with respect to the baseline value, which was removed from each value.
The simulations carried out with IJV diameters of 1 cm  and 0.7 cm  showed intracranial pressure increase for lower IJV diameter reductions. At the level of the RVs confluence to GV, pressure doubled its baseline value for bilateral IJVs obstruction of 70 % in the latter and 60 % in the former (Additional files 2 and 3) case.
In regards to the IJV upstream pressure, it doubled its initial value for IJV diameter reduction of at least 20 %.
When the total occlusion of one IJV was simulated, retrograde flow in the corresponding SS was computed whenever a partial diameter reduction of the other IJV was modeled.
The simulation of bilateral IJV occlusion greater than 90 % produced retrograde flow in the CS, IPS (Fig. 5b) and SPS. The simulations carried out with IJV diameters of 1 cm  and 0.7 cm  showed reflux of intracranial vessels for lower percentage of bilateral IJV diameter reduction (80 % in the former and of 70 % in the latter case).
The results obtained with IJV collaterals differed from those obtained without collaterals by less than 3.5 %. The reflux was observed in the same vessels as without collaterals.
Vertebral vein obstruction
The simulation of VV diameter reduction showed an increase in global high pressure in the cervical and vertebral areas, with higher values for severe VVs obstruction. Conversely, the intracranial area was less affected, where the pressure increases by less than 14 %. The CPa and CPp upstream, respectively the confluences of the BP and the OS to the cervical and vertebral areas, were highly affected by VV obstruction.
Retrograde flow in POS, POSl and POSr was observed for bilateral VV occlusion higher than 60 %.
The obstruction of the proximal azygos segments caused a pressure increment along the azygos. The pressure progressively increased as the azygos diameter was reduced. The obstruction effect was higher for the nodes near to the obstruction itself and lower for the distal parts. We reported (Fig. 6a) the pressure between AZY4 and AZY5, and AZY5 and AZY6 as representative of the former, and the one between AZY11 and AZY12 for the latter.
Retrograde flows were computed in the azygos itself and in the thoracic plexus: the higher the degree of diameter reduction, the higher the number of distal azygos and thoracic plexus segments affected by the reflux. Figure 6b shows the flow rate trend with progressive azygos diameter reduction for AZY1, AZY5, and AZY11, as representative of segments from proximal to distal azygos.
EchoColor Doppler results
ECD data of MS patients showed IJV abnormalities caused by stenosis, membranous obstruction or hypoplasia in 86 % of cases. Among the stenotic or hypoplastic cases, unilateral IJV occlusion was detected in 28 patients (21 of them had unilateral IJV but not VV occlusion), while bilateral in 37 (17 without VV occlusions). Intracranial reflux was observed in 21/28 subjects with monolateral IJV occlusion (17/21 considering those without VV occlusions), in 28/37 subjects with bilateral IJV occlusion (15/17 considering those without VV occlusions). The reflux was always located in CS, IPS and SPS, but never located in the RV.
In the present work we described a lumped parameter model of the cerebrospinal venous network based on a realistic anatomical paradigm, aimed at determining if an extracranial morphological change can cause an intracranial venous flow alteration.
Compared to previously published models of the physiologic and impaired cerebrospinal venous circulation [19, 21, 23, 25, 26], our model was designed with a detailed correspondence between its elements and anatomy. The accuracy of the model was evaluated comparing the flow estimates with that obtained with PC-MRI from a group of healthy subjects.
We applied the model to mimic obstruction of the main extracranial veins and to estimate the changes in intracranial pressure and flow. We compared the location of flow inversion simulated with IJV obstruction and what we observed with ECD in a group of 112 MS patients.
Comparison between model estimates and MRI measures
The physiologic flow rates computed by the model were realistic in terms of direction (Fig. 2) and amplitudes (Fig. 3). These agreements were not obvious since we a priori set the arterial input flow directions; the network connections; the venous diameters and lengths and the zero-level pressure at two points of the hydraulic system. We did not impose the venous flow directions or pressures, with the exception of the pressures in subclavian vein and the right atrium. In regards to the flow amplitude, the mean flow rates estimated with MRI were not significantly different from those predicted by the model, with the exception of ISS and ICVl.
The best correspondence between predicted and measured flow rates were found for those vessels whose anatomy can be approximated by an almost straight pipe of constant cross sectional area, i.e., where the basic assumption for Poiseuille's law holds, and whose dimension, shape and blood velocity allowed for the acquisition of good quality PC MRI data.
The significant difference between ISS (p=0.004) and ICVl (p=0.01) flow rates and those predicted by the model could be due to the small dimensions and low flow velocity of all these vessels and to the curve shape of the ISS. The same reasons could explain the trend for difference found between ICVr and both the SPS (p=0.06).
Regarding the SSS, its cross sectional dimensions increased anteriorly to posteriorly; as did its flow rate, consistently with the presence of the superior cerebral veins converging into it. We modeled the SSS as a single resistance and did not model the superior cerebral veins due to their high inter- subject variability in terms of number and size. For these matters the flow rate measured with the anterior, but not posterior, acquisition was significantly different from the model estimate. Furthermore, the SSS curve shape could explain the trend for difference between the model estimate and the posterior MRI measurement (p=0.055). The mean SSS flow rate measured with the posterior acquisition was similar to that reported by Gideon et al. . (443±124 ml/min, p=0.1) and Mehta et al.  (388±169 ml/min, p=0.9) on a group of 14 and 15 healthy volunteers, respectively. The model estimate was similar to the cited works. The same reasons could explain why there was a trend for the TS (p=0.09 and p=0.1for the left and the right side, respectively) for slower flow rates measured with MRI compared to that predicted by the model.
Another general source of mismatch between our model estimates and the MRI measurements involves the coupled vessels. We assumed a right-left symmetry of the venous network when designing the model, while in vivo asymmetry is described with a wide range of anatomic patterns [45–49]. The IJV, VV, and TS, seem to have right dominance regarding their size and flow [47, 49].
As a matter of principle, different anatomical connections could be additionally implemented with the same approach, even though this would require changing the equations associated to the network. The same model structure can instead be used for a variety of analyses, such as the study of asymmetric vessels, only by changing the dimensions of the vessels.
Venous obstruction simulation
The simulations obtained using different combinations of IJV diameters showed that the intracranial pressure is stable for a wide range of IJVs diameters and increases for severe IJV stenosis. Our finding is similar to that showed by the recent model of Gadda , where a mild effect on the venous sinuses pressure was produced by halving the IJVs conductances, and a strong effect by setting the IJVs conductances close to zero. The previously cited model analyzed the venous sinuses pressure as a single node of the network. Conversely, our model includes the main intracranial vessels, so we could simulate the pressures at their ends for progressive bilateral IJV constrictions. Even with different amplitudes, the relationship between venous pressure and IJV diameter reduction had the same shape for the different intracranial veins. The high pressures estimated at the confluence of RVs and ICVs to GV due to severe IJVs diameter reduction is a theoretical biomechanical finding supporting the hypothesis and the experimental findings of previous works [20, 50–52]. It has been suggested that such a high pressure can damage the brain tissue, with possible blood–brain-barrier rupture [50, 51] or permeability increase [13, 53]. In particular, since the Galenic system drains the periventricular white matter and the basal ganglia, its increase of pressure can damage those brain areas [51, 53]. The intracranial vein hypertension, especially that of the SSS (Fig. 4b), can also impact the cerebrospinal fluid flow hemodynamics. Indeed, the cerebrospinal fluid flow absorption is possible if cerebrospinal fluid pressure exceeds that of the venous sinuses by about 5–7 mmHg . For severe degrees of IJV stenosis, our model predicted a pressure increase in the SSS similar or even higher than the previously reported pressure gradient. A relationship between IJV stenosis, reflux and cerebrospinal fluid flow alterations was recently observed in different CNS disorders, like MS and AD [13, 14, 20, 54, 55]. The cerebrospinal fluid flow alterations were also recently reported to correlate with “dirty appearing white matter” in healthy subjects . The brain edema seen in cases of venous sinus thrombosis or stenosis could also be due to the decrease in cerebrospinal fluid drainage caused by an increase in venous pressure . In a recent review , Beggs explained the link between IJV stenosis, its hypertension and the impairment of cerebrospinal fluid flow absorption into the venous blood.
Other consequences of severe neck vein stenoses predicted by our model are the reduced intracranial venous flow rate, and even flow inversion. Schelling in 1986  suggested that refluent IJVs or intracranial veins could produce retrograde venous hypertension, i.e., from downstream to upstream the venous network. The resulting venous overload, engorgement, dilatation and leakage were proposed to be possible causes of MS lesions. The location and shape of such lesions were also suggested to be similar to those found through Schlesinger’s experiments with gelatin solution forced in the GV of human cadavers.
With the modelization used in this work, blood flow in any vein is determined by the pressure drop between its extremities, being directed toward the end with the lowest pressure. If the pressure difference between the vein ends reverts, the flow becomes retrograde. This general mechanism explains the cases of reflux found through our simulations, as explained for the intracranial, vertebral and azygos areas hereafter.
The upstream of one IJV is more influenced by its diameter change than intracranial vessels: its pressure begins to increase for small IJV diameters reduction. When both the IJVs are stenotic, our model predicted retrograde flow for the intracranial vein proximal to the IJVs themselves, in particular for the CS, IPS and SPS. Figure 5 represents how the flow inversion occurs in the IPSl. The pressures of its two ends were plotted over progressive reductions in IJV diameter (Fig. 5a): since they are superimposed on each other, it is clearly visible when one exceeds the other, i.e., for severe bilateral IJVs stenosis. Together with this pressure drop inversion, the IPSl flow inverted its direction (Fig. 5b). The preferred drainage path of CS, IPS and SPS in this case was the anterior cervical plexus (Fig. 7) instead of the IJVs (Fig. 2), since it offers lower resistance to flow. The presence of collateral vessels, as modeled in this work, did not prevent intracranial pressure increase and reflux. Indeed, one or more collateral vessels in parallel to a stenotic IJV have higher resistance to flow as compared to a patent IJV, even if the total cross sectional area is the same . If a clearer picture was needed, the model could be tuned for a single patient-specific case by setting ad hoc geometry, number, position and connections of collaterals. By doing this, it could be tested how a wider or higher number of collateral vessels would guarantee enough drainage and lower intracranial pressures, avoiding intracranial vessels reflux. Future work involving this kind of simulation could provide a better understanding about how collateral flow impacts intracranial pressure. Indeed, collateral number, dimensions and their anastomosis to the IJV are variable and are fundamental as an alternative pathway in clinical cases of IJV obstruction [27, 28].
The VV morphological obstructions had a low impact on the intracranial region, but a high impact on the pressure of the cervical and vertebral areas, besides on the areas connecting them to the intracranial vessels. When the pressure increase was higher in the upper portion of the VVs than in the intracranial region, a retrograde flow could occur in the posterior occipital sinuses. The blood of the last vessels was redirected to the IJVs instead of flowing to the CPp and VVs.
For a patent azygos, the pressure of its distal tracts is higher than that of the proximal ones: this is visible at Fig. 6a, for diameter reduction = 0 %. The flow is directed from its distal to its proximal tracts, to the superior vena cava, as shown in Fig. 2. The simulation of the proximal azygos obstruction produced asevere pressure increase (Fig. 6a, pressure between AZY4 and AZY5 and between AZY5 and AZY6). Conversely, the pressures of the distal districts were almost unaffected by the change (Fig. 6a, pressure between AZY11 and AZY12). When the obstruction caused an excess of pressure in the proximal compared to the distal part of an azygos portion, the pressure gradient changed its direction. The blood flow reverted its physiologic direction according to the pressure gradient. This can be observed in Fig. 6b, which shows how the azygos flow rate changes when proximal azygos obstruction is simulated. The AZY5 flow becomes retrograde when its distal pressure (between AZY5 and AZY6) exceeds the proximal one (between AZY4 and AZY5). For this particular case, this happens when the proximal azygos diameter is halved. The same can be observed for AZY11 flow, with less diameter reduction. A similar mechanism affected the thoracic plexus. The thoracic and azygos refluxes were drained to the right atrium through the inferior vena cava instead of the superior vena cava.
Jugular vein obstructions: comparison of model simulations and real cases
Comparing our ECD data with the model estimates in the most prevalent clinical situation (the IJV obstruction), we found agreements in intracranial reflux topography. Comparing our ECD data with the model estimates in the most prevalent clinical situation (the IJV obstruction), we found agreements in intracranial reflux topography.
A quantitative comparison could not be done since we did not measure the IJV diameter during ECD acquisitions. However, the results of the model simulations have to be read as a proof of concept: the percentages of IJVs diameter reductions corresponding to intracranial reflux, or to intracranial pressure increase of at least 200 %, are not intended to be considered absolute thresholds. Indeed, they depend on the chosen “normal” IJVs diameter and other vessel dimensions. Regarding venous dimensions, a high inter-subject variability is described in literature, with different physiological vessels dimensions and flow rates among the studies, and a high deviation from the mean within each study [30, 32, 44, 49, 58].
Furthermore, the venous reflux detected in clinical cases often involves only a portion of the vessel and can be intermittent . This fact can be explained by observing the occurrence of separated flow regions in a single vessel, as studied with a Lagrangian approach to fluid flow analysis [59, 60]. Conversely, with the present kind of modelization, a vessel can be classified as refluent if its flow rate is completely reverted. Local flow separation cannot be described through our approach. However, considering the relatively low values of Reynolds' number (maximum = 715 in the SSS), we may reasonably infer that such peculiar flow structures wouldn't significantly affect the vessel resistance calculations in our model.
Limits and perspectives
The results of our model and all the simulations provided by it have to be evaluated while taking into consideration the main limitations of our type of modelization: the model is steady-state and it is valid only for analyzing the supine position.
Regarding the subject position, a future work for studying a standing subject should consider the effect of atmospheric pressure on IJVs, which collapse, and the role of collateral pathways .
As to the steady-state limitation, this means that every pressure and flow of the network is not time-variant, so the modulation of the venous velocity due to the cardiac or the thoracic pump are neglected. A recent work  showed that the latter highly influences the venous return of the IJVs, especially for in the deep inspiration, because of the aspiration effect of the pleural cavity. However, the mean flow of the intracranial veins is less modulated by the respiratory phases. The minor effect of respiration and cardiac pulsations on intracranial vein flow was also shown through MRI with the analysis of the frequency components in the range of respiration and heart beat. The flattened flow rate waveform of cerebral sinus veins can be due to their fixed diameter and protective effect of the dura mater . Moreover, the IJV cross sectional area changes over the cardiac cycle . However, due to the high inter-subject variability in terms of IJV cross sectional areas and ranges of variation, the choice of considering an average constant value can be considered a reasonable approximation.
Another aspect to point out is that the refluxes observed in clinical cases could be due to the sum of the effects of many morphological alterations (i.e., other abnormalities besides IJV or VV or AZY obstruction), while we separately simulated the effects of the main routes obstruction, obtaining intracranial refluxes only for severe IJV stenoses. Nonetheless, it is possible to use this model to test the effect of many other morphological alterations. Owing to the linearity of the model, the effects of any combination of anatomical changes may be obtained through the superposition principle, i.e., by linearly combining the effects of single anatomical alterations. It is hence straightforward to adapt the model to patient-specific features.
In the context of a wide scientific interest for the link between venous flow and CNS disorders, the use of this model for the simulation of specific patterns (i.e., pervious vs stenotic vessels) or venous variants provides the estimate of their impact on pressures and flow in the main intracranial vessels.
Zivadinov R, Chung CP. Potential involvement of the extracranial venous system in central nervous system disorders and aging. BMC Med. 2013;11:260-7015-11-260.
D'haeseleer M, Cambron M, Vanopdenbosch L, De Keyser J. Vascular aspects of multiple sclerosis. Lancet Neurol. 2011;10:657–66.
Chung CP, Wang PN, Wu YH, Tsao YC, Sheng WY, Lin KN, et al. More severe white matter changes in the elderly with jugular venous reflux. Ann Neurol. 2011;69:553–9.
Chung CP, Hu HH. Pathogenesis of leukoaraiosis: role of jugular venous reflux. Med Hypotheses. 2010;75:85–90.
Chung CP, Hsu HY, Chao AC, Sheng WY, Soong BW, Hu HH. Transient global amnesia: cerebral venous outflow impairment-insight from the abnormal flow patterns of the internal jugular vein. Ultrasound Med Biol. 2007;33:1727–35.
Lochner P, Nedelmann M, Kaps M, Stolz E. Jugular valve incompetence in transient global amnesia. A problem revisited. J Neuroimaging. 2014;24:479–83.
Cejas C, Cisneros LF, Lagos R, Zuk C, Ameriso SF. Internal jugular vein valve incompetence is highly prevalent in transient global amnesia. Stroke. 2010;41:67–71.
Cheng CY, Chang FC, Chao AC, Chung CP, Hu HH. Internal jugular venous abnormalities in transient monocular blindness. BMC Neurol. 2013;13:94-2377-13-94.
Chung CP, Hsu HY, Chao AC, Chang FC, Sheng WY, Hu HH. Detection of intracranial venous reflux in patients of transient global amnesia. Neurology. 2006;66:1873–7.
Hsu HY, Chao AC, Chen YY, Yang FY, Chung CP, Sheng WY, et al. Reflux of jugular and retrobulbar venous flow in transient monocular blindness. Ann Neurol. 2008;63:247–53.
Doepp F, Valdueza JM, Schreiber SJ. Incompetence of internal jugular valve in patients with primary exertional headache: a risk factor? Cephalalgia. 2008;28:182–5.
Zamboni P, Galeotti R, Menegatti E, Malagoni AM, Tacconi G, Dall'Ara S, et al. Chronic cerebrospinal venous insufficiency in patients with multiple sclerosis. J Neurol Neurosurg Psychiatry. 2009;80:392–9.
Beggs C, Chung CP, Bergsland N, Wang PN, Shepherd S, Cheng CY, et al. Jugular venous reflux and brain parenchyma volumes in elderly patients with mild cognitive impairment and Alzheimer's disease. BMC Neurol. 2013;13:157-2377-13-157.
Chung CP, Beggs C, Wang PN, Bergsland N, Shepherd S, Cheng CY, et al. Jugular venous reflux and white matter abnormalities in Alzheimer's disease: a pilot study. J Alzheimers Dis. 2014;39:601–9.
Zivadinov R, Bastianello S, Dake M, Ferral H, Haacke EM, Haskal ZJ, et al. Recommendations for multimodal noninvasive and invasive screening for detection of extracranial venous abnormalities indicative of chronic cerebrospinal venous insufficiency: a position statement of the International Society for Neurovascular Disease. J Vasc Interv Radiol. 2014;25:1785–94.
Comi G, Battaglia MA, Bertolotto A, Del Sette M, Ghezzi A, Malferrari G, et al. Observational case–control study of the prevalence of chronic cerebrospinal venous insufficiency in multiple sclerosis: results from the CoSMo study. Mult Scler. 2013;19:1508–17.
Dolic K, Marr K, Valnarov V, Dwyer MG, Carl E, Hagemeier J, et al. Sensitivity and specificity for screening of chronic cerebrospinal venous insufficiency using a multimodal non-invasive imaging approach in patients with multiple sclerosis. Funct Neurol. 2011;26:205–14.
Dolic K, Siddiqui AH, Karmon Y, Marr K, Zivadinov R. The role of noninvasive and invasive diagnostic imaging techniques for detection of extra-cranial venous system anomalies and developmental variants. BMC Med. 2013;11:155-7015-11-155.
Zaniewski M, Simka M. Biophysics of venous return from the brain from the perspective of the pathophysiology of chronic cerebrospinal venous insufficiency. Rev Recent Clin Trials. 2012;7:88–92.
Beggs CB. Venous hemodynamics in neurological disorders: an analytical review with hydrodynamic analysis. BMC Med. 2013;11:142-7015-11-142.
Caiazzo A, Montecinos G, Müller L, Haacke EM, Toro E. Computational haemodynamics in stenotic internal jugular veins. J Math Biol. 2015;70(4):745–72.
Lakin WD, Stevens SA, Tranmer BI, Penar PL. A whole-body mathematical model for intracranial pressure dynamics. J Math Biol. 2003;46:347–83.
Mynard JP, Smolich JJ. One-dimensional haemodynamic modeling and wave dynamics in the entire adult circulation. Ann Biomed Eng. 2015;43(6):1443–60.
Ursino M, Lodi CA. A simple mathematical model of the interaction between intracranial pressure and cerebral hemodynamics. J Appl Physiol (1985). 1997;82:1256–69.
Ambarki K, Baledent O, Kongolo G, Bouzerar R, Fall S, Meyer ME. A new lumped-parameter model of cerebrospinal hydrodynamics during the cardiac cycle in healthy volunteers. IEEE Trans Biomed Eng. 2007;54:483–91.
Bergsneider M, Alwan AA, Falkson L, Rubinstein EH. The relationship of pulsatile cerebrospinal fluid flow to cerebral blood flow and intracranial pressure: a new theoretical model. Acta Neurochir Suppl. 1998;71:266–8.
Zamboni P, Sisini F, Menegatti E, Taibi A, Malagoni AM, Morovic S, et al. An ultrasound model to calculate the brain blood outflow through collateral vessels: a pilot study. BMC Neurol. 2013;13:81-2377-13-81.
Gadda G, Taibi A, Sisini F, Gambaccini M, Zamboni P, Ursino M. A new hemodynamic model for the study of cerebral venous outflow. Am J Physiol Heart Circ Physiol. 2015;308:217–31.
Gisolf J, van Lieshout JJ, van Heusden K, Pott F, Stok WJ, Karemaker JM. Human cerebral venous outflow pathway depends on posture and central venous pressure. J Physiol. 2004;560:317–27.
Valdueza JM, Schmierer K, Mehraein S, Einhaupl KM. Assessment of normal flow velocity in basal cerebral veins. A transcranial doppler ultrasound study. Stroke. 1996;27:1221–5.
Schaller B. Physiology of cerebral venous blood flow: from experimental data in animals to normal function in humans. Brain Res Brain Res Rev. 2004;46:243–60.
Ciuti G, Righi D, Forzoni L, Fabbri A, Pignone AM. Differences between internal jugular vein and vertebral vein flow examined in real time with the use of multigate ultrasound color Doppler. AJNR Am J Neuroradiol. 2013;34:2000–4.
Cho YI, Kensey KR. Effects of the non-Newtonian viscosity of blood on flows in a diseased arterial vessel. Part 1: Steady flows. Biorheology. 1991;28:241–62.
Newton DG, Potts NT. Radiology of the Skull and Brain. 1st ed. St Louis: C.V. Mosby; 1974.
Lambertini G. Anatomia umana. 1st ed. Padova (Italy): Piccin; 1978.
Rhoton Jr AL. The cerebral veins. Neurosurgery. 2002;51:159–205.
Felten DL, Shetty A. Atlante di neuroscienze di Netter. Milano (Italy): Elsevier srl; 2011.
Zamboni P, Consorti G, Galeotti R, Gianesini S, Menegatti E, Tacconi G, et al. Venous collateral circulation of the extracranial cerebrospinal outflow routes. Curr Neurovasc Res. 2009;6:204–12.
Doepp F, Hoffmann O, Schreiber S, Lammert I, Einhaupl KM, Valdueza JM. Venous collateral blood flow assessed by Doppler ultrasound after unilateral radical neck dissection. Ann Otol Rhinol Laryngol. 2001;110:1055–8.
Bastianello S, Romani A, Viselner G, Tibaldi EC, Giugni E, Altieri M, et al. Chronic cerebrospinal venous insufficiency in multiple sclerosis: clinical correlates from a multicentre study. BMC Neurol. 2011;11:132-2377-11-132.
Laganà MM, Preti MG, Forzoni L, D'Onofrio S, De Beni S, Barberio A, et al. Transcranial ultrasound and magnetic resonance image fusion with virtual navigator. IEEE Trans Multimedia. 2013;15:1039–48.
Forzoni L, D'Onofrio S, Farina M, Semplici P, Corsi M, Furia R, et al. Combined Ultrasound Technologies and Optimized Probe Design for Neck Veins Examination. Proceedings of the IASTED Biomedical Engineering BioMed Congress; 2012. p. 252–9.
Gideon P, Thomsen C, Gjerris F, Sorensen PS, Stahlberg F, Henriksen O. Measurement of blood flow in the superior sagittal sinus in healthy volunteers, and in patients with normal pressure hydrocephalus and idiopathic intracranial hypertension with phase-contrast cine MR imaging. Acta Radiol. 1996;37:171–6.
Mehta NR, Jones L, Kraut MA, Melhem ER. Physiologic variations in dural venous sinus flow on phase-contrast MR imaging. AJR Am J Roentgenol. 2000;175:221–5.
Park HK, Bae HG, Choi SK, Chang JC, Cho SJ, Byun BJ, et al. Morphological study of sinus flow in the confluence of sinuses. Clin Anat. 2008;21:294–300.
Saiki K, Tsurumoto T, Okamoto K, Wakebe T. Relation between bilateral differences in internal jugular vein caliber and flow patterns of dural venous sinuses. Anat Sci Int. 2013;88:141–50.
Manara R, Mardari R, Ermani M, Severino M, Santelli L, Carollo C. Transverse dural sinuses: incidence of anatomical variants and flow artefacts with 2D time-of-flight MR venography at 1 Tesla. Radiol Med. 2010;115:326–38.
Alper F, Kantarci M, Dane S, Gumustekin K, Onbas O, Durur I. Importance of anatomical asymmetries of transverse sinuses: an MR venographic study. Cerebrovasc Dis. 2004;18:236–9.
Stoquart-Elsankari S, Lehmann P, Villette A, Czosnyka M, Meyer ME, Deramond H, et al. A phase-contrast MRI study of physiologic cerebral venous flow. J Cereb Blood Flow Metab. 2009;29:1208–15.
Talbert DG. Raised venous pressure as a factor in multiple sclerosis. Med Hypotheses. 2008;70:1112–7.
Schlesinger B. The venous drainage of the brain, with special reference to the galenic system. Brain. 1939;62:274–91.
Schelling F. Damaging venous reflux into the skull or spine: relevance to multiple sclerosis. Med Hypotheses. 1986;21:141–8.
Wu X, Wu W, Wu J, Zhang HL. Diagnosis: thrombosis of the vein of Galen. Ann Saudi Med. 2014;34:450.
Zivadinov R, Magnano C, Galeotti R, Schirda C, Menegatti E, Weinstock-Guttman B, et al. Changes of cine cerebrospinal fluid dynamics in patients with multiple sclerosis treated with percutaneous transluminal angioplasty: a case–control study. J Vasc Interv Radiol. 2013;24:829–38.
Magnano C, Schirda C, Weinstock-Guttman B, Wack DS, Lindzen E, Hojnacki D, et al. Cine cerebrospinal fluid imaging in multiple sclerosis. J Magn Reson Imaging. 2012;36:825–34.
Beggs CB, Magnano C, Shepherd SJ, Belov P, Ramasamy D, Hagemeier J, et al. Dirty-appearing white matter in multiple sclerosis: volumetric MR imaging and magnetization transfer ratio histogram analysis. J Neuroimaging. 2015;24:1935–40.
Iencean SM, Poeata I, Iencean AS, Tascu A. Cerebral venous etiology of intracranial hypertension and differentiation from idiopathic intracranial hypertension. PubMed - NCBI. 2015;36:825–34.
Beards SC, Yule S, Kassner A, Jackson A. Anatomical variation of cerebral venous drainage: the theoretical effect on jugular bulb blood samples. Anaesthesia. 1998;53:627–33.
Simka M. Chronic cerebrospinal venous insufficiency: current perspectives. J Vasc Diagn. 2014;2:1–13.
Shadden SC, Taylor CA. Characterization of coherent structures in the cardiovascular system. Ann Biomed Eng. 2008;36:1152–62.
Zamboni P, Menegatti E, Pomidori L, Morovic S, Taibi A, Malagoni AM, et al. Does thoracic pump influence the cerebral venous return? J Appl Physiol. 2012;112:904–10.
Sisini F, Tessari M, Gadda G, Di Domenico G, Taibi A, Menegatti E, et al. An ultrasonographic technique to assess the jugular venous pulse: a proof of concept. Ultrasound Med Biol. 2015;41:1334–41.
The authors are grateful to Dr. Marco Petilli for helping in the figures preparation and to Dr. Paul Kokeny and Dr. Sean Sethi for reviewing the manuscript
The authors declare that they have no competing interests.
SM, LM, PC, AR, GF, and MML participated in the study design, projected and created the model. PC, AB, SM, LM, SV, and MML acquired and processed the echodoppler and MRI. SM, LM, MML performed the simulations with the model and the statistical analysis. SM, LM, PC, AR, GF, EV, SV and MML compared the model estimates with the MRI and echodoppler data, discussed the results. All authors gave substantial contribution or helped to draft the manuscript; read and approved the final manuscript.
Stefania Marcotti and Lara Marchetti contributed equally to this work.
Lengths and diameters of venous vessels included in the model. The plexus was considered an equivalent measure, due to the fact these districts are composed by a network-like structure. *: these vessels were subsequently divided in segments, as explained in the text. **: for IJVs also a physiological diameter of 1 cm  and 0.73 cm  were considered, to evaluate intracranial pressure and reflux onset related to smaller IJVs reported in literature.
Results obtained with model simulations of internal jugular veins (IJVs) diameter reduction and physiologic IJV diameter reported in Gisolf et al. . A. Pressure at the confluence of Rosenthal veins (RVs) to Galen vein (GV) vs internal jugular veins (IJVs) diameter reduction. Pressure doubled its initial value when both the IJV diameters reduced at least of 70 %. The result is symmetric for correspondent diameter reduction rates on the opposite side (not shown). B. Bars showing IJVl upstream pressure increase versus IJVs diameter reduction. Pressure value doubled its initial value when IJVl underwent a diameter reduction of at least 20 %.
Results obtained with model simulations of internal jugular veins (IJVs) diameter reduction and physiologic IJV diameter reported in Ciuti et al. . A. Pressure at the confluence of Rosenthal veins (RVs) to Galen vein (GV) vs internal jugular veins (IJVs) diameter reduction. Pressure doubled its initial value when both the IJV diameters reduced at least of 60 %. The result is symmetric for correspondent diameter reduction rates on the opposite side (not shown). B. Bars showing IJVl upstream pressure increase versus IJVs diameter reduction. Pressure value doubled its initial value when IJVl underwent a diameter reduction of at least 20 %.
About this article
Cite this article
Marcotti, S., Marchetti, L., Cecconi, P. et al. An anatomy-based lumped parameter model of cerebrospinal venous circulation: can an extracranial anatomical change impact intracranial hemodynamics?. BMC Neurol 15, 95 (2015) doi:10.1186/s12883-015-0352-y
- Blood flow
- Blood pressure
- Cerebrovascular circulation
- Nervous system