6 November 2017 Noninvasive in vivo optical characterization of blood flow and oxygen consumption in the superficial plexus of skin
Author Affiliations +
J. of Biomedical Optics, 22(11), 115002 (2017). doi:10.1117/1.JBO.22.11.115002
Assessing the metabolic activity of a tissue, whether normal, damaged, aged, or pathologic, is useful for diagnosis and evaluating the effects of drugs. This report describes a handheld optical fiber probe that contacts the skin, applies pressure to blanch the superficial vascular plexus of the skin, then releases the pressure to allow refill of the plexus. The optical probe uses white light spectroscopy to record the time dynamics of blanching and refilling. The magnitude and dynamics of changes in blood content and hemoglobin oxygen saturation yield an estimate of the oxygen consumption rate (OCR) in units of attomoles per cell per second. The average value of OCR on nine forearm sites on five subjects was 10 ± 5 (amol/cell/s). This low-cost, portable, rapid, noninvasive optical probe can characterize the OCR of a skin site to assess the metabolic activity of the epidermis or a superficial lesion.
Liasi, Samatham, and Jacques: Noninvasive in vivo optical characterization of blood flow and oxygen consumption in the superficial plexus of skin



Real-time monitoring of the cellular oxygen consumption rate [OCR, attomoles of O2 per cell per s (amol/cell/s)] can characterize the metabolic activity of a tissue, providing physiological information with clinical relevance. For instance, the metabolic rate may provide a basis for comparing normal and pathological states. Prior studies have used other methods such as respirometry of extracted mitochondria ex vivo and in vivo multiphoton microscopy to assess the OCR of tissues.12.3 This paper focuses on assessing the OCR of the skin’s epidermis. Monitoring the epidermal OCR may characterize normal, diseased, and aged skin, detect malignancy, and monitor the effects of drugs and treatments.

Different values of OCR for various cells and tissues have been reported in the literature. Based on data from in vitro studies, Wagner et al.4 cataloged a wide range of OCR from 1 to 1000  amolcell1s1, with a median value of 23. Balu et al.1 used multiphoton microscopy in vivo to image NADH fluorescence in the epidermis as a proxy for OCR. They reported an OCR of about 0.035  μmoles/106  cells/h in keratinocytes in the basal layer, which equals 9.7  amolcell1s1. They also reported a nonuniformity of OCR across the epidermal thickness, with more superficial keratinocytes having lower OCR.

This paper presents an alternate method using a topical optical fiber probe measurement at the skin surface to assess the blood volume and hemoglobin oxygen saturation before, during, and after manually applied pressure to blanch the superficial vascular plexus that provides oxygen to the epidermis. The refill time after release of pressure specifies the blood flow in the blanched superficial plexus. The magnitude and dynamics of changes in blood content and hemoglobin oxygen saturation yield an estimate of the epidermal OCR in units of amol cell1s1. The goal of this paper is to report a pilot study that illustrates the method and cites initial values for the in vivo epidermal OCR.

Previous work has used applied pressure to remove blood from a tissue for the purpose of improving the observation of other chromophores and tissue structures. For example, Ermakov and Gellermann detected carotenoids in the skin after using pressure to remove superficial blood.5 Such reports use pressure to remove the blood, which is regarded as an interfering factor. In this paper, we use pressure to identify the superficial vascular plexus, which is easily blanched by pressure. The spectroscopy of this blanchable plexus informs about the OCR of the epidermis.




Experimental Design

