Skip to main content

An ultrasound model to calculate the brain blood outflow through collateral vessels: a pilot study



The quantification of the flow returning from the head through the cervical veins and the collaterals of the internal jugular vein (IJV), is becoming of prominent interest in clinical practice. We developed a novel model to calculate the cerebral venous return, normalized to the arterial inflow, in the different segments of the IJV.


We assessed, by established Echo Colour Doppler (ECD) methodology, the head inflow (HBinF) defined as the sum of common carotids and vertebral arteries, as well as the cerebral flow (CBF) defined as the sum of internal carotid and vertebral arteries. We also assessed the head outflow (HBoutF) defined as the sum of the measurements at the junction of the IJV and the vertebral veins. In addition, we also calculated the collateral flow index (CFI) by estimating the flow which re-enters directly into the superior vena cava as the amount of blood extrapolated by the difference between the HBinF and the HBoutF. We preliminarily tested the model by comparing ten healthy controls (HC) with ten patients affected by chronic cerebral spinal venous insufficiency (CCSVI), a condition characterized by some blockages in the IJV which are bypassed by collateral circulation.


In HC the HBinF was 956+-105ml/min, whereas the HBoutF was > 90% of the HBinF, leading to a final CFI value of 1%. The last result shows that a very small amount of blood is drained by the collaterals. In upright we confirmed a reduction of the outflow through the IJV which increased CFI to 9%. When we applied the model to CCSVI, the HBinF was not significantly different from controls. In supine, the flow of CCSVI patients in the IJV junction was significantly lower (p < 0.001) while the correspondent CFI value significantly increased (61%, p < 0.0002).


Our preliminary application of the novel model in the clinical setting suggests the pivotal role of the collateral network in draining the blood into the superior vena cava under CCSVI condition.

Peer Review reports


There is general agreement in considering the internal jugular veins (IJVs) as the major route of cerebral outflow in the supine position, and the vertebral veins (VVs) as the major route of brain drainage in upright [14]. In a recent consensus, the IJV was subdivided into 3 segments: the segment J3 or higher, which is anatomically located at the carotid bifurcation and the mandibular angle; the middle segment or J2, related to the ipsilateral thyroid lobe; finally, the lower end or J1, corresponding to the confluence with the brachio-cephalic vein trunk [5].

In a recent paper we have shown that the flow tends physiologically to grow in volume from J3 to J1, both in basal conditions and under standardized conditions of activation of the thoracic pump [6]. The main question to be answered is why the IJV flow is increased from the skull to the chest. Our hypothesis is to consider the possibility of blood re-entry from jugular collaterals into the main trunk (in the following we will use the term re-entry to indicate when a collateral channel is flow tributary of the major truncal pathway). Furthermore, the increase of IJV flow along the extra-cranial segment could be related to the re-entry volume through collateral vessels draining not only the blood of face and neck soft tissue, but also a rate coming from the brain through extra-intra-cranial anastomosis. In fact, there is a never assessed quota of the head inflow that is conveyed into the IJV more caudally with respect to the J3 position, through intra- and extra-cranial anastomosis. We are aware of anatomical presence of intra- and extra-cranial connection [7] but their physiological contribution to brain circulation is completely unknown. To this aim we have developed an haemodynamic model which describes quantitatively the neck pathway of the cerebral venous return, normalized with respect to the arterial inflow. Flow parameters have been measured by means of established echo-colour Doppler (ECD) methodology. Finally, in the second phase of the research, we have performed preliminary measurements in normal subjects and in patients affected by Chronic Cerebrospinal Venous Insufficiency (CCSVI) [5].


First phase of the study

Total of eleven healthy volunteers were screened for CCSVI absence by means of established ECD criteria [5]. One out of eleven screened subjects presented with >2 ECD criteria positive for CCSVI and, consequently, was excluded from the healthy controls. Ten out of eleven entered the study (age ranging from 23 to 42 y.o., male:female ratio 3:2). This prospective study was in accordance with Ethical Standards of the Committee on Human Experimentation of the University of Ferrara. All the study participants were non-invasively investigated by means of ultrasonic scanning with an ECD machine (ESAOTE My-Lab 70, Genoa, Italy) at the same condition of room temperature (23° Celsius) and with all participants off of drugs influencing the venous tone. Measurements were all performed in the morning hours following recommendation to drink 500 ml after the wake, in order to have comparable conditions of hydration [5].

Protocol of ECD measurement

Subject positioning and condition of measurement

Each experiment was performed with the subject placed on a tilting chair in both supine and upright positions while breathing normally, by starting the examination in supine position. After changing position, an adaptation period of at least 2 minutes was allowed before any further measurement.

The examiner carefully observed the inclination of the patient’s neck and provided appropriate neck support to avoid neck flexion, hyperextension or rotation to the left or right, which could potentially compress the neck veins and consequently affect measurements.

