Optical oximetry of volume-oscillating vascular compartments: contributions from oscillatory blood flow

Abstract. We present a quantitative analysis of dynamic diffuse optical measurements to obtain oxygen saturation of hemoglobin in volume oscillating compartments. We used a phasor representation of oscillatory hemodynamics at the heart rate and respiration frequency to separate the oscillations of tissue concentrations of oxyhemoglobin (O) and deoxyhemoglobin (D) into components due to blood volume (subscript V) and blood flow (subscript F): O=OV+OF, D=DV+DF. This is achieved by setting the phase angle Arg(OF)−Arg(O), which can be estimated by a hemodynamic model that we recently developed. We found this angle to be −72  deg for the cardiac pulsation at 1 Hz, and −7  deg for paced breathing at 0.1 Hz. Setting this angle, we can obtain the oxygen saturation of hemoglobin of the volume-oscillating vascular compartment, SV=|OV|/(|OV|+|DV|). We demonstrate this approach with cerebral near-infrared spectroscopy measurements on healthy volunteers at rest (n=4) and during 0.1 Hz paced breathing (n=3) with a 24-channel system. Rest data at the cardiac frequency were used to calculate the arterial saturation, S(a); over all subjects and channels, we found 〈SV〉=〈S(a)〉=0.96±0.02. In the case of paced breathing, we found 〈SV〉=0.66±0.14, which reflects venous-dominated hemodynamics at the respiratory frequency.


Introduction 1.Arterial Saturation Measurements by Pulse
Oximetry: From the Finger to the Brain Pulse oximetry is a technique that is based on the optical measurement of blood volume changes (photoplethysmography) to measure the arterial saturation, S ðaÞ , in tissue noninvasively.It is the most used optical technique in clinical routines since the mid-1980s. 1,2The basis of pulse oximetry is to measure the optical spectral signature of the arterial pulsatile signal at the cardiac frequency.Typically, two wavelengths of light, one in the red and one in the near-infrared region, are used for illuminating an extremity, such as the fingertip, and measuring the transmitted light intensity.Every time the heart is pumping blood, arteries are expanding in diameter, hence arterial blood volume increases.Such an increase in blood volume results in a greater hemoglobin concentration in tissue and a higher absorption, hence a reduction in the detected light from the tissue.Since oxyhemoglobin and deoxyhemoglobin have distinguishably different wavelength-dependent absorption coefficients, the ratio of the pulsatile components of the measured intensities at the two wavelengths can be translated into oxygen saturation of hemoglobin in the arteries, i.e., the arterial saturation.The pulsatile component is a measure of the amplitude of intensity oscillations at the heartbeat frequency.Calculating S ðaÞ from pulsatile intensities at two wavelengths can be achieved by applying the modified Beer-Lambert law, 3 which yields the pulsatile components of oxy-and deoxyhemoglobin concentrations, ΔOðtÞ and ΔDðtÞ, respectively.However, the ratio of the differential path length factors (DPFs) at the two wavelengths needs to be known. 3In the absence of this information, commercial pulse oximeters are based on an empirical calibration. 3,4This is typically applied to healthy volunteers by calibrating optical measurements with actual S ðaÞ values obtained on drawn arterial blood, in protocols that involve altering S ðaÞ over a safe range by modulating the fraction of inspired oxygen (FiO 2 ).A key assumption in order to calculate the arterial saturation through the procedure described above is that the pulsatile intensity component originates only from the arterial blood volume.The error associated with finger pulse oximetry systems is typically reported by manufacturers to be 2%, evaluated by the standard deviation of the differences between pulse oximetry signals and actual concurrent measurements of S ðaÞ by drawing blood in healthy subjects. 4However, a standard deviation of 2% reflects an expected error of 4% (two standard deviations) or more in 5% of the examinations, which also corresponds to an error of 3% to 4% reported in clinical studies and even greater for S ðaÞ < 80%. 5,6eflectance-based pulse oximeters have also been proposed. 7,8Using a reflectance geometry on the subject's forehead has a considerable advantage over finger or toe sensors, since the extremities may be poorly perfused, which reduces the pulsatile signal to be measured.However, the pulsatile waveform measured on the forehead with a source-detector distance of about 1 cm has been found to be less reliable, or rather poorly understood, and readings of S ðaÞ tend to be underestimated. 3,9ne possible proposed explanation is that the sensor actually picks up a contribution from a venous pulsation, so that the resulting readings reflect a mixture of arterial and venous saturation. 9,10However, the reasoning for contributing venous readings as source of error was only based on the fact that the readings on the forehead have been underestimating S ðaÞ .It has further been found that the readings on the forehead improved when applying a mild pressure to the forehead sensor. 9,10It shall be pointed out though that venous pulsation has not been measured directly.Furthermore, contributions of venous pulsatile volume change to the measured optical signals can be considered negligible in comparison to arterial volume changes due to the fact that pulsatile arterial pressure is significantly greater than any pulsatile component in the central venous pressure.It is more likely, as argued in this article, that inaccurate readings of arterial saturation from pulsatile optical signals are due to pulsatile blood flow contributions, rather than pulsatile venous contributions.
While forehead sensors with small source-detector separations seem to eliminate possible errors by applying pressure on the sensor, cerebral measurements with large source detector distances (∼3 cm) cannot benefit from such maneuver.Measuring the cerebral concentration and oxygen saturation of hemoglobin at relatively large source-detector distances and two optical wavelengths is accomplished by the noninvasive technique of near-infrared spectroscopy (NIRS), which typically is used to report changes in oxy-and deoxyhemoglobin as a function of time.
While the contribution from a pulsatile venous blood volume change is questionable, it is known that capillary blood flow has a pulsatile component 11,12 that may also affect the measured signals.The oxygen saturation of hemoglobin in blood changes as the blood flows through the smaller arterioles and capillary bed due to oxygen diffusion to the surrounding tissue.Consequently, a pulsatile capillary blood flow may influence the reading of S ðaÞ by pulse oximetry, due to the flow-related change in the tissue concentrations of oxy-and deoxyhemoglobin.A contribution of blood flow changes to the measured optical signals would result in oscillations of oxy-and deoxyhemoglobin, which are out of phase with each other 13 and hence would influence the reading of S ðaÞ .We have repeatedly shown that oxy-and deoxyhemoglobin measurements on muscle tissue are in phase 14 at the respiratory frequency, indicating that blood volume changes are dominant and blood flow changes can be considered negligible.However, cerebral measurements have shown a consistent out of phase behavior between oxy-and deoxyhemoglobin, [14][15][16][17][18][19] indicating that blood flow changes are not a negligible contributor to the optical signals.The reason for this is that blood volume changes are smaller in the brain than other tissues due to the rigid enclosure of the skull, 20 damping the magnitude of vessel diameter changes.Since the amplitude of the pulsatile blood volume changes is smaller, it is to be expected that the influence of pulsatile blood flow on S ðaÞ measurements by pulse oximetry may be greater in the brain than other tissue locations, such as the finger, toe, or earlobe.As stated above, the underlying assumption to use pulsatile intensity changes for arterial saturation measurements is that the pulsatile intensity originates only from arterial blood volume. 21Any additional source of pulsatile intensity, be it from venous blood or pulsatile flow, invalidates the theory and negatively impacts the accuracy of arterial saturation measurements.