An optical fiber probe was used to spectrally measure the transport of light through in vivo forearm skin sites (Fig. 1). The probe consisted of a source fiber and a collection fiber with a center-to-center spacing of 1.88 mm (0.400-mm-core-dia fibers). The fibers were embedded in a flat plastic plate (three-dimensional-printed using polylactate) that allowed perpendicular fiber contact with the skin while distributing the pressure applied to the skin. The probe handpiece was a cylinder within which the optical probe assembly moved. As one manually pressed the probe against the skin, the outer cylinder pushed a spring against the inner probe assembly, hence applying pressure by the spring. The displacement of the spring was observable on a ruler attached to the inner assembly, which extended outside of the top of the cylinder to allow easy viewing. A reproducible pressure was applied by pressing the probe until the desired displacement was achieved. Repeated tests of the blanching force registered 610±10  g on a scale. The flat surface of the probe tip distributed the force, hence, the probe could also be placed lightly on the skin to monitor refill while avoiding blanching the superficial plexus of blood.

Fig. 1

The optical fiber probe placed on the forearm. Two optical fibers deliver and collect white light for measurement by a spectrometer. When the handpiece is pressed against the skin, pressure is applied by a spring whose displacement is monitored by reading the ruler. The inset shows the probe tip with two 0.400-mm-core-dia fiber separated center-to-center by 1.88 mm embedded in a broad circular plastic surface. The applied pressure is distributed over this circular surface, so the pressure is more easily controlled.


The source fiber delivered white light. The collection fiber returned light to a spectrometer (QE Pro, Ocean Optics Inc., Dunedin, Florida) which measured the light in the 300 to 1100 nm range. Light in the 500 to 1000 nm range yielded strong signals and was used for the analysis. The spectra were acquired continuously by a digital chart recorder (implemented using a MATLABTM program with a GUI interface). Spectra with light pressure were acquired for about 30 s to establish a baseline. Then, strong pressure was applied for about 10 s, which blanched the superficial blood volume, i.e., the superficial vascular plexus. Then light pressure was resumed, which allowed the superficial plexus to refill with blood.

Spectral analysis used the 500 to 1000 nm range of wavelengths. Light below 500 nm (blue light) is scattered by the epidermis and superficial papillary dermis before encountering significant blood content; therefore, spectra at such short wavelengths require a two-layer model for analysis. Light above 500 nm (green to near-infrared) penetrates into the dermis and samples the superficial plexus of microvasculature at about 200 to 400  μm depth. Hence, a simple diffusion theory that assumes uniform distribution of blood could be used (see Sec. 2.3).

The analysis of the spectrum at each time point specified the apparent average blood volume fraction (B) and the mixed arteriovenous oxygen saturation (S) (see Sec. 2.3). Hence, the time courses of B(t) and S(t) were specified. The values of B and S after applying pressure for about 10 s (t1) and after the refilling of the plexus (t2) were noted: t1 = at maximum blanching, t2 = after refill of the plexus upon release of pressure. The data at these time points were used to calculate the OCR.

The experiment was conducted on ventral forearm skin sites: nine measurements on five subjects. For each subject, the probe was placed on the forearm, and pressure was applied then released. Then, the probe was moved to a different site on the forearm and the experiment repeated, eventually measuring three to six sites. Measurements over a visible large blood vessel were avoided. Measurements were approved by the Institutional Internal Review Board at OHSU.


Calculation of OCR

Calculation of the OCR (amol/cell/s) was accomplished by the following:

  • 1. Blood flow


    J=(ΔB/Δt)(Bplexus/Bsp)Vplexus(cc/spercm2surface area)
    where ΔB/Δt=[B(t2)B(t1)]/(t2t1), change in apparent blood per unit time (s1); Bplexus/Bsp=2.5= ratio of blood in plexus versus blood observed spectrally (see below); Dplexus=0.0200  (cm) thickness of superficial vascular plexus (estimate); Vplexus=Dplexus(1  cm2)=0.0200 (ccpercm2 surface area).

  • 2. Oxygen extraction


    O2.ex=4ΔSJCHGb(molesO2/spercm2surface area),
    where 4=fourO2 molecules per hemoglobin molecule; f=Vv/Va, the ratio of venous volume to arterial volume (assumed); fa=1/(1+f), volume fraction of blood that is arterial; fv=f/(1+f), volume fraction of blood that is venous; Sa=0.97, arterial oxygen saturation (assumed); Sv=(SfaSa)/fv, venous oxygen saturation; ΔS=SaSv, difference in oxygen saturation between arterial and venous; CHGb=Cblood/MW, concentration (moles/cc); Cblood=0.150  g.HGb/cc, hemoglobin concentration of whole blood; MW=64,500  g/mole, molecular weight of hemoglobin.

  • 3. Number of cells in epidermis per cm2 surface area


    where Vepidermis=(0.0100  cm)(1  cm2) (ccpercm2 surface area); Vcell=πD3/6(cc/cell); D=10×104  cm diameter of spherical epidermal cell (approximation).

  • 4. OCR