We used a thick layer of ultrasonic gel as well as recommended maneuvers in order to reduce excessive pressure on the patient’s neck that may change the shape and size of the IJV [5].

Evaluation of Doppler venous haemodynamics

Total inflow and outflow volume per unit of time, namely the flow Q, was measured in both supine and sitting position for each volunteer.

Inflow has been measured at common carotid artery (CCA) just below the bifurcation, and at the proximal segment of both internal carotid (ICA), and external carotid (ECA) artery.

The vertebral artery (VA) was evaluated at V2 level for reproducibility reasons [5]. In addition, outflow was measured in sequence at J2, J3, J1 level of the IJVs and at C4-C5 level of the VVs [5].

Two different approaches were used to calculate inflow and outflow, depending on the different shapes of the cervical arteries and veins. The investigated arteries have almost circular cross sectional area (CSA), so the CSA in this case was calculated using the diameter measured in longitudinal aspect of the B-mode imaging. Therefore, the Doppler sample volume was placed in the artery with the sample aperture corresponding to the lumen, in order to perform flow measurements by means of uniform insonation techniques [810].

On the contrary, since the IJV exhibits an elliptical shape, the CSA and major axis were assessed in the transversal aspect of the B-mode imaging by manually tracking the boundaries of the lumen.

The mean velocity of the blood in the veins has been calculated in accordance with the assumed profile technique, namely by sampling such velocity at the point where its value is highest. To do this, a smaller sample volume of 0.5 mm has been adopted for three reasons: 1) the small sample volume assures that the Doppler angle is constant over the whole sample volume, 2) the use of constant sample volume simplifies very much the on-line work of the Doppler operator and thus assuring a more accurate measurement 3) the use of small sample volume minimizes the vessel wall artifacts [812].

In the segment J1 of the IJVs, when we observed either an absent or a turbulent flow, we calculated the velocity profile by finely sampling in five different positions from wall to wall.

For both arteries and veins, the duration of the acquired Doppler spectrum was 4 seconds. For the arteries we considered three cardiac cycles while for the veins one respiratory cycle.

Off-line assessment of Doppler haemodynamics

We carefully acquired images and traces as above described, trying to improve as much as we could the reliability of the Doppler assessment and of the variables determined by the operator (Angle, PRF, etc.). In particular, the actual Doppler angle has been always carefully checked off-line and the contour of the jugular cross section has been determined by observing the movement of the vessel wall during the respiratory cycle.

Actual measurements were carried out by EM while during the acquisition EM and FS agreed on the Doppler technique regarding angle, position of the SV, etc. Since measurements of both inflow and outflow took a long time, calculation of the haemodynamic parameters was performed off-line by using the stored images, in order to shorten the examination time and to avoid possible physiological changes.

The flow Q was calculated as Q = TAV × CSA, where TAV is the time average velocity of the blood when considering one respiratory cycle for the veins and three cardiac cycles for the arteries. TAV was calculated as TAV = TA V p × 1 η , where TAVp corresponds to the average velocity measured on the peak of the trace and η is the velocity factor [13] calculated following Vergara [14] and using the Womersley number [15].

Refinement of Doppler haemodynamics assessment

Off-line calculation permitted also to improve the accuracy of the derived parameters. Post-processing allowed us to record the minimum and maximum CSA during respiratory cycle by manual tracing. After that, the venous flow Q was determined by calculating the mean value of the CSA.

A second parameter needing accurate post-processing verification is the angle of the Doppler beam for the vessels under measurement (Doppler angle). Such parameter and the uncertainty of the operator in placing it usually affect the TAV assessment. In our off-line processing we managed to estimate the uncertainty of TAV measurements as described in [16]:

δ TAV = TAV × cos θ cos θ + ϵ - 1

where θ is the incident angle of the Doppler beam, and ϵ is the uncertainty of the operator. The uncertainty of the flow is given by: δQ = δTAV × CSA.

Parameters of head and brain circulation

All the measurements for the above mentioned arteries and veins have been taken on both right and left sides. In particular, the carotids have been measured in the CCA segment, in the ICA segment and in the ECA segment. In order to minimize the experimental error we assume that the total head blood inflow (HBinF) is:

HBinF = CCAs + ICAs + ECAs 2 + VAs

whereas the cerebral blood flow (CBF) was roughly assumed to be the sum of ICAs and VAs contribution and then calculated as the sum of ICAs and VAs flows [17]. The cerebral venous outflow (CVO) was calculated as the sum of the flow measured at level J3 of the IJVs and the flow measured in the VVs. The total head blood outflow (HBoutF) was calculated as the flow of both left and right IJVs at J1 plus the VVs flows.

Model of neck veins