Venous Saturation Measurements by
Spiroximetry or Venous-Pooling-Based Methods Using the same idea that forms the basis of pulse oximetry, i.e., to measure the oxygen saturation of hemoglobin of a specific vascular compartment by exploiting blood volume changes specific to that compartment, NIRS methods for measurements of venous saturation, S ðvÞ , have also been proposed.Microvascular S ðvÞ measurements are especially important for cerebral monitoring, since they provide indications about oxygen extraction (the oxygen extraction fraction) and cerebral metabolic demand.Instead of using pulsations at the cardiac frequency, which predominantly impact the arterial blood volume, it has been proposed that respiration-induced changes in the blood pressure and intrathoracic pressure mostly affect the venous compartment due to the much larger vascular compliance in veins with respect to arteries.3][24][25][26] This method, sometimes called spiroximetry, has been applied in animal models and in human skeletal muscle, 22 and also to cerebral measurements in healthy subjects 24,25 and ventilated patients. 23,24,26piroximetry is based on taking the ratio of amplitudes of oscillating signals (the envelope).Similar to pulse oximetry, spiroximetry is therefore only valid if volume changes are the sole source of signal changes.This implies that oxy-and deoxyhemoglobin are in phase with each other.It further implies and will be demonstrated here that blood flow changes at the respiration frequency will invalidate the assumptions made by spiroximetry.Measurements of oscillatory hemodynamics at the respiration rate are practically appealing for obtaining S ðvÞ since they rely on spontaneous respiration and do not require any external maneuvers.In addition to using oscillatory, respiration-induced hemodynamic changes, several other methods for eliciting venous volume changes have been proposed.For example, one method for cerebral measurements is based on a head-down tilting protocol; tilting the head down by 15 deg results in venous pooling and hence a venous volume change. 27In this approach, the venous saturation is calculated from the difference in the cerebral concentrations of oxy-and deoxyhemoglobin before and after head tilting.Another example, for measurements on fingers, involves eliciting venous volume changes by venous occlusion. 4,28Again, the difference in hemoglobin concentration before occlusion and during occlusion is used for the measurement of S ðvÞ .While these latter measurements have shown promise in body parts where venous occlusions are possible, applications to cerebral venous saturation are limited.For brain measurements, the head-tilting method has its own limitations, which includes some reliability issues in the induced hemodynamic changes, 27 a possible arterial contribution, and potential optical path length changes due to head tilting.
Similar to the case of pulse oximetry, these methods for venous saturation measurement also rely on the same assumption that a specific vascular compartment (venous in this case) dominates the blood volume changes that contribute to the optical signals.

Accounting for Contributions from Blood Flow Oscillations
While venous blood volume oscillations are indeed a major contributor to the optical signals at the respiration frequency, we have recently demonstrated that blood flow oscillations may also play a crucial role and need to be accounted for. 13,16,29his result follows from a hemodynamic model that we have recently introduced, which translates oscillatory changes in cerebral blood volume (CBV) and cerebral blood flow (CBF) into changes in the tissue concentrations of oxyhemoglobin, ΔO, and deoxyhemoglobin, ΔD.This model, combined with inducing hemodynamic oscillations at a set of controlled frequencies, has led to the new technique of coherent hemodynamics spectroscopy, 13,14,16,18 which aims at a quantitative assessment of the microcirculation.
In our previous applications of the hemodynamic model, which is an analytical model that takes into account dynamic autoregulation and the microvascular blood transit times, we have focused on low-frequency cerebral hemodynamics (of the order of 0.1 Hz or less) [13][14][15][16][17][18]29 for the assessment of blood flow, autoregulation, and blood volume in the microvasculature. Inthis work, we demonstrate that the pulsatile signals measured in the brain at the heart rate and at the respiration frequency can be quantitatively described by the model and that CBF changes should be taken into account to obtain reliable measurements of the oxygen saturation of hemoglobin of timevarying vascular compartments.We propose a method, based on our hemodynamic model, which allows for the separation of the contributions of pulsatile blood volume and blood flow to the measured oscillations of hemoglobin concentration.This method yields measurements of the oxygen saturation of hemoglobin of the volume-oscillating vascular compartment, be it arterial, venous, or a combination of both, that take into account the potentially confounding contributions from oscillatory blood flow.

Analytical Hemodynamic Model
We have recently introduced a hemodynamic model, which describes time-dependent expressions of the absolute tissue concentrations of OðtÞ, DðtÞ, and TðtÞ (with units of micromoles per liter of tissue) as a function of dynamic changes of CBV and CBF normalized to baseline (cbvðtÞ ¼ ΔCBVðtÞ∕CBV 0 , cbfðtÞ ¼ ΔCBFðtÞ∕CBF 0 ).Dynamic changes in the metabolic rate of oxygen, cmro 2 , are also considered by the model but they are neglected here because we are considering cases where cerebral hemodynamics is dominated by blood volume and blood flow changes.In general, the time evolution of hemoglobin concentrations in tissue depends on both cbv and cbf (again, we neglect cmro 2 here).The measurable, dynamic and hence time-dependent quantities of OðtÞ and DðtÞ can be written as the sum of dynamic blood flow and blood volume contributions, indicated with subscripts F and V, respectively: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 1 ; 6 3 ; 1 4 8 OðtÞ ¼ O F ðtÞ þ O V ðtÞ; (1) E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 2 ; 6 3 ; 1 1 9 DðtÞ ¼ D F ðtÞ þ D V ðtÞ: The subscripts "V" and "F" refer to dynamic changes in OðtÞ and DðtÞ resulting from dynamic changes in volume and flow, hence due to cbvðtÞ and cbfðtÞ.The hemodynamic model describes explicitly how given cbfðtÞ and cbvðtÞ yield O F ðtÞ, O V ðtÞ, D F ðtÞ, and D V ðtÞ.Specifically, the model derives temporal dynamics (time-domain description) of hemoglobin changes, as well as sinusoidal hemodynamic oscillations as a function of the angular frequency ω (frequency-domain description).Here, we adopt the frequency-domain version of the model for describing oscillating hemodynamic signals such as those related to cardiac pulsation and respiration.We adopted the phasor representation of oscillatory quantities, where phasors are indicated in boldface.By definition, phasors (boldface) are dynamic quantities.The model expressions for OðωÞ, DðωÞ, TðωÞ (i.e., the phasors that describe the oscillations of the oxy-, deoxy-, and total hemoglobin concentrations, respectively) as a function of cbvðωÞ and cbfðωÞ (i.e., the phasors that describe the oscillations of CBV and CBF) are as follows: 13,16 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 3 ; 3 2 6 ; 5 8 7 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 6 ; 3 2 6 ; 4 4 3 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 7 ; 3 correspond to the baseline contributions of blood volume in the arterial, capillary, and venous compartment, respectively.Boldface lowercase, quantities correspond to dynamic changes relative to baseline values.1][32][33][34][35] Note that the only measurable parameters in Eqs. ( 1)-( 9) are jOðωÞj, jDðωÞj, jTðωÞj, ArgðDÞ − ArgðOÞ, and ArgðOÞ − ArgðTÞ.Finally, we observe that one can define an overall blood volume phasor, cbvðωÞ, as follows: CBV 0 cbv ðvÞ ðωÞ; (10)   so that, as seen in Eq. ( 9), the phasor of total hemoglobin concentration, TðωÞ, has the same phase as cbvðωÞ.