The factor Bplexus/Bsp=2.5 is the ratio of blood in the plexus versus the apparent average blood deduced by diffusion theory from the observed spectrum. The spectral analysis uses a simple diffusion theory [see Eq. (7)] that assumes a homogeneous skin tissue with an overlying pigmented epidermis. Any observed blood absorption is assumed to be due to blood distributed evenly throughout a semi-infinite homogeneous skin. But the superficial vascular plexus has a higher local blood content (Bplexus) than the average (Bsp) and is located at about 200 to 400  μm depth (a 200-μm thickness assumed here). The spectrally observed average Bsp does not equal the true local Bplexus.

A Monte Carlo simulation (MCML6) was used to simulate the reflectance spectra (15 wavelengths from 500 to 1000 nm) when such a plexus was given different blood volume fractions (Bplexus=0.001 to 0.05). The spectra were then analyzed to specify an apparent Bsp, using the same diffusion theory analysis that was applied to the experimental data. There was a linear relationship between Bplexus and Bsp with a slope of 2.5. Hence, when spectral analysis of a measured spectrum specified Bsp to be 0.01, the true Bplexus was 0.025.

The relative sizes of the arterial and venous vascular volumes are an important parameter in the calculation of OCR. The ratio of volumes, f=Vvenous/Varterial, was assumed to be 1.0 in this paper. Equation (2), however, allows for f to deviate from 1.0.


Spectral Analysis

The measurements yielded a reflectance spectrum, Mskin(λ), a spectrum from a white standard (99% reflectance standard, SpectralonTM, Labsphere Inc., New Hampshire) with the probe held 1 cm above the standard, M99%(λ), and a spectrum with the probe held in air within a dark room, M0%(λ). The dark spectrum was subtracted from the other spectra, then a ratio was taken




where K=Cskin/C99%. This ratio cancels the wavelength dependence of the light source (S) and detector responsivity (D). The left side of Eq. (6) is known from measurements and R99%=0.99. The right side of Eq. (6) is fit by least squares fitting using the modified diffusion theory described in Eq. (7). After the fitting, the skin reflectance Rskin=MR/K is in units of (mm2). The factor K depends on the ratio of collection efficiency by the probe when held about 1 cm above the standard versus collection efficiency when the probe contacted the skin.

Figure 2(a) shows the raw reflectance spectra in units of (counts): M99%, M0%, and Mskin± applied pressure. Figure 2(b) shows the reflectance Rskin KR(λ)± applied pressure. The factor K equals 2.5 based on the average value fit for all acquired spectra on the five subjects.

Fig. 2

(a) The raw spectra Mskin±; pressure, M99%, and M0% (counts). The M99% is scaled to fit on graph. (b) The reflectance spectrum, where source and collection fiber contact the skin.


The spectrum KRskin was matched by least squares fitting using a simple light transport expression based on diffusion theory [Eq. (7)]. A predicted spectrum pKR was determined by specifying the fitting parameters: blood content (B), oxygen saturation (S), epidermal melanin (M), scattering (a), and a calibration factor (K). These parameters are defined below. The parameter a scales the magnitude of the reduced scattering coefficient, μs(λ), whose spectral shape is described below. The expression for pKR (counts) was