In order to analyze the results we propose a haemodynamic model (Figure 1) which includes the neck pathways of the cerebral venous return.

Figure 1

Model of the neck pathway concerning the cerebral venous return. a) Red tubes represent inflow arteries vessels while blue tubes represent outflow venous vessels. The dashed line are used to represents the collateral venous network. b) Representation of the flow in the vessel of right network of the neck. The direction of the arrows indicate the physiological direction of the flow.

As shown in Figure 1a the red tubes represent the inflow vessels (CCAs, ICAs, ECAs and VAs) while the blue ones represent the outflow vessels (VVs e IJVs). The CCA is divided in ICA and ECA. VAs and ICAs enter the brain compartment (Brain-C) and then the flow is normally drained by IJVs, VVs and collateral veins which are represented in the model by blue coloured tubes. The ECAs enter the facial and neck compartment (FN-C) and then are mainly drained by the collateral veins of the face and neck.

Vessel drawn with a continuous line are those evaluated by ECD in this study (CCAs, ICAs, ECAs, IJVs e VVs) while the ones with a dotted line are collateral veins which have been inserted in our model to account for the variations of the jugular flow.

In Figure 1b, flow directions are represented by a continuous arrow: QJ3, QJ2 and QJ1 are the measured flows in J3, J2 and J1 respectively, while Qvv is the measured flow in VV. Figure 1b also shows the collateral flows by means of a dotted arrow.

From top to bottom we now describe in detail the collateral flows of Figure 1b.

QC-D (Collateral-Distal) is the brain outflow which goes directly into the collateral network:

Q C - D = CBF - CVO

QFN, is the flow coming from the facial-neck compartment and going again into the collateral network:

Q FN = HBinF - CBF

Q23 is the collateral flow entering the IJV between J2 and J3:

Q 23 = Q J 2 - Q J 3

The above definition also applies to Q12:

Q 12 = Q J 1 - Q J 2

For both Q23 and Q12, we define a positive flow when it has a direction from a collateral towards the jugular. QC-P (Collateral Proximal) is the collateral outflow which goes directly into the caval system:

Q C - P = HBinF - HBoutF

Finally, all the flows defined above satisfy the continuity equation:

Q C - D + Q FN - Q 23 + Q 12 + Q C - P = 0

Calculated indexes

The above measured haemodynamic parameters also allow to extrapolate four indexesa:

  1. 1.

    Delta Cerebral Venous Outflow (DCVO), defined as:

    DCVO = Q J 1 s + Q VVs HBinF | Supine - Q J 1 s + Q VVs HBinF | Upright × 100

This index represents the normalized outflow difference between the supine and the upright position, as measured at the J1 level.

  1. 2.

    Distal Jugular and Vertebral Draining Index (DJVDI), defined as:

    DJVDI = CVO HBinF × 100

This index represents the percentage of the blood entering in the head that is drained directly from the IJVs at level J3 and from the VVs.

  1. 3.

    Collateral Flow Index (CFI), defined as:

    CFI = Q C - P HBinF × 100

This index represents the percentage of the blood entering in the head that is drained from collateral vessel instead to be drained from the IJVs or from the VVs.

  1. 4.

    Cerebral Collateral Draining Index (CCDI), defined as:

    CCDI = Q C - D CBF × 100

This index represents the percentage of the blood entering the brain that is drained from collateral vessels instead to be drained from the IJVs or from the VVs. The suffix ′s′ in VVs, J1s and J3s indicates that both left and right flow are considered.

Phase two of the study

We tested our model on a second population represented by ten patients (age ranging from 37 to 45 y.o., male:female ratio 5:5) affected by CCSVI. Such patients have been screened by the same ECD criteria among those affected by multiple sclerosis. All the selected patients showed a positivity of criterion 3, (i.e. presence of documented intra-luminal obstacles such as septa, membranes, webs, etc.) [5].

Statistical analysis

Data are expressed as mean ± sd. The haemodynamic parameters were analyzed either separately in the different jugular sides, or as a whole. Differences among the parameters assessed in both healthy volunteers and in CCSVI patients were tested by means of Wilcoxon-Mann–Whitney U-test; p value < 0.05 was considered significant.

Informed consent

The entire cohort of investigate subjects was informed about the methods and purpose of the experimental procedure and agreed to participate by signing an informed consent form. This study was in accordance with the Ethical Standards of the Committee on Human Experimentation of the University of Ferrara.


First phase of the study

Arterial inflow

The control subjects were successfully investigated. Calculated HBinF was 956 ± 105 ml/min, subdivided in 843 ± 200 ml/min in the CCAs, 462 ± 90 ml/min in the ICs, 255 ± 59 ml/min in the ECAs and 176 ± 72 ml/min in the VAs. Such values are similar to what was previously reported [17].

Venous outflow indexes

The calculated DCVO, DJDVI, CFI and CCDI indexes values are reported in Table 1 with their standard deviation and their meaning are discussed throughout the following text.

Table 1 Mean value (v) and standard deviation (sd) for Delta Cerebral Venous Outflow (DCDVI), Distal Jugular and Vertebral Draining Index (DJVDI), Collateral FLow Index (CFI) and Cerebral Collateral Draining Index (CCDI)

Venous outflow in supine posture

In Table 2 we report CSA, major axis, TAV, and Q respectively for right and left IJV, measured in supine; they increase from J3 to J1, and in J2 these values correspond to what was previously reported [18]. TAV increases significantly from J2 to J1, leading of course to increased Q values. The Q measured in left and right J2 plus VVs is about 11% less than the value reported by Doepp et al. [4]. However, this is coherent with a 14% higher CBF measured by the same authors in their normal subjects.

Table 2 Mean values (v) and standard deviation (sd) for cross sectional area (CSA), Major axis, TAV and flow (Q) in Jugular and Vertebral veins

In our sample, the rate of HBinF drained by the IJVs is 37% in J3, 55% in J2 and more than 90% in J1, respectively, and thus suggesting a re-entry of significant blood volume along the jugular vein through the collaterals.

It is worth noting that more than 90% of HBinF is drained by the IJVs in upright posture. Although there is evidence in the literature that VVs are the main draining route in this position, our finding refers to measurements in J1, a segment not previously investigated. Since this is a preliminary study that refers to a small sample size, it is important to investigate the current finding so as to determine the exact role of the gravitational gradient [1, 2] in the distribution changes of venous outflow from the brain.

In addition, our model permits to derive the volume of blood flowing into the collaterals of normal subjects, through the methodology reported above. As shown in Table 2, this is a consistent amount of blood never measured before: up to 350 ml/min for the collaterals entering between J2 and J3 and more than 500 ml/min for the collaterals entering between J1 and J2. However, the mean measured CFI was 1 ± 3%, clearly indicating that a very little fraction of blood flowing along the collaterals of normal subject bypasses the IJV and re-enters directly into the caval system.

The index DJDVI and CCDI were respectively 45 ± 17% and 33 ± 24%. The DJDVI reveals that for healthy controls in upright position, 45% of the mean HBinf is drained both by the IJVs at the J3 level and the VVs. Concerning the CCDI index, we found that about 33% of the CBF is drained through the collaterals. However, since the CFI is only 1%, this blood always flow into the jugulars.

Venous outflow in upright posture

In Table 2 we report CSA, major axis, TAV, and Q respectively for left and right IJV as measured in upright; TAV increases from J3 to J1, whereas CSA and major axis are apparently constant.

In our sample the rate of HBinF drained by the IJVs is 26% in J3, 33% in J2 and more than 90% in J1 and thus suggesting, also in upright, a re-entry of significant volume of blood along the jugular vein through the collaterals. As previously reported [3], we measured a significant reduction of the sum of the jugular and vertebral outflow in J2 when comparing the sitting with the supine position (mean 448 ml/min vs 600 ml/min).

Finally, the index DJDVI and CCDI were respectively 41 ± 10% and 39 ± 16%, while DCVO value was 5 ± 10%. In this case the DJDVI reveals that for healthy controls in upright position, 41% of the mean HBinf is drained both by the IJVs at the J3 level and the VVs. Concerning the CCDI index, we found that about 40% of the CBF is drained through the collaterals.

Second phase of the study

Arterial inflow

All the patients were successfully investigated. Calculated HBinF was 908 ± 90 ml/min subdivided in 758 ± 138 ml/min in the CCAs, 444 ± 123 ml/min in the ICs, 230 ± 83 ml/min in the ECAs and 192 ± 60 ml/min in the VAs.

Venous outflow in supine posture

In Table 2 we report CSA, major axis, TAV, and Q respectively for left and right IJV, measured in supine. Differently from what we measured in control subjects, Q and TAV increased from J3 to J2 but not from J2 to J1. This is confirmed by the rate of the HBinF drained in the different segments of the IJV, respectively 24% in J3 and 54% in J2, but dramatically reduced to 32% in J1. Since CFI is 61 ± 27%, our model permits to discover a significant volume of blood flowing in the collateral network rather than in the terminal segment of the IJV.

Finally, the index DJVDI and CCDI were respectively 33 ± 17% and 53 ± 23%. The high CCDI value shows that a significant fraction of the CBF is drained by the collaterals rather than the main routes (IJV e VV).

Venous outflow in upright posture

In Table 2 we report CSA, major axis, TAV, and Q respectively for left and right IJV, measured in sitting; TAV and Q increases from J3 to J1, whereas CSA and major axis are apparently constant. The Q measured in left and right J2 plus VVs is about 510 ml/min.