Derivation of the Basic System of Phasor Equations That Describes Cerebral Hemodynamics
We now make the key assumption that arterial and venous blood volumes oscillate in phase, which should be the case when blood volume oscillations are driven by systemic blood pressure changes.In other words, the phasors, cbv ðaÞ ðωÞ, cbv ðvÞ ðωÞ, and cbvðωÞ, have the same phase.Under these conditions, Eq. ( 10) directly translates into a condition for the magnitudes of the blood volume phasors: One may introduce a factor, ρ ðvÞ with 0 ≤ ρ ðvÞ ≤ 1, that specifies the venous fraction of the total oscillating blood volume as follows: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 2 ; 6 3 ; 3 9 2 ρ ðvÞ ¼ CBV ðvÞ 0 jcbv ðvÞ j CBV 0 jcbvj ; from which the arterial fraction of the total oscillating blood volume is ; t e m p : i n t r a l i n k -; e 0 1 3 ; 6 3 ; 3 2 3 ρ ðaÞ ¼ 1 − ρ ðvÞ ¼ CBV ðaÞ 0 jcbv ðaÞ j CBV 0 jcbvj : We point out that ρ ðvÞ and ρ ðaÞ depend on the nature of the oscillatory hemodynamics, and, in general, they also depend on ω.The condition of Eq. ( 11) and the definition of ρ ðvÞ in Eq. ( 12) allow us to express the blood volume components of the oxyand deoxyhemoglobin phasors of Eqs. ( 4) and ( 7) as follows: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 5 ; 6 3 ; 1 7 4 D V ðωÞ ¼ ctHb CBV 0 ½1 − ð1 − ρ ðvÞ ÞS ðaÞ − ρ ðvÞ S ðvÞ cbvðωÞ: Finally, by introducing the oxygen saturation of hemoglobin of the volume oscillating compartment S V ¼ jO V j∕jTj, one finds that it can be written as a weighted average of the arterial saturation (S ðaÞ ¼ jO ðaÞ V j∕jT ðaÞ j) and venous saturation (S ðvÞ ¼ jO ðvÞ V j∕jT ðvÞ j), with the weights given by the respective arterial and venous fractions of oscillating blood volume: ; t e m p : i n t r a l i n k -; e 0 1 6 ; 3 2 6 ; 7 5 2 With these definitions, the model equations for O V and D V become E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 7 ; 3 2 6 ; 6 7 4 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 8 ; 3 2 6 ; 6 4 4 We are now in a position to write the system of phasor equations that describes cerebral hemodynamics driven by blood volume and blood flow oscillations: ; t e m p : i n t r a l i n k -; e 0 1 9 ; 3 2 6 ; 5 8 5

> > <
> > : with the condition 0 ≤ S V ≤ 1.The first two equations in Eq. ( 19) repeat Eqs. ( 3) and ( 6), the third equation follows directly from Eqs. ( 17) and ( 18), and the fourth equation follows directly from Eqs. ( 5) and (8).By recalling that phasors are twodimensional vectors, Eq. ( 19) is a linear system of eight scalar equations in nine unknowns (the components of the phasors O V , O F , D V , D F and S V ), so that one needs to add one more condition to fully solve it.Two relevant scenarios are presented in the next two sections.

Setting the Oxygen Saturation of Hemoglobin of the Volume-Oscillating Vascular Compartment
One approach to fully solve the system of Eq. ( 19) is to set the value of S V , the oxygen saturation of hemoglobin of the volume oscillating vascular compartment.Figure 1 shows a representative phasor diagram for the general case in which blood flow oscillations significantly contribute to the hemoglobin concentration oscillations.Here, we take the phase of O as the phase reference [so that ArgðOÞ is zero by definition].Based on phasor algebra or a geometrical analysis of Fig. 1, one can derive the following solution to the linear system of Eq. ( 19): E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 2 ; 3 2 6 ; 2 2 2 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 3 ; 3 2 6 ; 1 9 7 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 4 ; 3 2 6 ; 1 7 1 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 5 ; 3 2 6 ; 1 4 6 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 6 ; 3 2 6 ; 1 2 0 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 7 ; 3 2 6 ; 9 5 D Fy ¼ S V D y ; where the subscripts x and y indicate the x and y components, respectively, of the corresponding phasor.Note that by taking the phase of O as the phase reference, O only has an x component (O y ¼ 0).By inspecting Fig. 1, one may appreciate that the phasor representation leads to an easy, graphical solution of the system of Eq. (19).In fact, one may start by drawing the phasors O, D, and T, which are the measured phasors in an experiment.Then, based on the assumed value of the oxygen saturation of hemoglobin of the time-varying vascular compartment (S V ), one can draw O V and D V along the direction of T, as given by O and D F can be drawn as the phasors from the tip of O V to the tip of O and from the tip of D V to the tip of D, respectively.

Setting the Angle Between O and O F
An alternative approach to setting S V , the case described in the previous section, is to set the angle between the O and O F phasors.To appreciate the meaning of this assumption, one needs to consider the blood flow phasor (cbf) and its phase relationship with the corresponding oxyhemoglobin concentration phasor (O F ).
Let us first consider the phasor cbf.As mentioned above, Eqs. ( 9) and (10) show that the measured total hemoglobin phasor T is in phase with the blood volume phasor cbv.Because of the high-pass nature of the cerebral autoregulation process that regulates CBF in response to blood pressure changes, [36][37][38] we consider the following high-pass relationship between cbf and cbv: 13 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 8 ; 6 3 ; 2 6 9 cbfðωÞ ¼ kH ðARÞ HP ðωÞcbvðωÞ; ( where k is the inverse of the modified Grubb's exponent, and H ðARÞ HP ðωÞ is the RC high-pass transfer function with cutoff frequency f ðARÞ c that describes the effect of autoregulation (f ðARÞ c ∼ 0.03 − 0.05 Hz).Equation (28) quantitatively describes the relationship between the cbf and cbv phasors in terms of a simple high-pass transfer function.While the specific quantitative aspects of this relationship may be improved by more sophisticated models of autoregulation (possibly considering also the relationship between blood volume and blood pressure, the latter being the quantity directly involved with autoregulation), the qualitative relationship between cbf and cbv is correctly described by Eq. (28).Specifically, at higher frequencies (say, at the heartbeat frequency), cbf and cbv are in phase, whereas at lower frequencies (say, at the frequency of normal breathing), cbf leads cbv.The specific angle by which cbf leads cbv depends on the autoregulation model used, and here we adopt the model of Eq. ( 28) based on an RC highpass transfer function with f ðARÞ c ¼ 0.03 Hz.Let us now consider the phasor O F , which describes the oxyhemoglobin concentration oscillations that are induced by blood flow oscillations.O F is not in phase with cbf because the effects of a change in blood flow on the concentration of oxyhemoglobin are mediated by the blood transit time in the microcirculation.This effect is quantitatively described by Eq. ( 5) of the hemodynamic model, where the low-pass transfer functions associated with capillary (H ðcÞ LP ) and venous (H ðvÞ LP ) flow account for the fact that O F lags cbf.The characteristic times of the capillary and venous transfer functions are on the order of t ðcÞ ∕e ¼ 0.37t ðcÞ and 0.28ðt ðcÞ þ t ðvÞ Þ, respectively. 13In this work, we use Eqs.( 5) and (28) to estimate the phase of O F (and therefore, the angle between O and O F ), but other methods may also be used on the basis of different hemodynamic and autoregulation models.
By setting the phase angle of O F with respect to O (we refer to this angle as ½ArgðO F Þ − ArgðOÞ ¼ ArgðO F Þ), one can immediately find an expression for S V .In fact, by definition, ðO Fy ∕O Fx Þ ¼ tan½ArgðO F Þ, and after replacing O Fx and O Fy with their expressions given in Eqs. ( 22) and ( 23), one finds E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 9 ; 3 2 6 ; 4 9 5 By replacing S V in Eqs. ( 20)-( 27) with the expression of Eq. ( 29), one finds the full solution to the system of Eq. (19).As in the case of the previous section, the phasor diagram of Fig. 1 can be used to yield a graphical solution of the system of Eq. (19).In fact, after assuming the angle between O and O F , i.e., the direction of the oxyhemoglobin phasor associated with blood flow oscillations, one finds the tips of O V and D V (and therefore, S V ¼ jO V j∕jTj) by intersecting the T phasor with the line parallel to O F that passes through the tip of O and D, respectively.
If the dynamics of the tissue concentrations of oxy-and deoxyhemoglobin are dominated by blood volume changes, with negligible contributions from blood flow changes, O and D are in phase with each other and only have x components, so that D y ¼ 0. Therefore, Eq. (29) shows that in this case, the oxygen saturation of hemoglobin of the volume oscillating compartment, S V , is simply given by the amplitude ratio of the oxyhemoglobin to total hemoglobin oscillations: One can see that Eq. ( 33) is in agreement with Eq. ( 30), where S V becomes S ðaÞ in this case where the volume-oscillating vascular compartment is the arterial one.If blood volume changes are the only source of signal contributions, oxy-and deoxyhemoglobin changes occur simultaneously, so that the phase difference Arg½D − Arg½O ¼ 0. Figure 2 illustrates the phasor diagram for this scenario of negligible blood flow contributions and blood volume contributions from solely arterial blood.
If, instead, blood flow oscillations yield non-negligible contributions to the measured hemoglobin concentrations, Eqs. ( 31)- (33) are not valid and one must take O F and D F , and hence, the cbf terms in the model equations [Eqs.( 5) and ( 8)], into account.Since cbfðωÞ introduces a phase shift (which is frequency dependent) between O and D, the phase difference Arg½D − Arg½O is no longer 0. From this fact, it follows that jO V j ≠ jOj, and therefore (case of non-negligible blood flow contributions), Equations ( 34) and ( 35) combine Eqs.(3)-( 5) and ( 6)-( 8), respectively, and specify that the oscillations are at the cardiac frequency (ω hr ).
In this case of non-negligible blood flow contributions: When O and D are not in phase with each other, it is important to point out that E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 3 7 ; 3 2 6 ; 6 7 2 In summary, measuring the arterial saturation on the basis of the rightmost term in Eq. ( 33) is only correct if blood flow contributions to the measured oscillations of oxy-and deoxyhemoglobin are negligible.

Measuring the Venous Saturation from Volume
3][24][25][26] The idea is that at the respiration frequency, venous blood volume oscillations may be the dominant source of oxy-and deoxyhemoglobin concentration dynamics.Using these hemodynamic oscillations for calculating S ðvÞ has led to a method of venous saturation measurements [16][17][18][19] for which the term spiroximetry was coined. 22As shown above for the arterial saturation, if blood flow contributions are negligible, then Arg½D − Arg½O ¼ 0. However, we have previously demonstrated, [13][14][15][16]29 and we will show again here that this is not the case for brain measurements, where blood flow contributions cannot be neglected. Hece, it is not correct to use the ratio jOj∕jTj to obtain S ðvÞ from data collected during paced or normal breathing, since O ≠ O V .While the value of arterial saturation S ðaÞ can often be reliably assumed since it is a robust systemic parameter that takes values in a narrow range for healthy individuals (leading to the approach of Sec 2.3), the value of venous saturation depends on the specific, local conditions of oxygen supply (through blood flow) and oxygen delivery (to accommodate metabolic demand) of the tissue considered.Therefore, the measurement of oxygen saturation of hemoglobin in a volume-oscillating compartment that comprises venous blood can be achieved according to the approach of Sec.2.4, where one assumes the phase of the flow component of the oxyhemoglobin concentration phasor (O F ) with respect to the total, measured phasor O (i.e., one assumes the value of ½ArgðO F Þ − ArgðOÞ).
In general, one should also consider the fact that arterial and venous blood may both contribute to the blood volume oscillations associated with the total hemoglobin concentration phasor: T ¼ T ðaÞ þ T ðvÞ .The relative venous and arterial contributions to the oscillatory blood volume are quantified by the factors ρ ðvÞ and ρ ðaÞ ¼ 1 − ρ ðvÞ introduced in Sec.2.2.If ρ ðvÞ ¼ 1, i.e., only venous blood contributes to T, then the oxygen saturation of hemoglobin of the volume-oscillating compartment (S V ) coincides with the venous saturation (S ðvÞ ).Otherwise, in general, one needs to estimate ρ ðvÞ to assess the contributions of arterial and venous saturation to S V .

Human Subjects
Ten healthy subjects were recruited for the study (six males and four females, age range: 25 to 49 years) without any history of neurological disorders or cardiovascular disease.The experimental protocol was approved by the Tufts University Institutional Review Board (IRB).Written informed consent was obtained from all participants prior to the study.Subjects were seated in a darkened room and instructed to breathe normally.The experimental protocol consisted of acquiring 10 min of resting baseline measurements.In addition, three of the subjects performed 6.5 min of paced breathing at a frequency of 0.1 Hz, as guided by a metronome.

Near-Infrared Spectroscopy Data Collection and Data Processing
Optical data were collected with a commercial NIRS optical topography system (ETG 4000, Hitachi Medical Corporation, Japan), operating at wavelengths of 700 and 830 nm.Optical probes were placed on the subject's head, in the area corresponding to the frontal lobe, and secured by a dedicated cap to hold them in place.The probe consisted of eight illumination fibers and eight collection fibers, arranged in a grid, as seen in Fig. 3.The source-detector distance for all measurement channels was 3 cm; there was a total of 24 source-detector pairs (or measurement channels), with each channel acquiring data with a 10-Hz sampling frequency.After removing slow temporal drifts from the NIRS data, optical intensity changes were translated into relative changes of ΔOðtÞ, ΔDðtÞ, and ΔTðtÞ (oxy-, deoxy, and total hemoglobin concentrations, respectively) by applying the modified Beer-Lambert law to each channel.A DPF ratio of DPF 700 ∕DPF 830 ¼ 1 (specifically: DPF 700 ¼ DPF 830 ¼ 6.51) was used unless otherwise stated.
The fast Fourier transform (FFT) was taken for each channel's intensity time trace, normalized to its baseline value (ΔIðtÞ∕I 0 , where ΔIðtÞ ¼ IðtÞ − I 0 ) for both wavelengths.The baseline value I 0 is defined as the mean value over the entire time trace.Only those channels for which the magnitude of the Fourier transform at the cardiac frequency (for the baseline measurements) or at 0.1 Hz (for the paced breathing measurements) was at least 5 times greater than the average magnitude of the Fourier transform at frequencies >4 Hz have been considered for further processing.The other channels were discarded because considered insensitive to cerebral hemodynamics and dominated by noise.In addition, only those subjects with at least 20 channels meeting this criterion were included for further processing.Four out of 10 subjects satisfied this requirement for baseline measurements, and 3 out of 3 subjects for paced breathing measurements.
In order to study the phase relationship, as well as the amplitudes, of the pulsatile signals, relative intensities at each wavelength, ΔI 700 ðtÞ∕I 0 700 and ΔI 830 ðtÞ∕I 0 830 , as well as ΔOðtÞ, ΔDðtÞ, and ΔTðtÞ, were band-pass filtered around the cardiac frequency, which was in the range 0.96 to 1.2 Hz in this study (filter width: AE0.2 Hz; filter order: n ¼ 214) for the baseline dataset, and around 0.1 Hz (filter width: AE0.02 Hz; filter order: n ¼ 536) for the paced breathing data with a linearphase band-pass filter (function "firpmord" in MATLAB, Mathworks, Natick, Massachusetts).A greater filter width was chosen for filtering at the cardiac frequency to accommodate possible changes in heart rate over time.The instantaneous amplitudes and phases of these signals have been obtained by applying the Hilbert transform to generate the complex analytic signal.Average phase differences, ArgðDÞ − ArgðOÞ and Argði 700 Þ − Argði 830 Þ, where boldface indicates the phasor representation of the band-pass filtered ΔOðtÞ, ΔDðtÞ, i 700 ¼ ΔI 700 ðtÞ∕I 0 700 , and i 830 ¼ ΔI 830 ðtÞ∕I 0 830 , and standard deviations were computed by means of circular statistics.In order to avoid any noise or instability introduced by the filter, the first and last 60 s of the filtered signals have been removed before computing averages and standard deviations.To guarantee measurements of regular hemodynamic oscillations in tissue rather than incoherent physiological fluctuations or motionrelated artifacts, we further applied the following selection criteria to the data: i.We only considered those signals where jDj > 0.015 μM and jOj > 0.015 μM, which are amplitude threshold values reported previously to discard frequency components of instrumental noise. 18i.We retained only those signals with a stable phase difference, namely σ ArgðDÞ−ArgðOÞ < 25 deg, with σ denoting the circular standard deviation, to measure coherent hemodynamic oscillations.iii.We discarded those channels that featured signals dominated by motion artifacts (synchronous with the heartbeat, respiration, or other physiological rhythms) or featuring a poor optical contact between the fiber optics and the scalp.

Cerebral Arterial Saturation from Hemodynamic Oscillations at the Cardiac Frequency
Out of the 10 subjects measured, four subjects met the criteria of having at least 20 channels with a signal-to-noise ratio ðSNRÞ > 5 at the cardiac frequency, as determined by the FFT of ΔIðtÞ∕I 0 , as described in Sec.2.8.Of those four subjects, three performed paced breathing at 0.1 Hz (subject nos. 1, 2, and 4).A representative case of filtered ΔI 700 ðtÞ∕I 0 700 and ΔI 830 ðtÞ∕I 0 830 at the cardiac frequency can be seen in Fig. 4(a) for subject no. 1, channel 18.We introduce a phasor notation for the optical intensity (i) that follows the same conventions used for the cbv and cbf phasors.The lower case intensity phasor, i, indicates the sinusoidal intensity oscillations relative to the mean value: ΔIðtÞ∕I 0 .An average phase difference of Argði 700 Þ − Argði 830 Þ ¼ 10 deg was found for this case (subject no. 1, channel 18).After converting the relative intensity changes into ΔOðtÞ, ΔDðtÞ, and ΔTðtÞ by using the modified Beer-Lambert law and a DPF ratio of 1 (i.e., the same DPF ¼ 6.51 at 700 and 830 nm), we found a phase difference ArgðDÞ−ArgðOÞ¼−269 deg, as seen in Fig. 4(b).Figures 4(c) and 4(d) show the average phase differences of the intensities at the two wavelengths and of the deoxy-and oxyhemoglobin concentrations, respectively, for all the channels.Error bars show the standard deviation over 10 min of instantaneous phase difference measurements.The black circle in Fig. 4(c) marks the channels with a phase difference that is significantly different from zero.Since oxy-and deoxyhemoglobin oscillations, as shown in Fig. 4(b), are not in phase with each other, ðjOj∕jTjÞ ≠ ½jOj∕ðjOj þ jDjÞ.Therefore, calculating S ðaÞ as jOj∕jTj or jOj∕ðjOj þ jDjÞ yields two different results.Such a difference is seen in Fig. 5 for subject no. 1 and for all measurement channels.The ratio jOj∕jTj is shown by the black circles and often takes values of S ðaÞ > 1, which is of course unreasonable.The ratio jOj∕ðjOj þ jDjÞ is shown by the gray squares and yields values of S ðaÞ < 0.98.
Figure 6 shows the results for channel 18 of subject no. 1 as a phasor diagram.The average phase difference for Arg½D − Arg½O ¼ 91 deg (or −269 deg).Based on Eqs. ( 20)-( 27) and the assumption S ðaÞ ¼ 0.98, we calculated O F and D F , as well as O V and D V .Arg½O F − Arg½O was found to be −67 deg.The phasor diagram clearly shows that the flow component of oxyhemoglobin oscillations is much smaller than the volume component (jO F j ≪ jO V j), while the magnitudes of D F and D V are comparable.For this particular case, jOj∕jTj yields 1.0, which overestimates the set value of arterial saturation (0.98), while jOj∕ðjOj þ jDjÞ yields 0.95, which underestimates it.
While the fact that Arg½D − Arg½O ≠ 0 is already indicative of CBF influences on the measured hemoglobin concentration oscillations, the specific value of this phase difference also depends on the ratio of the DPFs at the two wavelengths, DPF 700 ∕DPF 830 .Since the measurements were performed with a CW domain system, we did not measure the DPFs.We have therefore evaluated the effect of the DPF ratio on Arg½D − Arg½O as well as Arg½O F − Arg½O (data not shown).We found that the DPF ratio would need to extend outside literature values (0.9 to 1.1) 38 in order to result in Arg½D − Arg½O ¼ 0 for all channels.Hence, the DPF ratio cannot account for the phase difference between D and O measured in most channels, but CBF effects do, as we show next.
We have used our hemodynamic model, specifically Eqs. ( 34) and ( 35), to estimate Arg½D − Arg½O and Arg½O F − Arg½O.As described in Eqs. ( 34) and ( 35), we have assumed the venous volume oscillation to be negligible.Furthermore, we have set the frequency of oscillations to ω∕ð2πÞ ¼ 1 Hz to represent a typical cardiac frequency.At a frequency of 1 Hz, which is significantly higher than the cutoff frequency of the RC high-pass filter used to model cerebral autoregulation, the transfer function in Eq. ( 4) becomes a factor of about 1, thus simplifying Eq. ( 28) to E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 3 8 ; 3 2 6 ; 4 4 4 We calculated the phase differences Arg½D − Arg½O and Arg½O F − Arg½O for a physiological range of values of the model parameters.It shall be pointed out that these phase differences are independent of the value of T 0 , the cerebral hemoglobin concentration, which has a typical value of 57.5 μM. 39,40Specifically, the model parameters were iterated within their physiological ranges and histograms of phase differences were obtained.The model parameters, and their considered physiological ranges, 13,16 were capillary blood transit time t ðcÞ ¼ 0.4 to 2 s, venous blood transit time t ðvÞ ¼ 2 to 7 s, inverse of the modified Grubb's exponent k ¼ 2 to 5, rate constant of oxygen diffusion α ¼ 0.2 to 1.4 s −1 , and arterial saturation S ðaÞ ¼ 0.94 to 1.The histograms in Fig. 7 show the computed phase differences using the model equations.The wider spread in the experimental data may be attributed to errors in the phase measurements, especially in channels with lower SNR, and to the set DPF ratio, which may vary across the measurement channels.
Instead of setting S ðaÞ ¼ 0.98 in order to obtain Arg½O F − Arg½O, we also calculated S ðaÞ by setting Arg½O F − Arg½O based on the model prediction, as described in Sec.2.4.For this, we set Arg½O F − Arg½O ¼ −72 deg, which is the average phase difference found by using the model.Using this phase difference, the average arterial saturation over all subjects and channels was found to be S ðaÞ ¼ 0.96 AE 0.02.In comparison, without taking into account blood flow contributions, we found ðjOj∕jTjÞ ¼ 0.97 AE 0.02 and ½jOj∕ðjOj þ jDjÞ ¼ 0.93 AE 0.09.We also considered the spread of the model predictions for Arg½O F − Arg½O, namely Arg½O F − Arg½O ¼ −72 AE 11 deg, resulting in the range from −83 to −61 deg (mean AE standard deviation).Correspondingly, we found S ðaÞ ¼ 0.96 AE 0.02, for Arg½O F − Arg½O ¼ −83 deg and S ðaÞ ¼ 0.96 AE 0.02, for Arg½O F − Arg½O ¼ −61 deg, with errors given by the standard deviation over all channels and subjects.Therefore, the exact value of Arg½O F − Arg½O does not play a crucial role in the determination of the value of S ðaÞ .As described above (Sec.3.1) in the case of data filtered at the cardiac frequency, the model was used to predict Arg½D − Arg½O and Arg½O F − Arg½O for the data filtered at the respiration frequency of 0.1 Hz (as applied for 1 Hz as well).The model parameters were iterated as described above, but setting S ðaÞ ¼ 98%.While the autoregulation high-pass filter did not influence the phase of the hemoglobin concentration oscillations at 1 Hz (which is much higher than the filter cutoff frequency), it does at 0.1 Hz (which is much closer to the filter cutoff frequency).The cutoff frequency of autoregulation has been set to f ðARÞ c ¼ 0.03 Hz, which is a typical value we found previously. 17While venous volume oscillations could be neglected at the cardiac frequency, both arterial and venous volume oscillations have been taken into account for the paced breathing protocol at 0.1 Hz.The reason for taking arterial volume changes into account is that it is known that mean arterial blood pressure (MAP), which is directly linked to arterial blood volume changes, is oscillating with paced breathing. 19Specifically, histograms of the phase differences Arg½D − Arg½O and Arg½O F − Arg½O were generated for three scenarios: (a) T ¼ T ðaÞ , i.e., assuming no venous volume oscillations; (b) T ¼ T ðaÞ þ T ðvÞ , i.e., assuming both arterial and venous contributions to hemoglobin concentration oscillations; (c) T ¼ T ðvÞ , i.e., assuming no arterial volume oscillations.It shall be pointed out that these three scenarios refer to the amplitudes of the oscillations (changes) of the arterial and venous blood volumes, not to the relative contributions of arterial blood versus venous blood.For the case of T ¼ T ðaÞ þ T ðvÞ , we have assumed that ρ ðvÞ ¼ 1∕2; hence, T ðaÞ ¼ T ðvÞ .In order to achieve T ðaÞ ¼ T ðvÞ with the hemodynamic model, we set to the same value the steady state blood volumes of the two vascular compartments, specifically CBV ðaÞ 0 ∕CBV 0 ¼ CBV ðvÞ 0 ∕CBV 0 ¼ 0.2 based on literature, 41 and the oscillatory changes in the blood volumes; hence, cbv ðaÞ ¼ cbv ðvÞ .The corresponding histograms are seen in Fig. 9, with the three colors corresponding to the three different scenarios (a), (b), and (c) defined above.

Cerebral Oxygen Saturation of Hemoglobin of the Vascular Compartment Oscillating at the Respiration Frequency
While Arg½D − Arg½O varies considerably over the different scenarios considered (Fig. 9, left), Arg½O F − Arg½O is strongly peaked between −9 and −4 deg, regardless of which compartment is the source of volume changes and parameters of the model (Fig. 9, right).
We have assumed a value of −7 deg for Arg½O F − Arg½O, which is within the narrow range −9 to −4 deg predicted by the model.Using this value, we have found S V from Eq. ( 29) and phasors O V , O F , D V , and D F from Eqs. ( 20)- (27).Figure 10 shows the solution as a phasor diagram for the case of subject no. 1, channel 18, where one can see the flow and volume contributions to the hemoglobin concentration phasors.The oxygen saturation of hemoglobin of the volume-oscillating vascular compartment was found to be S V ¼ 0.67.The measured angle between D and O was Arg½D − Arg½O ¼ 160 deg (or −200 deg).Based on the model predictions reported in Fig. 9, this phase difference between D and O is close to the peak value for the case T ¼ T ðaÞ þ T ðvÞ (green histogram) and falls on the tails of the distributions for T ¼ T ðaÞ and T ¼ T ðvÞ .This result suggests that both arterial and venous blood may contribute to the volume oscillating compartment, in agreement with the observation that MAP oscillates at the respiration frequency, 19 which may induce arterial volume oscillations in the brain.
Using all 24 channels from the three subjects investigated with the paced breathing protocol, we calculated jOj∕jTj and jOj∕ðjOj þ jDjÞ.The results are reported in the histograms of Fig. 11.The average of jOj∕jTj is 1.1 AE 0.14, and the average of jOj∕ðjOj þ jDjÞ is 0.88 AE 0.05.As discussed above in Sec.2.4 [see Eq. ( 37)], out-of-phase oscillations of oxy-and deoxyhemoglobin concentrations results in jO þ Dj ≠ jOj þ jDj, and hence the difference between these ratios.
By assuming the value of −7 deg for Arg½O F − Arg½O on the basis of model predictions, we calculated the oxygen saturation of hemoglobin of the volume-oscillating vascular compartment, S V for all subjects and channels.The resulting histogram of S V ¼ jO V j∕jTj is reported in Fig. 11 and features a mean and standard deviation of hS V i ¼ 0.66 AE 0.14.Knowledge or an estimation of the factor ρ ðvÞ and the arterial saturation S ðaÞ allows for the determination of the venous saturation S ðvÞ from the volume-oscillating compartment oxygen saturation of hemoglobin S V .For example, assuming ρ ðvÞ ¼ 1∕2, hence T ðaÞ ¼ T ðvÞ , and S ðaÞ ¼ 0.98, the measured value S V ¼ 0.66 AE 0.14 translates into a venous saturation S ðvÞ ¼ 0.35 AE 0.28 [on the basis of Eq. ( 16)].Of course, a greater value of ρ ðvÞ would yield a higher venous saturation.The average phase difference for all channels and subject was found to be Arg½O F − Arg½O ¼ −221 AE 16 deg.This phase difference falls within the green curve in Fig. 9, where T ðaÞ ¼ T ðvÞ .However, in general, it is not expected that the relative proportion of arterial and venous volume is constant across all channels and subjects, so that the assumption of ρ ðvÞ and its constant value   across all subjects and measurement channels is affected by some level of uncertainty, which affects the mean value and variance of S ðvÞ .