where K=calibration constant (counts mm2); Tepi=exp(Mμa.melLepi), the attenuation due to the pigmented epidermis; M=volume fraction of melanosomes in 60-μm-thick forearm epidermis; μa.mel=absorption coefficient of interior of cutaneous melanosome;7 Lepi3× epidermal.thickness =0.180  mm; μa is the absorption coefficient (mm1)
where μa.oxy=absorption spectrum of oxygenated whole blood;8 μa.deoxy=absorption spectrum of deoxygenated whole blood;8 μs=reduced scattering coefficient (mm1) of skin, based on literature;9 a=μs at λ=500  nm, which scales the shape of μs(λ) such that; μs(λ)=a[0.32(λ/500  nm)4+0.68(λ/500  nm)0.91] = Rayleigh plus Mie scattering;9 getRr(μa,μs,r,ri) = local reflectance (mm2) (diffusion theory, Farrell et al.10); r=1.88  mm center-to-center distance between source and collection fibers. nr=1.4/1.5, ratio of the approximate refractive indices, nskin/nprobe.

Hence, the spectral fitting yielded values for B,S,M,a, and K. After an initial analysis of all spectra on the five subjects, a common value of K was determined to be 2.5. The same M99% was used to analyze all the spectra. Subsequent analysis fit for B,S,M, and a only. The water content was assumed to be W=0.65 for all subjects. Figure 3 shows the spectral fitting of the nine spectra from five subjects.

Fig. 3

Measurements of reflectance spectra, KRskin (cm2), from the forearms of five subjects (red), and the diffusion theory fit to the spectra (black). Subject 2 was the most strongly pigmented subject. (calibration K=2.5).




First, the calculation of OCR is illustrated by applying the calculation to one of the pressure-release data from one subject in Fig. 4. The calculation steps are listed in Table 1. The mean ±standard deviation values of OCR for the six pressure releases on that subject yielded an OCR=9.1±2.7  amol/cell/s, which is close to the value of 9.7  amol/cell/s reported by Balu et al.1

Fig. 4

The time course of signals. (a) MR(t), showing how applied pressure blanches the skin, causing an increase in reflectance. (b) B(t), showing how pressure reduces the blood content, i.e., blanches the skin, and release of pressure yields an initial linear rise in B, specifying the slope dB/dt. (c) S(t), showing the modulation of the oxygen saturation of hemoglobin in the mixed arteriovenous vasculature. In all three figures, the arrows point to bold black lines that show the same period t1 to t2, during which there is linear refilling of the plexus.


Table 1

Example calculation from the fourth pressure release in Fig. 3.

B(t1)=B measured while blanched0.00121
B(t2)=B measured after refill0.00313
S(t2)=S measured after refill0.655
dB/dt=[B(t2)B(t1)]/(t2t1)5.99×104  s1
Bratio=Bplexus/Bsp from Monte Carlo sims.2.5
Vplexus0.02  cc/percm2 skin surface
J=(dB/dt)BratioVplexus2.99×105  cc/s
Cblood = conc. of HGb in whole blood0.150  moles/cc
MWHGb64,500  moles/cc
CHGb=Cblood/MW2.33×106  moles/cc
Sa assumed0.97
Dcell0.0010 cm
OCR=(ΔO2/ncells)10189.2  amol/cell/s

Second, the analysis of the nine skin sites on five subjects is summarized in Table 2. The means and standard deviations of the OCR on each of the nine skin sites are listed. The average coefficient of variation for measurements on a site was 34%, which is site-to-site variation on the forearm during one sitting of a subject. The range of mean OCR values extended from 0.9 to 18.0  amol/cell/s. The average mean OCR value was 10±5  amol/cell/s (n=9 skin sites). There was significant sitting-to-sitting variation on subjects 1, 2, and 5, who were measured more than once.

Table 2

Summary of OCR for nine sites on five subjects.

