Print version ISSN 1405-9940
Arch. Cardiol. Méx. vol.82 n.2 México Apr./Jun. 2012
Investigación clínica experimental
Relation of the baroreflex mechanism with the photoplethysmographic volume in healthy humans during orthostatism
Relación del mecanismo barorreflejo con el volumen fotopletismográfico en humanos sanos durante ortostatismo
Paola MartínezGarcía,1 Claudia Lerma,2 Óscar Infante.2
1 Graduate Studies in Medical Physics, Institute of Physics, UNAM. Mexico City, Mexico.
2 Department of Electromechanical Instrumentation, National Institute of Cardiology Ignacio Chávez. Mexico City, Mexico.
Oscar Infante Vázquez.
Departament of Electromechanical Instrumentation,
Instituto Nacional de Cardiología Ignacio Chávez.
Juan Badiano # 1, Sección XVI, Tlalpan. Z.P. 14080. Mexico City, Mexico.
Telephone: +52 (55) 5573 2911, ext 1386. Fax. +52 (55) 5573 0926.
Email address: email@example.com
Received on March 14, 2011.
Accepted on February 2, 2012.
Objective: Our aim was to evaluate the effect of the baroreflex mechanism upon peripheral blood volume during sympathetic stimulation by orthostatism.
Methods: Nineteen clinically healthy volunteers were included (12 men), 28.4 ± 6.2 years old. Blood pressure was monitored with a Finometer and blood volume with a photoplethysmograph during supine position and orthostatism (15 minutes each), in order to obtain systolic blood pressure (SBP), diastolic blood pressure (DBP), systolic volume (SysV), diastolic volume (DiaV), and inter beat intervals (IBI) measurements. Baroreflex sensitivity index (IBI/SBP) and baroreflex effect on blood volume (IBI/SysV) were estimated by the sequence method. The pertinence of using only systolic values was tested by linear regression analysis of systolic versus diastolic measurements.
Results: More than 70% of DBP and DiaV variations can be explained by SBP and SysV, respectively (p<0.001), with coherence >0.5 in frequencies between 0.04 and 0.15 Hz. IBI/SBP and IBI/SysV were linearly correlated (R>0.4) and both decreased during orthostatism (p<0.05).
Conclusion: The sequence method showed a strong baroreflex effect upon peripheral blood volume that became more apparent during sympathetic stimulation with orthostatism. This approach could be clinically useful for the evaluation of blood volume regulation for many diseases such as diabetes mellitus and heart failure, and during therapeutic interventions such as hemodialysis.
Keywords: Baroreflex mechanism; Blood volume regulation; Orthostatic challenge; Mexico.
Objetivo: Evaluar el efecto del mecanismo barorreflejo sobre el volumen sanguíneo periférico durante estimulación inducida por ortostatismo.
Métodos: Se incluyeron 19 voluntarios sanos (12 hombres), con edad de 28.4 ± 6.2 años. La presión arterial se midió con un Finometer y el volumen sanguíneo con un fotopletismógrafo, ambos durante posiciones supina y ortostatismo activo (15 minutos cada una), para obtener los valores de presión arterial sistólica (PAS), presión arterial diastólica (PAD), volumen sistólico (VS), volumen diastólico (VD) e intervalo inter pulso (IIP). Se estimó la sensibilidad barorrefleja (IIP/PAS) y el efecto barorreflejo sobre el volumen sanguíneo (IIP/VS) mediante el método de secuencias. La pertinencia de usar sólo variables sistólicas, se evaluó mediante análisis de regresión lineal de las mediciones sistólicas versus las diastólicas.
Resultados: Más de 70% de las variaciones de presión arterial diastólica y volumen diastólico pueden ser explicadas mediante presión arterial sistólica y volumen sistólico, respectivamente (p<0.001), con coherencia >0.5 en frecuencias entre 0.04 y 0.15 Hz. IIP/PAS y IIP/VS tuvieron correlación positiva (R>0.4) y ambos disminuyeron durante ortostatismo (p<0.05).
Conclusiones: El método de secuencias demostró un importante efecto barorreflejo sobre el volumen sanguíneo periférico que se hizo más notable durante estimulación simpática con ortostatismo. Este enfoque podría ser clínicamente útil en la evaluación de la regulación del volumen sanguíneo en distintas enfermedades como diabetes mellitus o falla cardiaca, y durante intervenciones terapéuticas como la hemodiálisis.
Palabras clave: Mecanismo barorreflejo; Regulación de volumen sanguíneo; Reto ortostático; México.
The baroreflex mechanism allows an adequate shortterm control of blood pressure (BP), in response to dailylife stimuli such as orthostatism.13 The neural and cardiovascular components of the baroreflex mechanism are complex, both anatomically and physiologically.4,5 However, the basic baroreflex function is often referred as a negative feedback mechanism, as it is illustrated in Figure 1. The reference variable (i.e. the BP) is "monitored" by the central nervous system through the baroreceptors, and in response to changes of BP, the central nervous system modulates the autonomic activity towards the heart and blood vessels, with the aim of increasing or decreasing BP by changes in heart rate (or inter beat interval, IBI) and peripheral arteries resistance.2,6 The baroreflex sensitivity is then defined as the magnitude of the change in IBI with respect to a change in BP,7 which normally takes effect between 1 to 3 seconds of the BP change (baroreflex delay).811
Impairment of the baroreflex mechanism, evidenced by decreased baroreflex sensitivity, occurs in many diseases such as essential arterial hypertension,1,12 diabetes mellitus,13 myocardial infarction14 and chronic renal failure.15 Two examples of clinical implications are: 1) the association of low baroreflex sensitivity with a high risk of cardiac mortality in patients with myocardial infarction,14 and 2) the observation of reduced baroreflex sensitivity during hypotension in patients with chronic renal failure prone to intrahemodialysis hypotension.16
Several noninvasive methods for baroreflex sensitivity evaluation have been developed,7,17,18 that have contributed to the understanding of the pathophysiology of the baroreflex mechanism. However, most studies are focused on the heart rate as the controlled variable, and little is known about the effect that the baroreflex has upon the peripheral blood volume. Considering that sympathetic regulation of peripheral resistance is one of the regulated variables necessary to control BP,6 we hypothesize that, given that the blood volume regulation has an important role on the cardiovascular control through the baroreflex mechanism19,20 this contribution must be seen in photoplethysmographic recordings. According to the representation proposed in Figure 1, correlated changes in systolic blood volume (SysV) and IBI (highlighted by a question mark) could be associated to sympathetic control on both variables. In that case, and considering no other confounding variables, the spectral coherence between both signals would be significant within the frequency band of 0.04 to 0.15 Hz. Furthermore, the effect of the baroreflex mechanism has been studied by systolic blood pressure (SBP) in most methods, but the equivalence of using systolic measurements or diastolic values has not been tested.
The aim of this work was to evaluate the effect of the baroreflex mechanism upon the peripheral blood volume as it is seen by photoplehysmography. Linear regression analysis was applied to the systolic and diastolic measurements of BP and photoplethysmographic volume (PV), to test the pertinence of using only the systolic values in the subsequent analysis. We applied a noninvasive method for assessment of baroreflex sensitivity to measurements of SBP, SysV and IBI, both in supine position and during active orthostatism.
The study included 19 volunteers (12 men and 7 women) with age 28.4 ± 6.2 years old and body mass index 26.5 ± 4.6 kg/m2. All participants were evaluated as healthy by a physician based on the medical history and physical examination only. The participants were not smokers, and were asked to avoid intake of medications or stimulants such as caffeine during 12 hours before the study. The study protocol was completely noninvasive, it complies with the ethical standards laid down in the 1964 "Declaration of Helsinki", and informed consent was obtained from each participant.
BP was recorded with a Finometer (Finapress Medical Systems),21 and PV was obtained by a previously validated photoplethysmograph.22 The PV signal was incorporated into a channel of the Finometer, in order to capture and digitize simultaneously both signals (BP and PV) with the same equipment. The photoplethysmographic signal is associated with blood volume changes in the microvascular bed of tissue.23 Since PV is not a direct quantitative measurement, it is represented in arbitrary units (a.u.). The sampling frequency was predefined by the Finometer equipment as 200 samples per second, which complies with the international recommendations for analysis of variability of heart rate and BP.24
BP measurements with the Finometer21 were calibrated for each patient by compressive sphygmomanometry. The pressure sensor (thimble) was secured on the middle finger and the PV sensor on the index finger, both in the left hand, with the arm resting on a sling to reduce hand movements. The first recording started after five minutes on supine position, and lasted 15 minutes. Then the patient changed to active orthostatism. The second recording started after five minutes of position change and lasted 15 minutes. All recordings were obtained between 11:00 am and 1:00 pm, and participants were instructed to breath normally at all recording times (spontaneous breathing).
BP and PV signals were analyzed with techniques implemented in Matlab (Mathworks®, Inc), as briefly described below.
The first step was the identification of each heartbeat and its corresponding values of SBP, diastolic blood pressure (DBP), SysV, diastolic blood volume (DiaV) and IBI. The BP signal was filtered by a second derivative algorithm,25 it was rectified to positive values and a threshold was applied in order to find the peaks that correspond to SBP (Figure 2A). Correct identification of the SBP for all heartbeats was verified manually. Then, the minimum value between consecutive SBP values was measured to obtain the DBP. The PV signal was filtered with a moving average filter with bandwidth of 0.02 to 10Hz, in order to improve the signal quality for better identification of SysV and DiaV in each heartbeat. Filtering the PV signal did not affect the frequency band of interest for SysV and DiaV (known as lowfrequency band, LF, from 0.04 to 0.15Hz).24,26 Then, the occurrence time for SBP and DBP of each heartbeat was used to locate the corresponding SysV and DiaV values from the PV signal. IBI was obtained by the time difference of SBP between consecutive heartbeats.
The systolic variables (SBP and SysV) were compared against their corresponding diastolic variables (DBP and DiaV) by bivariate linear regression analysis. For each recording, the prediction equation of simple linear model and the determination coefficient (R2) was obtained using the least squares method27 and comparing the diastolic signals (xaxis or independent variable) versus the systolic signals (yaxis or predicted variable) (Figure 2B). The statistical significance (p) of each R was evaluated for each regression line according to the size of the recording.
The value of R was considered as statistically significant if p<0.01, in this linear regression analysis.
In order to test the frequency range of the linear correlation between the systolic and diastolic variables, a spectral coherence analysis was performed. Since the sampling time all variables is not constant, they were resampled at 3 Hz to preserve their harmonics with frequencies up to 0.4 Hz.28 In order to maintain the same resolution in frequency (Sf) for all recordings, the systolic and diastolic resampled signals of BP and PV in supine position and orthostatism of each subject were cutoff to 2550 samples. Spectral leakage was decreased by fitting a straight line that was subtracted from the resampled signal. The power spectral density (PSD) of each patient was obtained by the Welch periodogram28 with a 50% overlap and a Hanning window of 300 points (i.e., window size of 100 s) (Figure 2C). The spectral coherence between diastolic and systolic signals of the BP and PV was calculated using Equation 1, where f is the frequency in Hertz, PSDxy denotes the crossed PSD between the diastolic signal (x) and the systolic signal (y), PSDx refers to the PSD of the diastolic signal (x), and PSDy is the PSD of the systolic signal (y).
Equation 1: Cxy(f) = |PSDxy(f)|2 / (PSDx(f) · PSDy(f)) A value of spectral coherence Cxy>0.5 was considered significant, and this was considered of interest within the LF band because blood vessels have no parasympathetic innervations (Figure 2C).4,6,26,29
The baroreflex effect estimation was performed as follows. Considering that arteries receive only sympathetic innervations,4,29 and the spectral content of interest falls within the LF band,24 a low pass filter with moving average cutoff frequency of 0.15 Hz was applied to SBP, SysV and IBI. The number of points needed to filter each signal was calculated with the sampling rate (approximated by the average IBI) divided by the cutoff frequency (0.15 Hz). Using the sequences method, each sequence comprised data points where three or more consecutive beats increase or decrease in SBP and IBI concurrently (Figure 3).7,30 The sequences were plotted (IBI vs SBP), and the least squares method was applied to obtain the slope of the linear fit for each sequence (Figure 3). The baroreflex sensitivity was then estimated as the average of all slopes, either in a positive (+IBI/+SBP) or negative (IBI/SBP) direction. The same procedure was applied to SysV and IBI signals, obtaining the baroreflex effect on SysV as the averaged values of all positive sequences (+IBI/+SysV) and negative sequences (IBI/SysV).
Results are given as mean ± standard deviation. Statistical analysis was performed using SPSS ver. 15 (StatSoft®, Inc.). Normal distributions of all variables were verified by KolmogorovSmirnov tests. The variables were compared by pairedt tests (supine position versus orthostatism). The linear correlation between the baroreflex sensitivity and the baroreflex effect on SysV was estimated by Pearson's correlation coefficient (R), considered statistically significant if p<0.05 and R>0.335 in each case (n=35).
The results of the linear regression analysis between systolic and diastolic measurements are shown in Table 1. In all cases (n=19) the correlation coefficient (R) was statistically significant (p<0.001) according to the number of data for each signal (between 900 and 1200 data points), with R2>0.7. Therefore, a significant amount of DBP variation can be explained by SBP and also a significant amount of variation of DiaV can be explained by SysV.
The spectral coherence analysis of diastolic versus systolic variables is summarized in Figure 4. From all cases, 15 subjects had a spectral coherence greater than 0.5 in the frequency band of interest (LF, 0.04 to 0.15 Hz), which confirms a strong linear association between the systolic and diastolic variables for each frequency within LF.
The effect of the orthostatic challenge on all evaluated variables is summarized in Table 2, including the baroreflex sensitivity index (i.e. IBI/SBP) and the estimated baroreflex effect upon PV (IBI/SysV). Compared to the supine position, orthostatism induced higher heart rate (i.e. shorter IBI), higher BP (both SBP and DBP), lower PV (both SysV and DiaV), lower baroreflex sensitivity and lower baroreflex effect upon PV.
Figure 5 shows the correlation coefficient (R) between the baroreflex sensitivity (IBI/SBP) and the baroreflex effect on SysV (IBI/SysV), for all subjects in supine position and orthostatism. In both type of sequences (positive and negative), a significant correlation between these two estimators was confirmed (p<0.05).
Orthostatism is a common activity in daily life that represents an important challenge for the cardiovascular control system. Normal response to orthostatism includes decreased baroreflex sensitivity,30 which is attributed to vagal withdrawal and sympathetic activation towards the heart,4,6 and sympathetic activation towards blood vessels.31 Our results are concordant with those observations, and showed also that peripheral blood volume experiences important change due to orthostatism, manifested by a mean decreased of both SysV and DiaV (which has already been reported)31 and by a decreased baroreflex effect on blood volume (IBI/SysV index). Our results support the proposal that the correlation of changes in Sys and IBI are mediated by the sympathetic nervous system,8,32 since both variables were significantly coherent at all frequencies within the frequency band of 0.04 to 0.15 Hz. Nevertheless, non all changes in SysV and DiaV observed by photopletysmography should be ascribed to sympathetically modulated vasoconstriction or vasodilation. Many other factors (such as respiration, nitric oxide, endothelin, hormones, and the reninangiontensinaldosteron sytem) could influence the variability of SysV and DiaV. Such factors could have local, regional or even central origin and some of them could be independent of the baroreflex mechanism. PV is a reflex of stroke volume, but its waveform and variability is conditioned by changes in the vessels properties (capacitance, inductance and resistance), which are dynamically affected by all mentioned factors including the sympathetic modulation. In this pilot study, none of these factors were controlled, which could explain the high dispersion within the indexes IBI/SBP and IBI/SysV (Figure 5). In spite of this, the current work demonstrated a significant linear correlation between these indexes.
Abnormal responses of the baroreflex reaction to orthostatism have been observed in patients prone to hypotension and syncope2 and in patients with hypertension either essential1,12 or secondary to chronic renal failure.15,16 In both extremes of abnormal BP, the main invoked culprit is the autonomic nervous system: syncope is preceded by increase of vagal tone and decrease of sympathetic activity (towards the heart, which produces bradycardia) and decreased sympathetic activity towards blood vessels2,8 (which decreases arterial resistance producing vasodilatation), while hypertension is associated to increased sympathetic activity both to the heart and blood vessels.12 Since variations of sympathetic activity to the arteries can cause significant changes in their blood volume, it is surprising that few studies have considered the variations in blood volume themselves as a key factor to explain the abnormal response to orthostatism31 and other stimuli such as hemorrhage,33 vasoconstriction induced by smoking,34 respiratory variations or apnea,35 local cold exposure,36 sympathetic blockade by epidural anaesthesia37 or even carotid baroreceptors denervation and thoracic sympathectomy.38 Although here we presented only the normal response of blood volume to one particular stimulus (orthostatism), we consider that this novel approach of estimation of the baroreflex effect on blood volume could be very relevant in further studies with different pathologies. The advantages of this approach include its noninvasiveness, the evaluation of the cardiovascular dynamics on a beattobeat basis, and the requirement of relatively shorttime measurements.
The main contribution of this work is the evidence of a strong baroreflex effect upon the peripheral blood volume, which was more apparent during stimulation of the cardiovascular system with an orthostatic challenge. This approach could be clinically useful for the evaluation of blood volume regulation for many diseases such as diabetes mellitus and heart failure, and during therapeutic interventions such as hemodialysis.
The authors thank Brayans Becerra for the technical assistance during the recording of signals. During the development of this work, P. MartínezGarcía was supported by CONACyT with a scholarship for graduate studies (Masters Degree on Medical Physics).
1. Lucini D, Pagani M, Mela GS, et al. Sympathetic restrain of baroreflex control of heart period in normotensive and hypertensive subjects. Clin Sci (Lond) 1994;86:547556. [ Links ]
2. Mano T, Iwase S. Sympathetic nerve activity in hypotension and orthostatic intolerance. Acta Physiol Scand 2003;177:359365. [ Links ]
3. Timmers HJ, Wieling W, Karemaker JM, et al. Denervation of carotid baro and chemoreceptors in humans. J Physiol 2003;553:311. [ Links ]
4. Estañol B, PorrasBetancourt M, SanchezTorres G, et al. Neural control of the peripheral circulation and blood pressure. Arch Cardiol Mex 2009;79:109116. [ Links ]
5. Persson PB. Modulation of cardiovascular control mechanisms and their interaction. Physiol Rev 1996;76:193244. [ Links ]
6. Heesch CM. Reflexes that control cardiovascular function. Am J Physiol 1999;277:S234S243. [ Links ]
7. Parati G, Di Rienzo M, Mancia G. How to measure baroreflex sensitivity: from the cardiovascular laboratory to daily life. J Hypertens 2000;18:719. [ Links ]
8. Gulli G, Claydon VE, Cooper VL, et al. RR intervalblood pressure interaction in subjects with different tolerance to orthostatic stress. Exp Physiol 2005;90:367375. [ Links ]
9. Gulli G, Cooper VL, Claydon VE, et al. Prolonged latency in the baroreflex mediated vascular resistence response in subjects with postural related syncope. Clin Auton Res 2005;15:207212. [ Links ]
10. Westerhof BE, Gisolf J, Stok WJ, et al. Timedomain crosscorrelation baroreflex sensitivity: performance on the EUROBAVAR data set. J Hypertens 2004;22:12591263. [ Links ]
11. Westerhof BE, Gisolf J, Karemaker JM, et al. Time course analysis of baroreflex sensitivity during postural stress. Am J Physiol 2006;291:H2864H2874. [ Links ]
12. Scholze A, Burkert A, Mardanzai K, et al. Increased arterial vascular tone during the night in patients with essential hypertension. J Hum Hypertens 2007;21:6067. [ Links ]
13. Buchs A, Slovik Y, Rapoport M, et al. Rightleft correlation of the sympathetically induced fluctuations of photoplethysmographic signal in diabetic and nondiabetic subjects. Med Biol Eng Comput 2005;43:252257. [ Links ]
14. La Rovere MT, Bigger JT Jr., Marcus FI, et al. Baroreflex sensitivity and heartrate variability in prediction of total cardiac mortality after myocardial infarction. ATRAMI (Autonomic Tone and Reflexes After Myocardial Infarction) Investigators. Lancet 1998;351:478484. [ Links ]
15. Chesterton LJ, Selby NM, Burton JO, et al. Categorization of the hemodynamic response to hemodialysis: the importance of baroreflex sensitivity. Hemodial Int 2010;14:1828. [ Links ]
16. Sapoznikov D, Backenroth R, Rubinger D. Baroreflex sensitivity and sympathovagal balance during intradialytic hypotensive episodes. J Hypertens 2010;28:314324. [ Links ]
17. CarrascoSosa S, GaitanGonzalez MJ, GonzalezCamarena R, et al. Baroreflex sensitivity assessment and heart rate variability: relation to maneuver and technique. Eur J Appl Physiol 2005;95:265275. [ Links ]
18. Mancia G, Grassi G, Bertinieri G, et al. Arterial baroreceptor control of blood pressure in man. J Auton Nerv Syst 1984;11:115124. [ Links ]
19. Fadel PJ. Arterial baroreflex control of the peripheral vasculature in humans: rest and exercise. Med Sci Sports Exerc 2008;40:20552062. [ Links ]
20. Convertino VA. Effects of exercise and inactivity on intravascular volume and cardiovascular control mechanisms. Acta Astrounaut 1992;27:123129. [ Links ]
21. Langewouters GJ, Settels JJ, Roelandt R, et al. Why use Finapres or Portapres rather than intraarterial or intermittent noninvasive techniques of blood pressure measurement? J Med Eng Technol 1998;22:3743. [ Links ]
22. Infante O, Cortes S. Relación entre el pulso fotopletismográfico digital y la presión arterial incruenta latido a latido ante el reto ortostático. IFMBE Proceedings 2007;18:5155. [ Links ]
23. Allen J. Photoplethysmography and its application in clinical physiological measurement. Physiol Meas 2007;38:139. [ Links ]
24. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Heart rate variability. Standards of measurement, physiological interpretation, and clinical use. Eur Heart J 1996;17:354381. [ Links ]
25. Infante O, Valenzuela S, Polo S. Algoritmo que utiliza la segunda derivada para identificar el complejo QRS en tiempo real. Rev Mex Ing Biomed 1992;13:2332. [ Links ]
26. Malpas SC. Neural influences on cardiovascular variability: possibilities and pitfalls. Am J Physiol 2002;282:H6H20. [ Links ]
27. VargasSabadías A. Regresión y correlación. En: Servicio de publicaciones de la Universidad de Castilla, ed. Estadística descriptiva e inferential. 2 ed. Murcia, La Mancha. COMPOBELL S.L. 1995. 183218. [ Links ]
28. Welch PD. The use of fast Fourier transform for the estimation of power spectra: a method base on time averagingover short, modified periodograms. IEEE Trans Audio and Electroaccoust 1967;15:7073. [ Links ]
29. Parati G, Castiglioni P, Di Rienzo M, et al. Sequential spectral analysis of 24hour blood pressure and pulse interval in humans. Hypertension 1990;16:414421. [ Links ]
30. Kardor A, Rudas L, Simon J, et al. Effect of postural changes on arterial baroreflex sensitivity assessed by the spontaneous sequence method and Valsalva manoeuvre in healthy subjects. Clin Auton Res 1997;7:143148. [ Links ]
31. Foo JY, Lim CS, Wilson SJ. Photoplethysmographic assessment of hemodynamic variations using pulsatile tissue blood volume. Angiology 2008;59:745752. [ Links ]
32. Hermosillo AG, Marquez MF, JaureguiRenaud K, et al. Orthostatic hypotension, 2001. Cardiol Rev 2001;9:339347. [ Links ]
33. Middleton PM, Chan GS, O'Lone E, et al. Changes in left ventricular ejection time and pulse transit time derived from finger photoplethysmogram and electrocardiogram during moderate haemorrhage. Clin Physiol Funct Imaging 2009;29:163169. [ Links ]
34. Guerrisi M, Vannucci I, Toschi N. Differential response of peripheral arterial compliancerelated indices to a vasoconstrictive stimulus. Physiol Meas 2009;30:81100. [ Links ]
35. Palada I, Eterovic D, Obad A, et al. Spleen and cardiovascular function during short apneas in divers. J Appl Physiol 2007;103:19581963. [ Links ]
36. Raamat R, Jagomagi K, Talts J. Calibrated photoplethysmographic estimation of digital pulse volume and arterial compliance. Clin Physiol Funct Imaging 2007;27:354362. [ Links ]
37. Babchenko A, Davidson E, Ginosar Y, et al. Photoplethysmographic measurement of changes in total and pulsatile tissue blood volume, following sympathetic blockade. Physiol Meas 2001;22:389396. [ Links ]
38. Nitzan M, Babchenko A, Shemesh D, et al. Influence of thoracic sympathectomy on cardiac induced oscillations in tissue blood volume. Med Biol Eng Comput 2001;39:579583. [ Links ]