The rate of HBinF drained in the different segments of the IJV is 32% in J3, 41% in J2 and more than 80% in J1. Besides, we found a consistent amount of blood, more than 500 ml/min, for both the collaterals entering between J2 and J3 and the collaterals entering between J1 and J2. The CFI is 33 ± 31%, so considerably lower than the supine position and thus indicating a reduction of blood circulating into the collateral network when the drainage occurs in favour of gravity.

Finally, the index DCVO DJDVI and CCDI were -42 ± 80%, 41 ± 24%, and 40 ± 31%, respectively.

Comparison between healthy controls and CCSVI patients according to the model

The HC cohort was compared to the CCSVI one. It is worth noting that HBinF and CBF did not show significant differences among the groups (p > 0.14 and p > 0.95 respectively), hence, permitting a more focused comparison of the differences of cerebral venous return between the two groups. From this point of view, the main difference is the flow in J1 which, for the CCSVI patients, is about 70% less than the healthy controls (p < 0.001). Consequently, in the latter we found a significant higher CFI (p < 0.0002), clearly indicating the level of activation of the collateral network in the latter group.

Comparison in supine position

The above results are the consequence of the significant flow differences measured between the two groups in the supine position. While in J3 the flow Q showed simply a trend (p = 0.07), in J1 both Q and the CFI dramatically decreased (p < 0.000002). The latter result depends on the fact that the CFI index for healthy controls is separated by two standard deviations from the CFI of the MS patients (see Table 1).

Comparison in upright position

By turning the subjects in sitting posture, we did not find out significant differences in the control group by comparing the flow in the two postures. The major limitation is linked with the small sample and the big sd. To the contrary, by turning the CCSVI patients from supine to upright there is a drop in the the jugular flow in J1.


First phase of the study

In the first part of the study we tested the model on a HC cohort based on medical history and a controversial US CCSVI screening [5, 1927]. However, a recent meta-analysis clearly shows that the majority of HC are not affected by CCSVI [28]. Finally, also MRI data, more objective and less operator dependent with respect to US, are still controversial because there are confirmatory and not confirmatory studies [2527, 29].

Our measurements of the inflow are definitively comparable with previously published data [17]. Same result was found for the evaluation of the outflow, because the Q value assessed in J3, J2 and in the VVs are similar to the values previously reported [3, 4, 19].

The novelty of the present study is the application of a complete model which takes into account the haemodynamics of cerebral venous return normalized to the HBinF. Our model, for the first time, also includes J1 and haemodynamic analysis of collaterals.

Furthermore, we confirm that the flow in the IJV increases from the jaw to the chest [6], with consequent increased rate of the initial HBinF which is drained by the three considered segments. This is likely due to the re-entry of the collaterals into the main outflow route, as demonstrated by the calculated part flowing in the collateral network. However, we underline that in HC only 1% of the HBinF was not measured in the final amount of the HBoutF, thus indicating that a very small amount of blood volume in physiology re-enters through the collaterals into the caval system by skipping the IJV.

However, even if our model is complete and permits to better detail the modality of drainage from the brain, we would discuss some potential shortcomings linked to the proposed experimental setting. The first observation is that, due to the work of the valve leaflets, in J1, also in physiologic condition, the flow is turbulent. Flow turbulences may potentially affect the measurements of the TAV, so resulting in an overestimated assessment of Q in J1 [30]. This means that the assessment is less precise with respect to a straight venous segment having an ideal laminar flow and this issue will be subject of future work.

A second limitation in the ECD assessment of Q in J1 is linked to the technical feasibility to place a steering angle different from 0° in the lumen, especially when insonating the supra-clavicular fossa in subjects with more pronounced clavicular bone [31].

Moreover, we estimated an uncertainty of about 5 degrees in ϵ when the operator places the sample volume into the J1 lumen. This uncertainty has been estimated by asking the operator to recursively assess the correct Doppler angle so as to evaluate the standard deviation of the mean. The variability of such technical aspect may potentially lead to an overestimation of the TAV, finally affecting the Q up to 20% [16].

Finally, statistical comparisons were not adjusted for demographics and vascular risk factors, but this analysis is beyond the aim of the present study where we tested the feasibility of the proposed model in a limited sample size.

Second phase of the study

Once we developed the above described model, we tested its potential utility in the clinical setting by performing preliminary measurements in CCSVI condition. The main finding of the second phase is the significantly higher fraction of blood flowing in the collateral network of the CCSVI patient with respect to the HC. Our model permits to extrapolate that about 60% of the initial HBinF is transported directly to the caval system, significantly higher than 1% of CFI assessed in HC. This quantity dramatically increases because does not include only the flow drained in the soft tissue of the face and neck, but likely a high rate of blood transported by the IJV. Such a vision is clearly supported by two measured parameters. The former is the decreased Q passing from J2 to J1 (about 40%). The latter is the negative flow measured in C1-2 in the same population, which indicates the inversion of flow direction in the upper collaterals and it is likely due to the increased resistance exhibited by the terminal jugular vein. This anomalous behaviour could be the consequence of the intra-luminal obstacles detected in J1 at the time of ECD screening.