SiteSubjectMean±std (n pressure releases)
119.1±2.7 (n=6).
2112.9±2.3 (n=4)
3115.8±5.0 (n=3)
420.9±0.2 (n=3)
529.9±2.5 (n=4)
6311.4±5.7 (n=5)
744.1±2.0 (n=3)
859.1±5.8 (n=3)
9518.0±4.3 (n=3)

Figure 5 places the OCR data of this report on a histogram drawn using the tabulated data in the review of OCR values for various tissues and cell cultures by Wagner et al.4 A couple of the skin site data are on the edge of the histogram but most are clustered in the middle. Hence, the skin OCR values are in the more common range found by various investigators.

Fig. 5

The distribution of OCR (amol/cell/s) for various cells and tissues, as compiled in the review by Wagner et al. 2011. The red lines indicate the OCR of the nine skin sites measured in this report.




Topical placement of an optical fiber probe on the skin followed the dynamics of blanching by applied pressure then refilling upon release of the pressure. A spectrogram of reflectance spectra versus time was acquired and analyzed by a simple diffusion theory model to yield the time course of blood volume fraction (B) and mixed arteriovenous oxygen saturation of hemoglobin (S). The B(t) and S(t) recordings specified the blanched value at t1 and the refilled values at t2, which yielded values for the blanchable blood volume and the blood flow to this volume, which is presumed to be the superficial plexus just under the epidermis. This plexus exclusively provides oxygen for the epidermis and does not supply significant oxygen to the dermis. So, the oxygen delivery of the plexus is an indicator of the metabolic consumption of oxygen by the epidermis.

Balu et al.1 used two-photon NADH fluorescence to quantify the gradient of OCR across the epidermis. More work is needed to understand the difference between the average OCR measured by the optical fiber measurement and the OCR gradient that they observed. One of the assumptions in the fiber probe analysis is that the rate of cellular oxygen extraction from the blood is equal to the cellular OCR, which assumes that all of the oxygen used by the epidermis comes from the blood. On the other hand, some studies have reported a contribution of atmospheric oxygen to the epidermal oxygen supply.1112.13

Gaining a better understanding of epidermal oxygen consumption may prove useful in differentiating normal from pathological states. The optical fiber probe is a low-cost portable system that offers a convenient, rapid, noninvasive assessment of the metabolic activity of the epidermis, or any pathologic growth near the skin surface. Such a probe may prove useful in clinical dermatology, which was the motivation for this study.


The authors have no conflicts of interests to disclose.


The authors would like to thank Dr. Sancy Leachman for arranging Dr. Liasi’s participation in the research project at OHSU.


1. M. Balu et al., “In vivo multiphoton NADH fluorescence reveals depth-dependent keratinocyte metabolism in human skin,” Biophys. J. 104, 258–267 (2013). http://dx.doi.org/10.1016/j.bpj.2012.11.3809 Google Scholar

2. M. F. Forni et al., “Bioenergetic profiling in the skin,” Exp. Dermatol. 25(2) 147–148 (2016).EXDEEY0906-6705 http://dx.doi.org/10.1111/exd.12856 Google Scholar

3. D. Pouli et al., “Imaging mitochondrial dynamics in human skin reveals depth-dependent hypoxia and malignant potential for diagnosis,” Sci. Transl. Med. 8, 367ra169 (2016).STMCBQ1946-6234 http://dx.doi.org/10.1126/scitranslmed.aag2202 Google Scholar

4. B. A. Wagner, S. Venkataraman and G. R. Buettner, “The rate of oxygen utilization by cells,” Free Radical Biol. Med. 51(3), 700–712 (2011).FRBMEH0891-5849 http://dx.doi.org/10.1016/j.freeradbiomed.2011.05.024 Google Scholar

5. I. V. Ermakov and W. Gellermann, “Dermal carotenoid measurements via pressure mediated reflection spectroscopy,” J. Biophotonics 5(7), 559–570 (2012). http://dx.doi.org/10.1002/jbio.201100122 Google Scholar