Oxygen Saturation of Hemoglobin of Volume-Oscillating Vascular Compartments
Pulse oximetry and spiroximetry are both based on the assumption that the measured oscillatory changes are based on volume changes alone.In case of pulse oximetry, the arterial compartment is dominant, whereas spiroximetry relies on the respiratory signals being dominated by the venous compartment.Both methods are then based on taking the ratio between the magnitudes of oxyhemoglobin oscillations to total hemoglobin oscillations to obtain a measure of the volume oscillating oxygen saturation of hemoglobin. 4,22,26It shall be pointed out that in pulse oximetry technically, the ratio of intensity changes is typically used, which however can be translated into changes in oxy-and deoxyhemoglobin. 4Both methods are only valid if the underlying source of signal changes is based on blood volume changes and if blood flow contributions can be neglected.While this may be the case in tissue other than the brain, 1,4,14,22 neglecting blood flow changes in the brain is questionable.Indeed, studies applying spiroximetry to measure cerebral venous saturation have commented on the fact that measured hemodynamic oscillations are out of phase, 23,25,26 attributing this fact to blood flow changes and invalidating the assumption of spiroximetry.We have shown here that the idea of exploiting blood volume oscillations for measuring S ðaÞ and S ðvÞ by simply taking a ratio of the magnitudes of ΔO∕ΔT, (jOj∕jTj), is not strictly correct and can be at best only an approximation when measuring oscillations of cerebral concentrations of oxy-and deoxyhemoglobin.Even at the cardiac frequency, where arterial volume changes are indeed the dominant source of optical signal changes, CBF oscillations may provide some contributions as well.This is best seen in Figs. 4 and 8, where we have shown that oxy-and deoxyhemoglobin concentrations oscillate out of phase with each other.This phase difference cannot be explained by wrong assumptions in the DPF ratio at the two wavelengths.Of course, the issue is to identify those cases where neglecting the blood flow contributions to the measured optical signals is acceptable and those where it is not: in this article, we have presented a quantitative tool to address this issue.
In the case of hemodynamic oscillations at the cardiac frequency, we have assumed S ðaÞ ¼ 0.98 in order to calculate O V and O F using Eqs.( 20)- (27).However, our hemodynamic model could also be used to first estimate Arg½O F − Arg½O (a phase difference that we have found to be peaked at −72 deg in Fig. 7 for 1 Hz oscillations) and then determine S ðaÞ using Eqs.( 20)- (27) or Eqs. ( 27)- (31).By setting a value of −72 deg for Arg½O F − Arg½O, we have found S ðaÞ ¼ 0.96 AE 0.02.This latter procedure is even more relevant in the case of hemodynamic oscillations at the respiration frequency, as observed in the paced breathing protocol reported here.In this case, the lower paced-breathing frequency of 0.1 Hz (versus ∼1 Hz for the cardiac frequency) results in blood flow oscillations that lead blood volume oscillations, as described by the high-pass relationship of Eq. ( 28), and a smaller phase lag of O F versus O.Such a phase difference, Arg½O F − Arg½O, depends on the assumed cutoff frequency of autoregulation.In this work, we have set the autoregulation cutoff frequency, f ðARÞ c , at 0.03 Hz, and found that Arg½O F − Arg½O is centered at −7 deg (see Fig. 9) for 0.1 Hz oscillations.While this cutoff frequency is a good approximation based on previous results, 17,18  On the basis of these results, considering also that the width of the distribution of Arg½O F − Arg½O does not increase by more than 5 deg over the autoregulation cutoff frequencies considered, we do not expect the uncertainty in f ðARÞ c to significantly impact the results for S V .

Venous Saturation Measurements
3][24][25][26] Based on spiroximetry, as well as low frequency oscillation ratios, Leung et al. 23 found values of venous saturation between 0.55 and 0.71.Lynch et al., 24 using ventilated pulmonary hypertension pediatric patients, reported a comparison to superior vena cava (SVC) saturation.The found values were in good agreement (59.6 for spiroximetry versus 0.64 for SVC).However, using healthy adults and regular respiration, the average value found was 0.79 AE 0.07, which is in good agreement  11 Results of the analysis of the data during paced breathing for all subjects and measurement channels, demonstrating the difference between the ratios of jOj∕jTj, jOj∕ðjOj þ jDjÞ, and reporting the measured oxygen saturation of hemoglobin of the volumeoscillating vascular compartment.The mean values of the measured quantities were jOj∕jTj ¼ 1.1 AE 0.14, jOj∕ðjOj þ jDjÞ ¼ 0.88 AE 0.05, and hS V i ¼ 0.66 AE 0.14.
with our value of hjOj∕ðjOj þ jDjÞi ¼ 0.88 AE 0.05.Both these values are higher than jugular venous saturation, which is typically reported to be between 0.55 and 0.75. 42Wolf et al. 26 also reported cerebral venous saturation based on spiroximetry, with reported values of 0.73 AE 0.09, also greater than expected values.The same group did also mention an out of phase behavior between oxy-and deoxyhemoglobin oscillations, acknowledging the fact that blood flow changes may influence the results.Based on these reported values and results presented here, jOj∕ðjOj þ jDjÞ will overestimate the oxygen saturation of hemoglobin, due to uncorrected blood flow influences.The results presented here are demonstrating that it is possible to correct for such influences by disentangling O V and O F .With this, we have found the volume oscillating compartment oxygen saturation of hemoglobin to be hS V i ¼ 0.66 AE 0.14, which falls within the typically measured oxygen saturation of hemoglobin in the jugular vein, 42 indicating that the respiratory signals are indeed dominated by venous blood.
However, an arterial influence on the respiratory signals may also be possible.Once the oxygen saturation of hemoglobin of a volume-oscillating vascular compartment is determined, one is left with its interpretation on the basis of the arterial, capillary, and venous composition of the oscillatory blood volume.1][32][33][34][35] The question to tackle then pertains to the arterial and venous contributions to blood volume oscillations.
We have evaluated three different scenarios of possible blood volume oscillations, specifically T ¼ T ðaÞ , T ¼ T ðaÞ þ T ðvÞ , and T ¼ T ðvÞ .For T ¼ T ðaÞ þ T ðvÞ , the underlying assumption was that arterial volume changes and venous volume changes occur synchronously (so that the phasors T ðaÞ and T ðvÞ are in phase with each other).If, however, arterial and venous volume changes are not in phase with each other, Arg½T ðaÞ − Arg½T ðvÞ ≠ 0 deg, it follows that jT ðaÞ j þ jT ðvÞ j ≠ jT ðaÞ þ T ðvÞ j, and separating arterial from venous volume oscillations would not be possible in the absence of additional information.
In order to separate the arterial and venous volume contributions to the optical signals, the factor ρ ðvÞ (which represents the ratio of the venous-to-total oscillatory blood volume) must be known.A direct measurement of ρ ðvÞ may not be possible, even though specially designed physiological maneuvers may be used to estimate it.Of course, there may be cases where it can be reasonably assumed that the volume-oscillating vascular compartment is solely arterial (ρ ðvÞ ¼ 0) or venous (ρ ðvÞ ¼ 1).We stress again, however, that even in these special cases where the source of blood volume oscillations are known, one may still need to account for oscillatory blood flow contributions to the optical signals.This latter point was the main objective of this work.