As an example, we applied the proposed model to compare HC subjects with CCSVI ones having same age and gender. Comparing Figure 2a with Figure 2b it is apparent that the flow at J3 and J2 are comparable, as well as the amount of blood flowing in the collaterals.

Figure 2

Comparison between the neck pathway of a) HC subject and b) CCSVI subject. The numbers in parenthesis refer to the blood flow in ml/min. The dashed arrows in red colour have an opposite direction.

What is dramatically changed is the value of Q in J1, where we assessed in the control subjects a further increase which leads to an overall amount of about 90% of the HBinF. On the contrary, a flow decrease of about 50% is apparent when passing from J2 to J1 in the CCSVI case. Finally, also the application of the model in venous pathology clearly shows the same limitations described above.

Comparison in supine position

We proposed four novel parameters in order to characterize the cerebral venous return, but only the CFI showed significant differences between the two cohorts in our study. CFI expresses the blood flowing into collateral network rather than the main outflow routes. It is really interesting that CCDI, which represents the rate of the inflow going into the collaterals at J3 level, is not significantly different in the two cohorts. This result may be linked to the limited flow assessed in the IJV at J1 level.

Comparison in upright position

The main finding is represented by the considerable drop of IJV flow measured in J1 when changing the position of MS patients, as also previously assessed by Doepp [19] and Monti [21].


We developed a new model that permits a detailed ECD quantification of the cerebral venous return, including an estimation of the amount of blood flowing from the collaterals to the caval system or to the IJV. The preliminary application of the model seems to indicate how a significant rate of the head inflow is drained by the collateral network rather than by the IJV in the CCSVI condition. This may help the interpretation of several findings assessed with different techniques, where it was not possible to assess the outflow contribution of the collateral network, as well as the rate of the inflow going in the main venous paths. For instance, the higher flow in the collateral network may explain the longer cerebral circulation time measured by means of contrast-enhanced US, as well as the slower discharge and increased resistance measured in MS [3234]. Our preliminary report needs to be further corroborated by reproducibility analysis, wider number of subjects and pathological conditions, and possibly, by a multicenter design. This may lead to a further advancement for the circulatory quantification of the CCSVI condition in the clinical setting via ultrasonography.


aThe named indexes are subject to copyright.



Brain compartment


cerebral blood flow


Common carotid artery


Collateral Flow Index


Cross sectional area


Cerebral Collateral Draining Index


Chronic Cerebrospinal venous insufficiency


Cerebral venous outflow


Delta Cerebral Venous Outflow


Distal, Jugular Vertebral Draining Index


External carotid artery


Echo coulor Doppler


Head blood in-flow


Head blood out-flow


Healthy control


Internal carotid


Internal jugular veins


Total of patients


Facial and neck compartment


Sample volume


Time average velocity


Vertebral artery


Vertebral veins.


  1. 1.

    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-327. 10.1113/jphysiol.2004.070409.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Schaller B: Physiology of cerebral venous blood flow: from experimental data in animals to normal function in humans. Brain Res Rev. 2004, 46: 243-260. 10.1016/j.brainresrev.2004.04.005.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Valdueza JM, von Münster T, Hoffman O, Schreiber S, Einhäupl KM: Postural dependency of the cerebral venous outflow. Lancet. 2000, 355: 200-201.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Doepp F, Schreiber SJ, von Münster T, Rademacher J, Klingebiel R, Valdueza JM: How does the blood leave the brain? A systematic ultrasound analysis of cerebral venous drainage patterns. Neuroradiology. 2004, 46: 565-570.

    Article  PubMed  Google Scholar 

  5. 5.

    Zamboni P, Morovic S, Menegatti E, Viselner G, Nicolaides AN, The Intersociety Faculty: Screening for chronic cerebrospinal venous insufficiency (CCSVI) using ultrasound. Recommendation for a protocol. Int Angiol. 2011, 30: 571-597.

    CAS  PubMed  Google Scholar 

  6. 6.

    Zamboni P, Menegatti E, Pomidori L, Morovic S, Taibi A, Malagoni AM, Cogo AL, Gambaccini M: Does thoracic pump influence the cerebral venous return?. J Appl Physiol. 2012, 112: 904-910. 10.1152/japplphysiol.00712.2011.

    Article  PubMed  Google Scholar 

  7. 7.

    Zamboni P, Consorti G, Galeotti R, Gianesini S, Menegatti E, Tacconi G, Carinci F: Venous collateral circulation of the extracranial cerebrospinal outflow routes. Curr Neurovasc Res. 2009, 6: 204-212. 10.2174/156720209788970054.

    Article  PubMed  Google Scholar 

  8. 8.

    Evans DH: On the measurement of the mean velocity of blood flow over the cardiac cycle using Doppler ultrasound. Ultrasound Med Biol. 1985, 11: 735-741. 10.1016/0301-5629(85)90107-3.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Maulik D: Doppler ultrasound in obstetrics and gynaecology. 1996, Heidelberg: Springer

    Google Scholar 

  10. 10.

    Hoskins P, Thrush A, Martin K, Whittingam T: Diagnostic Ultrasound: Physics and Equipment Diagnostic, Greenwich Medical Media. 2010, New York, NY: Cambridge University Press

    Google Scholar 

  11. 11.

    Ponzini R, Vergara C, Rizzo G, Veneziani A, Roghi A, Vanzulli A, Parodi O, Redaelli A: Womersley number-based estimates of blood flow rate in Doppler analysis: in vivo validation by means of phase-contrast MRI. IEEE Trans Biomed Eng. 2010, 57: 1807-1815.

    Article  PubMed  Google Scholar 

  12. 12.

    Li S, Hoskins PR, Anderson T, McDicken WN: Measurement of mean velocity during pulsatile flow using time-averaged maximum frequency of Doppler ultrasound waveforms. Ultrasound Med Biol. 1993, 19: 105-113. 10.1016/0301-5629(93)90002-6.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Zhong Z, Song H, Chui TY, Petrig BL, Burns SA: Noninvasive measurements and analysis of blood velocity profiles in human retinal vessels. Invest Ophthalmol Vis Sci. 2011, 52: 4151-4157. 10.1167/iovs.10-6940.

    Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Vergara C, Ponzini R, Veneziani A, Redaelli A, Neglia D, Parodi O: Womersley number-based estimation of flow rate with Doppler ultrasound: sensitivity analysis and first clinical application. Comput Methods Programs Biomed. 2010, 98 (2): 151-160. 10.1016/j.cmpb.2009.09.013. May

    Article  PubMed  Google Scholar 

  15. 15.

    Womersley JR: Flow in the larger arteries and its relation to the oscillating pressure. J Physiol. 1954, 124 (2): 31-32.

    CAS  PubMed  Google Scholar 

  16. 16.

    Yamamoto M, Carrillo J, Insunza A, Mari G, Ville Y: Error Introduced into Velocity Measurements by Inappropriate Doppler Angle Assignment. Ultrasound Obstet Gynecol. 2006, 28 (6): 853-854. 10.1002/uog.3847.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Schoning M, Walter J, Scheel P: Estimation of Cerebral Blood Flow Through Color Duplex Sonography of the Carotid and Vertebral Arteries in Healthy Adults. Stroke. 2004, 25: 17-22.

    Article  Google Scholar 

  18. 18.

    Mortensen JD, Talbot S, Burkart JA: Cross sectional Internal Diameters of Human Cervical and Femoral Blood Vessels: Relationship to Subject’s Sex, Age, Body Size. Anat Rec. 1990, 226 (1): 115-124. 10.1002/ar.1092260114.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Doepp F, Paul F, Valdueza JM, Schmierer K, Schreiber SJ: No cerebrocervical venous congestion in patients with multiple sclerosis. Ann Neurol. 2010, 68: 173-183.

    Article  PubMed  Google Scholar 

  20. 20.

    Bastianello S, Bergamaschi F, Viselner G, et al: An international multicenter observatory on prevalence of CCSVI in MS. BMC Neurol. 2011, in press

    Google Scholar 

  21. 21.

    Monti L, Menci E, Ulivelli M, Cerase A, Bartalini S, Piu P, Marotti N, Leonini S, Galluzzi P, Romano DG, Casasco AE, Venturi C: Quantitative ColourDopplerSonography evaluation of cerebral venous outflow: a comparative study between patients with multiple sclerosis and controls. PLoSOne. 2011, 6 (9): e25012-

    CAS  Article  Google Scholar 

  22. 22.

    Mayer CA, Waltraud P, Matthias WL, Nedelmann M, Bechmann I, Ziemann , Steinmetz : The perfect crime? CCSVI not leaving a trace in MS. J Neurol Neurosurg Psychiatry. 2011, 82: 436-440. 10.1136/jnnp.2010.231613.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Baracchini C, Perini P, Calabrese M, Causin F, Rinaldi F, Gallo P: No evidence of chronic cerebrospinal venous insufficiency at multiple sclerosis onset. Ann Neurol. 2011, 69: 90-99. 10.1002/ana.22228.

    Article  PubMed  Google Scholar 

  24. 24.

    Blinkenberg M, Akeson P, Sillesen H, Lövgaard S, Sellebjerg F, Paulson OB, Siebner HR, Sørensen PS: Chronic cerebrospinal venous insufficiency and venous stenoses in multiple sclerosis. Acta Neurol Scand. 2012, 126 (6): 421-427. 10.1111/j.1600-0404.2012.01671.x.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    McAuliffe W, Kermode AG: Mystery of Chronic Cerebrospinal Venous Insufficiency: Identical Venographic and Ultrasound Findings in Patients with MS and Controls. AJNR Am J Neuroradiol. 2013, Jan 31. [Epub ahead of print]

    Google Scholar 

  26. 26.

    Utriainen D, Trifan G, Sethi S, Elias S, Hewett J, Feng W, Haacke EM: Magnetic resonance imaging signatures of vascular pathology in multiple sclerosis. Neurol Res. 2012, 34 (8): 780-92. 10.1179/1743132812Y.0000000078. Oct

    Article  PubMed  Google Scholar 

  27. 27.

    Utriainen D, Feng W, Elias S, Latif Z, Hubbard D, Haacke EM: Using magnetic resonance imaging as a means to study chronic cerebral spinal venous insufficiency in multiple sclerosis patients. Tech Vasc Interv Radiol. 2012, 15 (2): 101-12. 10.1053/j.tvir.2012.02.003. Jun

    Article  PubMed  Google Scholar 

  28. 28.

    Laupacis A, Lillie E, Dueck A, Straus S, Perrier L, Burton JM, Aviv R, Thorpe K, Feasby T, Spears J: Association between chronic cerebrospinal venous insufficiency and multiple sclerosis: a meta-analysis. CMAJ. 2011, 183 (16): E1203-12.

    Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Zivadinov R, Schirda C, Dwyer MG, Haacke ME, Weinstock-Guttman B, Menegatti E, Heininen-Brown M, Magnano C, Malagoni AM, Wack DS, Hojnacki D, Kennedy C, Carl E, Bergsland N, Hussein S, Poloni G, Bartolomei I, Salvi F, Zamboni P: Chronic cerebrospinal venous insufficiency and iron deposition on susceptibility-weighted imaging in patients with multiple sclerosis: a pilot case-control study. Int Angiol. 2010, 29 (2): 158-75.

    CAS  PubMed  Google Scholar 

  30. 30.

    Pozniak MA, Zagzebski JA, Scanlan KA: Spectral Doppler and Color Doppler Artifacts. Radiographics. 1992, 12: 35-44.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Steinman AH, Yu AC, Johnston KW, Cobbold RS: Effects of beam steering in pulsed-wave ultrasound velocity estimation. Ultrasound Med Biol. 2005, 31 (8): 1073-1082. 10.1016/j.ultrasmedbio.2005.04.005.

    Article  PubMed  Google Scholar 

  32. 32.

    Mancini M, Morra VB, Di Donato O, et al: Multiple sclerosis: cerebral circulation time. Radiol. 2012, 262 (3): 947-955. 10.1148/radiol.11111239.

    Article  Google Scholar 

  33. 33.

    Beggs C, Shepherd S, Zamboni P: Cerebral venous outflow resistance and interpretation of cervical plethysmography data with respect to the diagnosis of chronic cerebrospinal venous insufficiency. Phlebology. 2012, Oct 11. [Epub ahead of print]

    Google Scholar 

  34. 34.

    Zamboni P, Menegatti E, Conforti P, Shepherd S, Tessari M, Beggs C: Assessment of cerebral venous return by a novel plethysmography method. J Vasc Surg. 2012, 56 (3): 677-685. 10.1016/j.jvs.2012.01.074.

    Article  PubMed  Google Scholar 

Pre-publication history

  1. The pre-publication history for this paper can be accessed here:

Download references


This study was partially supported by the Italian Ministry of Education, University and Research (MIUR Programme PRIN 2010–2011), Grant no. 2010XE5L2R.

Author information



Corresponding author

Correspondence to Paolo Zamboni.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

PZ conceived the study, contributes to develop the model and wrote the paper. FS developed the model and the equations, collected and analyzed the raw data, performed the statistical analysis and wrote the paper. EM performed the Echo Colour Doppler screening and collected the data. AT wrote the paper and revised it critically. AMM and SM contribute to collect and analyze the data. MG provided scientific supervision and founded the study. All authors participated in the design study. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Zamboni, P., Sisini, F., Menegatti, E. et al. An ultrasound model to calculate the brain blood outflow through collateral vessels: a pilot study. BMC Neurol 13, 81 (2013).

Download citation


  • Chronic cerebro-spinal venous insufficiency
  • Internal jugular vein
  • IJV
  • Echo colour doppler
  • Model
  • Ultrasound
  • Haemodynamics
  • Cerebral outflow