6. L. Wang, S. L. Jacques and L. Zheng, “MCML—Monte Carlo modeling of photon transport in multi-layered tissues,” Comput. Meth. Programs Biomed. 47(2), 131–146 (1995) http://dx.doi.org/10.1016/0169-2607(95)01640-F Google Scholar

7. S. L. Jacques and D. J. McAuliffe, “The melanosome: threshold temperature for explosive vaporization and internal absorption coefficient during pulsed laser irradiation,” Photochem. Photobiol. 53(6), 769–775 (1991).PHCBAP0031-8655 http://dx.doi.org/10.1111/php.1991.53.issue-6 Google Scholar

8. S. A. Prahl, “Hemoglobin,”  http://omlc.ogi.edu/spectra/hemoglobin/index.html (28 May 2017). Google Scholar

9. S. L. Jacques, “Optical properties of biological tissues: a review,” Phys. Med. Biol. 58(11), R37–R61 (2013).PHMBA70031-9155 http://dx.doi.org/10.1088/0031-9155/58/11/R37 Google Scholar

10. T. J. Farrell, M. S. Patterson and B. Wilson, “A diffusion theory model of spatially resolved, steady-state diffuse reflectance for the noninvasive determination of tissue optical properties in vivo,” Med. Phys. 19(4), 879–888 (1992).MPHYA60094-2405 http://dx.doi.org/10.1118/1.596777 Google Scholar

11. G. Ronquist et al., “Human epidermal energy metabolism is functionally anaerobic,” Exp. Dermatol. 12(5), 572–579 (2003).EXDEEY0906-6705 http://dx.doi.org/10.1034/j.1600-0625.2003.00018.x Google Scholar

12. M. Stuecker et al., “The transepidermal oxygen flux from the environment is in balance with the capillary oxygen supply,” J. Invest. Dermatol. 114(3), 533–540 (2000).JIDEAE0022-202X http://dx.doi.org/10.1046/j.1523-1747.2000.00911.x Google Scholar

13. M. Stuecker et al., “The cutaneous uptake of atmospheric oxygen contributes significantly to the oxygen supply of human dermis and epidermis,” J. Physiol. 538(3), 985–994 (2002).JPHYA70022-3751 http://dx.doi.org/10.1113/jphysiol.2001.013067 Google Scholar


Faezeh Talebi Liasi is a medical doctor and an alumni of the University of Washington School of Medicine. She worked on this project at the Oregon Health and Science University (OHSU) while she was still in medical school. Her research interests include optical characterization of tissue and its use in diagnostics.

Ravikant Samatham received his MS degree in mechanical engineering from the University of Nevada, Reno, in 2004 and a PhD in biomedical engineering from OHSU in 2012. He is now a senior research scientist in the Department of Dermatology at OHSU. His research focuses are second harmonic generation microscopy, optical spectroscopy, and optical imaging.

Steven Jacques received his BS degree in biology from MIT in 1972, his MSEE degree in electrical engineering and computer science, and his PhD in biophysics and medical physics from the University of California-Berkeley in 1984. He has served at Massachusetts General Hospital/Harvard Medical School, the University of Texas M.D. Anderson Cancer Center, the OHSU, and now at Tufts University. His research focuses on light propagation and interactions with tissues.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Faezeh Talebi Liasi, Ravikant Samatham, Steven L. Jacques, "Noninvasive in vivo optical characterization of blood flow and oxygen consumption in the superficial plexus of skin," Journal of Biomedical Optics 22(11), 115002 (6 November 2017). http://dx.doi.org/10.1117/1.JBO.22.11.115002 Submission: Received 23 June 2017; Accepted 16 October 2017
Submission: Received 23 June 2017; Accepted 16 October 2017


Optical character recognition



Tissue optics

In vivo imaging

Blood circulation


Back to Top