Conclusion
We have presented an analytical method to translate dynamic NIRS measurements of cerebral hemoglobin concentrations [ΔOðtÞ and ΔDðtÞ] into the oxygen saturation of hemoglobin of the volume-oscillating vascular compartment, S V , without relying on the assumption that blood volume changes are the only source of the dynamics of the measured optical signals.This result requires taking into account the potential contributions of blood flow dynamics, which may indeed yield non-negligible contributions to the optical signals.These considerations are particularly significant in cerebral NIRS measurements.The Jana M. Kainerstorfer is an assistant professor of biomedical engineering at Carnegie Mellon University.Her research includes noninvasive optical imaging of tissue for disease detection and/or treatment monitoring, with a particular emphasis on diffuse optical spectroscopic imaging (with applications to the brain, breast, skin, and muscle).Her research spans from instrument development to novel data analysis tools all the way to clinical applications.Her research is published in more than 20 peer-reviewed scientific articles.
Angelo Sassaroli is research assistant professor of biomedical engineering at Tufts University.His research interests are within the fields of near-infrared spectroscopy and diffuse optical tomography for monitoring cerebral hemodynamics and breast cancer.He is mainly interested in developing analytical/computational methods for photon migration in tissue.He has coauthored more than 60 peer-reviewed scientific publications.
Sergio Fantini is a professor of biomedical engineering at Tufts University.His research interests are in diffuse optical spectroscopy and imaging of biological tissues, with applications to the human brain (functional imaging and assessment of cerebral hemodynamics), breast (cancer detection and monitoring response to therapy), and skeletal muscle (quantitative tissue oximetry).His research resulted in 11 patents and about 200 scientific publications, and he coauthored with Irving Bigio the textbook Quantitative Biomedical Optics.

Journal of
e m p : i n t r a l i n k -; e 0

Fig. 2
Fig.2Phasor diagram for the oxy-and deoxyhemoglobin phasors (O and D) at the cardiac pulsation frequency, with a negligible blood flow contribution.For this ideal case, O and D are in phase and the arterial saturation can be calculated as jOj∕jO þ Dj.

Fig. 3
Fig.3Optical helmet configuration.Light gray circles correspond to the illumination fibers, dark gray circles correspond to the detection fibers.The squares indicate the measurement channels, numbered from 1 to 24, which correspond to the tissue regions between each source-detector pair used for NIRS measurements.

Fig. 4 Fig. 5
Fig. 4 Representative dataset for subject no. 1 at the cardiac frequency.(a) Band-pass filtered time traces of ΔI∕I 0 at both wavelengths for channel 18.(b) The corresponding time traces of ΔO, ΔD, and ΔT with ArgðDÞ − ArgðOÞ ¼ −269 deg.Average phase differences between (c) the intensities at two wavelengths, and (d) the deoxy-and oxyhemoglobin concentrations.The circle in panel (c) encloses those channels where the phase differences are significantly different from zero.
Arg½D − Arg½O ranges from −342 to −215 deg [Fig.7(a) light gray], with a peak at −255 deg.Arg½O F − Arg½O shows a much smaller spread of values, ranging from −91 to −34 deg [Fig.7(b) light gray], with a peak at −72 deg.The dark gray histograms show the experimental data.In order to obtain the histogram of experimental Arg½O F − Arg½O, S ðaÞ was assumed to be 0.98.The model prediction for Arg½D − Arg½O and Arg½O F − Arg½O fall right into the range of the measured values.

Figure 8
Figure 8 shows a representative dataset for subject no. 1 during paced breathing at a frequency of 0.1 Hz.The filtered time traces of ΔI 700 ðtÞ∕I 0 700 and ΔI 830 ðtÞ∕I 0 830 , as well as ΔOðtÞ, ΔDðtÞ, and ΔTðtÞ at the respiration frequency of 0.1 Hz are shown.An average phase difference of Argði 700 Þ − Argði 830 Þ ¼ 99 AE 17 deg was found for channel 18.The phase difference Arg½D − Arg½O was found to be 200 AE 10 deg.Figures 8(c) and 8(d) show the phase difference of the intensity phasors and of D and O, respectively, for all channels of subject no.1.As described above (Sec.3.1) in the case of data filtered at the cardiac frequency, the model was used to predict Arg½D − Arg½O and Arg½O F − Arg½O for the data filtered at the respiration frequency of 0.1 Hz (as applied for 1 Hz as well).The model parameters were iterated as described above, but setting S ðaÞ ¼ 98%.While the autoregulation high-pass filter did not influence the phase of the hemoglobin concentration oscillations at 1 Hz (which is much higher than the filter cutoff frequency),

Fig. 6 Fig. 7
Fig.6Representative phasor diagram for subject no. 1, channel 18, at the cardiac frequency.Based on assuming S ðaÞ ¼ 0.98, O F and D F , as well as O V and D V could be calculated as shown.

Fig. 8
Fig. 8 Representative dataset for subject no. 1 during paced breathing at 0.1 Hz.(a) Band-pass filtered time traces of ΔI∕I 0 at both wavelengths for channel 18.(b) The corresponding time traces of ΔO, ΔD, and ΔT .Average phase differences between (c) the intensities at two wavelengths, and (d) the deoxy-and oxyhemoglobin concentrations.

Fig. 9
Fig. 9 Histograms of model predictions for Arg½D − Arg½O and Arg½O F − Arg½O at the paced breathing frequency of 0.1 Hz.Three different scenarios were simulated, with T ¼ T ðaÞ (red), T ¼ T ðaÞ þ T ðv Þ , with T ðaÞ ¼ T ðvÞ (green), and T ¼ T ðvÞ (blue).While Arg½D − Arg½O varies over a wide range of possible values, Arg½O F − Arg½O peaks in a small range of values between −9 and −4 deg.

Fig. 10
Fig. 10 Phasor diagram for subject no. 1, channel 18, at the paced breathing frequency of 0.1 Hz, obtained by assuming Arg½O F − Arg½O ¼ −7 deg, as estimated by the hemodynamic model.The resulting oxygen saturation of hemoglobin of the volume oscillating compartment is S V ¼ 0.67.
the exact value was unknown here.Varying f ðARÞ c from 0.01 to 0.05 Hz, we have found peak values of Arg½O F − Arg½O at −13 deg for f ðARÞ c ¼ 0.01 Hz and −1 deg for f ðARÞ c ¼ 0.05 Hz.
Fig.11Results of the analysis of the data during paced breathing for all subjects and measurement channels, demonstrating the difference between the ratios of jOj∕jTj, jOj∕ðjOj þ jDjÞ, and reporting the measured oxygen saturation of hemoglobin of the volumeoscillating vascular compartment.The mean values of the measured quantities were jOj∕jTj ¼ 1.1 AE 0.14, jOj∕ðjOj þ jDjÞ ¼ 0.88 AE 0.05, and hS V i ¼ 0.66 AE 0.14.