|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Departments of
1 Physiology
2 Internal Medicine, Academic Medical Center, Cardiovascular Research Institute Amsterdam, the Netherlands
3
Systems and Control Group, Faculty of Mechanical Engineering, Delft University of Technology, the Netherlands
4
The Copenhagen Muscle Research Center, Rigshospitalet, University of Copenhagen, Denmark
| Abstract |
|---|
|
|
|---|
(Received 21 June 2004;
accepted after revision 26 July 2004;
first published online 29 July 2004)
Corresponding author J. Gisolf: Department of Physiology, Room M01-07, Academic Medical Center, University of Amsterdam, PO Box 22700, 1100 DE Amsterdam, the Netherlands. Email: j.gisolf{at}amc.uva.nl
| Introduction |
|---|
|
|
|---|
We developed a mathematical model of the cerebral venous outflow tract to study the distribution of cerebral blood flow to the internal jugular veins and to an alternative pathway (vertebral venous plexus). Input to the model are beat-to-beat measurements of central venous pressure (CVP) and cerebral blood flow velocity (CBFV), with and without an increased CVP by a Valsalva manoeuvre, in the supine and standing position. For the internal jugular veins we implemented a nonlinear pressurevolume relationship based on measurements by Braakman et al. (1989). Knowing the jugular veins to be collapsed in the upright position and re-opened during positive pressure breathing, we hypothesized the cerebral outflow pathway to be dependent on posture and CVP. We expected the internal jugular veins to be the major drain for the brain in the supine position; in the standing position these veins are likely to be collapsed and to allow little blood flow. Furthermore we expected a 40 mmHg increase in CVP in the standing position to be sufficient to reopen the jugular veins and facilitate jugular blood flow for as long as CVP was raised. Ultrasound imaging of the internal jugular veins verified model outcome.
| Methods |
|---|
|
|
|---|
To assess the cerebral venous outflow distribution over the internal jugular veins and the alternate pathway (the cervical vertebral venous plexus) we developed a beat-to-beat model programmed in Simulink of MATLAB (Release 5.2, The MathWorks, Natick, MA, USA). A detailed description of the mathematical model and a parameter sensitivity analysis of the jugular vein model properties are given in the Appendix. Input data to the model are beat-to-beat measurements of CVP and CBFV variations. Lassen (1974) reported an average whole brain blood flow of 50 ml (100 g)1 min1 and considering a brain weight of 1500 g (Nelson, 1982), we assumed a supine cerebral blood flow of 750 ml min1 for all subjects. CBFV (blood velocity rather than blood flow) was set to this level at baseline in the supine position and used to track variations in cerebral blood flow. A schematic representation of the modelled flow distribution to the internal jugular veins and the vertebral venous plexus is shown in Fig. 1. The blood flow out of the brain is assumed to be equal to the blood flow (velocity) into the brain as measured in the middle cerebral artery, as the brain has little pooling capacity (Monro-Kellie doctrine) (Kellie, 1822). Pulsatility is omitted by using the beat-averaged mean CBFV. Per heartbeat, the model computes the flow distribution over both jugular veins, which are modelled as being identical, and the vertebral venous plexus, which is modelled as an invariable resistance. To implement the effect of a hydrostatic pressure gradient in the standing position, the model internal jugular vein was subdivided into a chain of 10 segments with a length of 1.5 cm each. The 9 most cranial segments each contain a variable resistor and capacitor, representing resistance to flow and the volume accommodating properties, respectively, while the most caudal segment contains a variable resistor but no capacitor (Fig. 1). In this caudal segment a jugular valve is implemented as a switch to a very high resistor (see Appendix), limiting flow reversal. The effects of gravity are implemented in the model as follows. The capacitance and resistance in each segment of the jugular veins are a function of the prevailing pressure in the segment (Braakman et al. 1989). The pressurevolume relation is highly nonlinear, with switch-like properties: at low (transmural) pressure, vessel volume is low; at high pressure the vessel volume approaches a maximal value (Fig. 2). In the standing position the pressure in each segment is corrected for hydrostatic pressure, and thus volume and resistance are computed from this level-corrected pressure. For example, in the standing position the most cranial part of the internal jugular vein is approximately 27 cm above heart level. This corresponds with a hydrostatic pressure correction of 21 mmHg. The height-corrected (low) transmural segment pressure will result in a low volume and a high resistance, in other words collapse of the vessel, and subsequent shunting of blood to the venous plexus. Approximate venous plexus resistance is set at 0.068 mmHg s ml1 (as derived from previous studies, see Appendix). CVP measurements are input to the model; in the standing position a substantial increase in CVP, such as occurs during a Valsalva manoeuvre, will result in increased filling of the internal jugular veins.
|
|
The physiological data we used to test our model consist of CBFV (in the middle cerebral artery) and CVP measurements in 10 healthy young subjects (4 women, age 27 ± 6 years, height 180 ± 10 cm, weight 76 ± 10 kg) who participated in the study of Pott et al. (2000) for which informed consent had been obtained from all participants, and which was approved by the ethics committee of Copenhagen (KF 01120/96) and was performed in accordance with the guidelines laid down in the Declaration of Helsinki. The procedures are fully described elsewhere (Pott et al. 2000); summarizing, after instrumentation the subjects rested in the supine position for 30 min. After a test run of a Valsalva manoeuvre, subjects rested for another 5 min and performed three Valsalva manoeuvres (40 mmHg expiratory pressure for 15 s), each followed by 3 min of recovery. Subjects were then asked to stand up in a relaxed position, and after 5 min they performed three Valsalva manoeuvres, each followed by 3 min of recovery. We analysed the first Valsalva in each body position.
The CBFV in the middle cerebral artery was measured using Transcranial Doppler (Multidop X, DWL, Sipplingen, Germany). Mean CBFV was computed as the integral of the maximal frequency shifts over one beat divided by the corresponding beat interval. Finger arterial pressure was measured with a Finapres model 5 (Biomedical Instrumentation, TNO-BMI). For CVP measurement, a catheter (1.7 mm i.d., 16 gauge) was placed in the superior caval vein through the basilic vein. CVP was recorded from a transducer (Bentley, Uden, the Netherlands) fastened to the subject in the midaxillary line at the level of the right atrium and connected to a monitor (8041, Simonsen & Weel, Copenhagen, Denmark). Measurements were computed after analog to digital conversion at a sampling rate of 100 Hz. Finger arterial pressure, CVP and CBFV were averaged beat-to-beat, and CBFV was further processed to approximate variations in total cerebral blood flow.
Experiments
To verify the model outcome of internal jugular vein cross-sectional area before and during the straining phase of the Valsalva manoeuvre in the supine and standing position, the following protocol was carried out in six healthy young subjects (3 women, age 31 ± 2 years, height 177 ± 3 cm, weight 70 ± 4 kg). Signed informed consent was obtained from all participants. The study was approved by the ethics committee of the Academic Medical Center (MEC 00/243) and performed in accordance with the guidelines laid down in the Declaration of Helsinki. Valsalva manoeuvres were conducted as described for Data Set; this was practiced 1 day prior to the experiments. Subjects rested in the supine position for 5 min before they performed a Valsalva manoeuvre. They were then asked to stand up and remain standing for 5 min, and perform a Valsalva manoeuvre in the upright position. Upper arm blood pressure was measured at 23 min into the supine and upright periods (Omron M5-I, oscillometric blood pressure monitor).
The cross-sectional area of the right internal jugular vein was imaged using ultrasound (Acuson Aspen 7.0). The ultrasound probe was placed on the neck approximately at the level of the laryngeal prominence, so that the probe was perpendicular to the vessel and the location was marked. The imaging depth was 4 cm, and the gain 65 dB. To avoid compressing the vein, care was taken to exert minimal pressure with the probe. The image of the venous lumen was frozen on the screen of the ultrasound unit and saved on optical disk. The cross-sectional lumen area of the internal jugular vein was determined off-line from the ultrasound image. Ultrasound measurements were conducted 1 min prior to (baseline) and during the straining phase of the Valsalva manoeuvre (10 s from the start) in the supine and in the standing position. Where pulsations in the jugular vein were detected, the image at mid-point of the pulsation was stored.
Statistical analysis
Data are presented as means ± S.E.M. Differences were tested by using paired t test; unequal variance (heteroscedastic) t test or sign test where appropriate. Pearson correlation coefficient was calculated for the correlation between the model estimates of the internal jugular vein cross-sectional areas and the ultrasound-determined areas. Group averages per event (baseline and Valsalva manoeuvre in the supine and standing position) were used. Significance was set at P = 0.05.
| Results |
|---|
|
|
|---|
The hemodynamic response to each Valsalva phase as defined from changes in blood pressure is described elsewhere (Pott et al. 2000); Table 1 gives the CVP and CBFV during baseline, phase IIa, IIb and IV. By definition (Sharpey-Schafer, 1965), in phase IIa (straining phase) of the Valsalva manoeuvre mean arterial pressure, pulse pressure and stroke volume decrease. This is followed by partial recovery of mean arterial pressure and heart rate toward the end of the strain, in phase IIb. Phase IV represents the arterial pressure overshoot after release of the strain. Cerebral venous outflow distribution over the internal jugular veins and the vertebral venous plexus was simulated using beat-to-beat variation in CVP and CBFV as input to the model (see Appendix for model description and settings). Figure 3 shows a representative example of CBFV and CVP measurements as well as simulated cerebral outflow distribution during supine and standing Valsalva manoeuvres; Fig. 4 shows simulation results averaged for all subjects.
|
|
|
Simulations of baseline flow and pressure in the supine position showed a cerebral venous outflow predominantly via the internal jugular veins (summed) and little flow via the vertebral venous plexus (Fig. 4). The pressure in the most cranial and most caudal jugular segments is given in Table 2; the pressure drop over the jugular veins in the supine position was only 0.2 mmHg. During phase IIb of the Valsalva manoeuvre, internal jugular vein volume and pressure increased (P < 0.01); the jugular veins remained the predominant cerebral outflow pathway.
|
Model simulation of baseline in the standing position showed a reduced internal jugular vein flow (Fig. 4; P < 0.01) and an increased venous plexus flow compared with supine baseline; the vertebral venous plexus was the main cerebral venous outflow pathway at baseline in the standing position. With a straining-induced increase in CVP (Valsalva), internal jugular vein volume and blood flow increased in phase IIb (P < 0.01 when compared with standing baseline); the jugular veins were the main cerebral venous outflow pathway during the straining phase of the Valsalva manoeuvre. Venous plexus blood flow decreased during straining in the standing position. Model estimates of the cross-sectional area of the internal jugular vein in the middle segment (Segment 5) during baseline and Valsalva manoeuvre phase IIa/b in the supine and standing position are given in Table 1.
Experiments
Upper arm arterial pressures, systolic and diastolic, were 120 ± 4 and 74 ± 1 mmHg, respectively, in the supine position and 120 ± 3 and 84 ± 3 mmHg, respectively, in the standing position. Table 1 gives the ultrasound-determined right internal jugular vein cross-sectional areas at the level of the laryngeal prominence. In the supine position the cross-sectional area was greater during a Valsalva manoeuvre compared with baseline (P < 0.05). In the standing position cross-sectional areas had diminished 6-fold, whereas a Valsalva manoeuvre while standing re-opened the vein (Fig. 5 shows a representative example). The ultrasound-determined right internal jugular veins areas, averaged per phase (baseline and 10 s from the start of a Valsalva manoeuvre approximating phase IIb) and per body position, correlated well with model predictions (R2 = 0.97).
|
| Discussion |
|---|
|
|
|---|
Collapsed internal jugular veins in upright humans imply an increased cerebral venous outflow resistance, and dependency on the vertebral venous plexus for cerebral venous return; this is illustrated in Fig. 6. During a Valsalva manoeuvre, elevation of CVP re-opens the jugular veins, which implies a reduction of cerebral outflow resistance. The important implication is that everyday events involving straining, such as actively standing up, or coughing and sneezing, affect the cerebral outflow pathway and total outflow resistance. Evaluation of the total cerebral resistance is outside the scope of this study; for this, changes in intracranial pressure should be taken into account. The brain is enclosed in a rigid compartment and there is limited space for volume expansion; an increased tissue pressure may induce passive collapse of the (intracranial) venous outflow vessels with the development of a venous outflow resistance (Kongstad & Grände, 1999).
|
Incidentally, advances in high-resolution two- and three-dimensional computed tomography have allowed assessment of cranial venous outflow pattern in a very different subject group; this technique is used to assess endocranial capacity of remnants of the early hominids. An enlarged occipital-marginal sinus system in evolving bipedal hominids has been suggested to be related to factors involved with more efficient blood flow to the vertebral venous plexus (Conroy et al. 1990).
Model considerations and study limitations
The vertebral venous plexus is rudimentary in our model; it consists of a single non-varying resistance. Considering the multitude of anastomoses in this complicated network of plexus in humans, a more realistic approach would be a resistor and a capacitor with significant pooling properties. Adding a capacitor to the vertebral venous plexus in the model would result in charging of the capacitor during increased CVP, featuring caudal-to-cranial blood flow. The direction of flow would be determined by the placement of the capacitor; for example, significant extracranial venous pooling properties would result in caudal-to-cranial blood flow during straining. Although simplification of the venous plexus is a model limitation, the present study focuses on the effect of hydrostatics on internal jugular vein collapse and subsequent shunting of blood to the alternative pathway, and a rudimentary model seems sufficient to test our hypothesis. Further investigation into the vertebral venous plexus requires refinement of the model as well as venous plexus blood flow measurements.
The internal jugular veins are modelled as being identical. Asymmetry of the jugular veins volume under baseline conditions would result in unequal distribution of flow: a greater venous outflow through the larger jugular vein. The collapse pressure of the veins would not be affected, however (see Appendix), and the overall distribution of jugular veins vs. alternative venous pathway would not be significantly altered. The nonlinear properties of the internal jugular veins are supported by a recent study addressing the absence of a siphon to support cerebral blood flow in standing humans (Dawson et al. 2004); collapse of jugular veins at tilt angles greater than 3035° are reported. Valdueza et al. (2000) reported a large internal jugular flow decrease already at 15° tilt (from the horizontal position). A 10° head-down position increases internal jugular vein cross-sectional area (Schreiber et al. 2002).
For the individual model simulations, subject height was not taken into account. The rationale to take the distance from the heart to the internal jugular veins as equal for all subjects was the nonlinear characteristic of the veins: for a height difference of some 30 cm between subjects, the heart-to-neck distance might differ only a few centimetres. At a CVP of around zero at heart level, in the standing position the hydrostatic pressure correction at neck level will be sufficient to ensure collapse of the veins, irrespective of subject height.
Conclusions
In conclusion, in humans cerebral outflow pathways include internal jugular veins and an alternative route, the vertebral venous plexus system. Results of mathematical modelling suggest that whereas internal jugular veins are the major drain for the brain in the supine position, in the standing position they are liable to collapse and cerebral venous blood is returned via the alternative pathway. During increased CVP in the standing position the internal jugular veins are re-opened and these veins are again the primary pathway for cerebral venous return. Ultrasound images of the internal jugular vein cross-sectional area verify model outcome.
| Appendix |
|---|
|
|
|---|
The venous plexus resistance (Rven plex) is indirectly derived from measurements by Cirovic et al. (2003). Using their measured jugular vein flow (Qjug) and resistance (Rjug) in the supine and sitting position in the following equations:
|
| (A1) |
|
| (A2) |
|
| (A3) |
Model equations
The structure of the cerebral venous outflow model is shown in Fig. 1. Inputs to the model are central venous pressure (CVP) and variations in cerebral blood flow (CBFV), which is computed by scaling beat-to-beat cerebral blood flow velocity (CBFV) data to a physiological range as described in the Methods section. Extra-vascular pressure in the neck is assumed to be atmospheric pressure. CBFV is distributed over two identical jugular veins and the cervical vertebral venous plexus.
Flow through the venous plexus is computed from the resistance in the venous plexus, the pressure at the beginning of the cerebral outflow tract (P) and the central venous pressure:
|
| (A4) |
The jugular vein is arbitrarily divided into 10 segments of equal length. Nine segments consist of a variable resistor and capacitor (Fig. 1), while the most caudal segment consists of a resistor only. For each segment X of the jugular vein, where X = 1 is the most cranial and X = 10 the most caudal segment, the hydrostatic pressure correction can be calculated as
ghX, where hX is the height difference between segment X and the heart,
is the density of blood and g is the gravitational acceleration. The height-corrected pressure (PHC) is only used for calculation of the resistance (RX), capacity (CX) and volume (VX) of each segment (eqns A57) using equations presented by Braakman et al. (1989) as illustrated in Fig. 2:
|
| (A5) |
|
| (A6) |
|
| (A7) |
is the viscosity of blood, P0 and V0 are parameters related to vascular tone and volume range, respectively, and A is the slope of the PV relation at V0.
The blood flow out of segment X (Qout,X) is computed from the inflow (Qin,X) and the volume change (dVX/dt):
|
| (A8) |
The pressure (non-height corrected) before the resistance of each segment, PX1 (Fig. 1), is calculated from Qin, RX and the pressure (PX) after the resistance:
|
| (A9) |
Parameter sensitivity analysis
To determine the sensitivity of the model to the input variables, we computed the pressurevolume and the pressureresistance relation using the input variables as were used for the model simulations, and by varying the parameters determining the pressurevolume relation by 10%. For P0, the zero-flow pressure intercept, which was set at 0, we computed 5 and +5 mmHg as well as 0. The results are shown in Fig. 2, where the resistance and volume are those in a jugular vein segment, as a function of the height-corrected pressure in the segment. The figure illustrates that P0 determines the volume of the vein segment at a certain segment pressure. This pressure is uncorrected for pressure in the neck; neck-suction for example will decrease the surrounding pressure in the neck, and the pressure in the segment can be corrected to obtain transmural pressure. Figure 2 also shows the influence of hydrostatic pressure differences. In transition from supine to standing, the reduction in segment pressure due to hydrostatics greatly influences the volume and resistance of the segment at pressures below 0; at high (central venous) pressure, such as occurs during straining, the height-corrected pressure in the vessel will be in the positive range regardless of posture.
|
| References |
|---|
|
|
|---|
Batson OV (1944). Anatomical problems concerned in the study of cerebral blood flow. Fed Proc 3, 139144.
Braakman R, Sipkema P & Westerhof N (1989). A dynamic nonlinear lumped parameter model for skeletal muscle circulation. Ann Biomed Eng 17, 593616.[CrossRef][Medline]
Chaynes P, Verdie JC, Moscovici J, Zadeh J, Vaysse P & Becue J (1998). Microsurgical anatomy of the internal vertebral venous plexuses. Surg Radiol Anat 20, 4751.[Medline]
Chou CH, Chao AC & Hu HH (2002). Ultrasonographic evaluation of vertebral venous valves. AJNR Am J Neuroradiol 23, 14181420.
Cirovic S, Walsh C, Fraser WD & Gulino A (2003). The effect of posture and positive pressure breathing on the hemodynamics of the internal jugular vein. Aviat Space Environ Med 74, 125131.[Medline]
Conroy GC, Vannier MW & Tobias PV (1990). Endocranial features of Australopithecus africanus revealed by 2- and 3-D computed tomography. Science 247, 838841.
Dawson EA, Secher NH, Dalsgaard MK, Ogoh S, Yoshiga CC, Gonzalez-Alonso J, Steensberg A & Raven PB (2004). Standing up to the challenge of standing: a siphon does not support cerebral blood flow in humans. Am J Physiol Regul Integr Comp Physiol (ahead of print).
Epstein HM, Linde HW, Crampton AR, Ciric IS & Eckenhoff JE (1970). The vertebral venous plexus as a major cerebral venous outflow tract. Anesthesiology 32, 332337.[Medline]
Holt JP (1941). The collaps factor in the measurement of venous pressure. Am J Physiol 134, 292299.
Kellie G (1822). On death from cold, and on congestion of the brain. Transactions Med-Chirurg Soc Edinburgh 1, 84169.
Kongstad L & Grände PO (1999). The role of arterial and venous pressure for Volume regulation of an organ enclosed in a rigid compartment with application to the injured brain. Acta Anaesthesiol Scand 43, 501508.[CrossRef][Medline]
Lassen NA (1974). Control of cerebral circulation in health and disease. Circ Res 34, 749760.
Nelson GE (1982). Fundamental Concepts of Biology. New York: Wiley, p. 262.
Pott F, van Lieshout JJ, Ide K, Madsen P & Secher NH (2000). Middle cerebral artery blood velocity during a valsalva maneuver in the standing position. J Appl Physiol 88, 15451550.
Schreiber SJ, Lambert UK, Doepp F & Valdueza JM (2002). Effects of prolonged head-down tilt on internal jugular vein cross-sectional area. Br J Anaesth 89, 769771.
Sharpey-Schafer EP (1965). Effect of respiratory acts on the circulation. In Handbook of Physiology. Circulation. Sect. 2, Vol. III. Chapt. 52, pp. 18751886. Washington, DC.
Toung TJK, Aizawa H & Traystman RJ (2000). Effects of positive end-expiratory pressure ventilation on cerebral venous pressure with head elevation in dogs. J Appl Physiol 88, 655661.
Valdueza JM, von Munster T, Hoffman O, Schreiber S & Einhaupl KM (2000). Postural dependency of the cerebral venous outflow. Lancet 355, 200201.[Medline]
Zippel KC, Lillywhite HB & Mladinich CR (2001). New vascular system in reptiles: anatomy and postural hemodynamics of the vertebral venous plexus in snakes. J Morph 250, 173184.[CrossRef][Medline]
Zouaoui A & Hidden G (1989). The cervical vertebral venous plexus, a drainage route for the brain. Surg Radiol Anat 11, 7980.[CrossRef][Medline]
| Acknowledgements |
|---|
This article has been cited by other articles:
![]() |
C. D. Etz, M. Luehr, F. A. Kari, C. A. Bodian, D. Smego, K. A. Plestis, and R. B. Griepp Paraplegia after extensive thoracic and thoracoabdominal aortic aneurysm repair: does critical spinal cord ischemia occur postoperatively? J. Thorac. Cardiovasc. Surg., February 1, 2008; 135(2): 324 - 330. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. H. Secher, T. Seifert, and J. J. Van Lieshout Cerebral blood flow and metabolism during exercise: implications for fatigue J Appl Physiol, January 1, 2008; 104(1): 306 - 314. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Rasmussen, P. Plomgaard, R. Krogh-Madsen, Y.-S. Kim, J. J. van Lieshout, N. H. Secher, and B. Quistorff MCA Vmean and the arterial lactate-to-pyruvate ratio correlate during rhythmic handgrip J Appl Physiol, November 1, 2006; 101(5): 1406 - 1411. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Gisolf, A. Gisolf, J. J. van Lieshout, and J. M. Karemaker The siphon controversy: an integration of concepts and the brain as baffle Am J Physiol Regulatory Integrative Comp Physiol, August 1, 2005; 289(2): R627 - R629. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |