# Skin Temperature Influence on Transcutaneous Carbon Dioxide (CO 2) Conductivity and Skin Blood Flow in Healthy Human Subjects at the Arm and Wrist

Emmanuel Dervieux, François Guerrero, Wilfried Uhring, Marie-Agnès Giroux-Metgès, Michael Theron

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Emmanuel Dervieux, François Guerrero, Wilfried Uhring, Marie-Agnès Giroux-Metgès, Michael Theron. Skin Temperature Influence on Transcutaneous Carbon Dioxide (CO 2) Conductivity and Skin Blood Flow in Healthy Human Subjects at the Arm and Wrist. Frontiers in Physiology, 2023, 14, pp.1293752. 10.3389/fphys.2023.1293752. hal-04301290

HAL Id: hal-04301290

https://hal.science/hal-04301290v1

Submitted on 23 Jan 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN ACCESS

EDITED BY
Xudong Ge,
University of Maryland, Baltimore
County, United States

REVIEWED BY
Venkatesh Srinivasan,
University of Maryland, Baltimore
County, United States
Kenta Litani,
Tokyo Medical and Dental University,
Japan

*CORRESPONDENCE
Emmanuel Dervieux,
emmanuel.dervieux@biosency.com

RECEIVED 13 September 2023

ACCEPTED 05 December 2023

PUBLISHED 23 January 2024

**CITATION**
Dervieux E, Guerrero F, Uhring W, Giroux-Metgès M-A and Théron M (2024). Skin temperature influence on transcutaneous carbon dioxide (CO₂) conductivity and skin blood flow in healthy human subjects at the arm and wrist.
*Front. Physiol.* **14**:1293752.
doi: 10.3389/fphys.2023.1293752

**COPYRIGHT**
© 2024 Dervieux, Guerrero, Uhring, Giroux-Metgès and Théron. This is an open-access article distributed under the terms of the [Creative Commons Attribution License (CC BY)](https://creativecommons.org/licenses/by/4.0/). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

# Skin temperature influence on transcutaneous carbon dioxide (CO₂) conductivity and skin blood flow in healthy human subjects at the arm and wrist

Emmanuel Dervieux¹,²,³*, François Guerrero², Wilfried Uhring³, Marie-Agnès Giroux-Metgès²,⁴ and Michaël Théron²

¹Biosency, Cesson-Sévigné, France, ²EA4324-ORPHY, Univ Brest, Brest, France, ³ICube, University of Strasbourg and CNRS, Strasbourg, France, ⁴Explorations Fonctionnelles Respiratoires, Centre Hospitalier Régional et Universitaire de Brest, Brest, France

**Objective:** present transcutaneous carbon dioxide (CO₂)—tcpCO₂—monitors suffer from limitations which hamper their widespread use, and call for a new tcpCO₂ measurement technique. However, the progress in this area is hindered by the lack of knowledge in transcutaneous CO₂ diffusion. To address this knowledge gap, this study focuses on investigating the influence of skin temperature on two key skin properties: CO₂ permeability and skin blood flow.

**Methods:** a monocentric prospective exploratory study including 40 healthy adults was undertaken. Each subject experienced a 90 min visit split into five 18 min sessions at different skin temperatures—Non-Heated (NH), 35, 38, 41, and 44°C. At each temperature, custom sensors measured transcutaneous CO₂ conductivity and exhalation rate at the arm and wrist, while Laser Doppler Flowmetry (LDF) assessed skin blood flow at the arm.

**Results:** the three studied metrics sharply increased with rising skin temperature. Mean values increased from the NH situation up to 44°C from 4.03 up to 8.88 and from 2.94 up to 8.11 m·s⁻¹ for skin conductivity, and from 80.4 up to 177.5 and from 58.7 up to 162.3 cm³·m⁻²·h⁻¹ for exhalation rate at the arm and wrist, respectively. Likewise, skin blood flow increased elevenfold for the same temperature increase. Of note, all metrics already augmented significantly in the 35–38°C skin temperature range, which may be reached without active heating—i.e. only using a warm clothing.

**Conclusion:** these results are extremely encouraging for the development of next-generation tcpCO₂ sensors. Indeed, the moderate increase (x 2) in skin conductivity from NH to 44°C tends to indicate that heating the skin is not critical from a response time point of view, i.e. little to no skin heating would only result in a doubled sensor response time in the worst case, compared to a maximal heating at 44°C. Crucially, a skin temperature within the 35–38°C range already sharply increases the skin blood flow, suggesting that tcpCO₂ correlates well with the arterial paCO₂ even at such low skin temperatures. These two conclusions
further strengthen the viability of non-heated tcpCO₂ sensors, thereby paving
the way for the development of wearable transcutaneous capnometers.

KEYWORDS

carbon dioxide, ptCO₂, tcpCO₂, transcutaneous, exhalation rate, diffusion, blood flow, hyperaemia

# 1 Introduction

Due to its clinical significance, continuous monitoring of the
arterial CO₂ partial pressure—paCO₂—is of paramount importance
in medical practice, especially for patients presenting severe
respiratory disorders (Wagner, 2015). The gold standard to get a
single paCO₂ reading consists in an arterial puncture followed by
a gaseous analysis of the collected blood sample. Unfortunately,
this procedure—tersely referred to as the “blood gases” in
clinical settings—is both painful and risky (Scheer et al., 2002),
requiring trained personnel as well as expensive blood gas
analysers. It also calls for a quick analysis of the blood samples
following their collection, which adds stress to hospital logistics
(Nanji and Whitlow, 1984). These major drawbacks led to the
development of transcutaneous CO₂ monitors, which consist
in a Stow-Severinghaus electrode—mainly a pH-meter bathing
in a bicarbonate solution—heated in the 41-44°C range, and
placed against a patient's skin (Severinghaus and Astrup, 1986;
Huttmann et al., 2014). This electrode measures a transcutaneous
CO₂ partial pressure—the tcpCO₂—which correlates well with the
paCO₂ if the skin is heated above at least 38°C (Wimberley et al.,
1985).

Yet, these monitors also suffer from several weaknesses: i) their important drift requires a recalibration with an appropriate gas mixture every 8 hours, at most (Bendjelid et al., 2005), ii) the thin membrane covering the electrode is fragile, and needs to be replaced every 2 weeks or so (Lermuzeaux et al., 2016), iii) the heating power required by the electrode to maintain the skin in the above-mentioned temperature range is about 100-200 mW<fnref n="1" />, thus precluding its use in a battery-powered wearable stand-alone device and iv) their elevated price tag prevents their widespread use, would it be in a clinical or home-based setting. Consequently, the development of an alternative to the existing tcpCO₂ monitors appears mandatory, and has been an active research field in the last decades (Dervieux et al., 2022; Section 4.1).

Recently, in a review article aiming at encompassing the diversity
of CO₂ measurement techniques with a focus on biomedical
applications, we divided the issue of developing such an alternative
tcpCO₂ monitor into three research areas (Dervieux et al., 2022):

1. Due to the above-mentioned drift, and high cost of the Stow-Severinghaus electrode, an alternative CO₂ *measurement technique* is needed.

2. Then, in order to dimension the sensor-to-be, it is essential to accurately know the CO₂ *exhalation rate* through the skin, as the latter directly influences the response time of the sensor.

3. Finally, it is mandatory that the tcpCO₂ and paCO₂ are in good agreement at the skin temperature considered for measurement. *Id est*, that the *correlation* between tcpCO₂ and paCO₂ is satisfactory at the latter temperature.

Regarding the first point, it appeared to us that, among the many
technologies reviewed, a polymer patch embedding a CO₂-sensitive
fluorophore would be particularly advisable. Interestingly, this trail
has recently been followed by Cascales et al. (2022) or Tufan and
Guler (2022) with some success, although no in vivo experiment
have been conducted to date. This point will be the object of future
studies and is not developed any further in this paper.

The second and third points, on the contrary, are at the very heart
of the present study. Starting with the exhalation rate, the main issue
with data available in the literature—see Table 3—is that the skin
temperature is only mentioned once—by Eöry (1984)—and never
accurately regulated when this parameter is measured (Fitzgerald,
1957). Even though Levshankov et al. (1983) crafted a heating
device, they do not report the temperature setpoint that they used.
Thus, the present study aims at filling this gap by measuring the
influence of skin temperature on the transcutaneous CO₂ exhalation
rate.

While the tcpCO₂/paCO₂ correlation is excellent at—or
above—42°C (Conway et al., 2018), scarce are the authors
who investigated lower skin temperatures, with none but
Wimberley et al. (1985) experimenting with temperatures as low
as 38°C. The reason for heating the skin in the first place is to trigger
a local reactive hyperaemia (Roustit and Cracowski, 2012). By doing
so, the subcutaneous tissues are flushed with fresh arterial blood, and
their gaseous content thus gets closer to the arterial one (Koch, 1965;
Rooth et al., 1987; Zavorsky et al., 2007). While temperatures in the
42-44°C range have often been used to trigger maximal hyperaemia,
lower temperatures have been seldom explored (Hodges et al.,
2016), and we thus took advantage of our exhalation rate
measurements to measure the skin blood flow at lower temperatures
simultaneously.

Then, the reader should bear in mind the importance of
skin temperature for designing a new kind of tcpCO₂ sensor.
Indeed—ideally—such a sensor should heat the skin as little as
possible for two main reasons: i) heating the skin is uncomfortable
for the patient and can wake them up in the case of night time
monitoring, and ii) it consumes a significant amount of power,
which precludes using such sensor in a wearable, as mentioned
above. Yet, such an unnoticeable and wearable tcpCO₂ sensor
would be highly desirable in a telemonitoring context for home
use, reducing the need for hospital visits. Indeed, if the positive
impact of tcpCO₂ telemonitoring is yet to be demonstrated—for the
obvious reason that the corresponding wearable tcpCO₂ monitor
does not exist at the time being—several clinical trials demonstrate
the beneficial contribution of telemedicine—a.k.a. telehealth—on

<footnote n="1">1 Typical value for a TCM4 tcpCO₂ monitor (Radiometer, Denmark), as measured by the authors.</footnote>
<figure><img src="image_4.png" /><figcaption>**FIGURE 1**
General outline of the rate sensor and its peripherals. See the text for further explanations.</figcaption></figure>

both patient's outcome and costs of admission in a variety of conditions (Steventon et al., 2012; Yun et al., 2018; Kruse et al., 2019). Additionally, the outbreak of contagious pandemics—such as COVID-19 (Garfan et al., 2021)—and the rapid development of the health wearable market (Dunn et al., 2018; Yetisen et al., 2018; Chung et al., 2019; Dagher et al., 2020) may also promote the use of telemonitoring in medical practice. For these reasons, not heating the skin while measuring tcpCO₂ would be highly desirable.

The following study focuses on measuring the transcutaneous CO₂ exhalation rate and cutaneous micro-circulation on the full NH-44°C skin temperature range. Two measurement sites were investigated: the dorsal side of the wrist and the lateral aspect of the upper arm, while the skin blood flow was only measured at the upper arm. Additionally, a strong emphasis was placed on the transcutaneous CO₂ conductivity, which may be preferred to the well-known exhalation rate because of its more intrinsic nature—in particular, the latter conductivity does not depend on the ambient CO₂ level, nor on the subject's paCO₂, as opposed to the exhalation rate, which is influenced by both.

# 2 Materials

## 2.1 The transcutaneous CO₂ rate sensor

A custom transcutaneous CO₂ diffusion rate sensor—hereafter simply denoted as “the sensor”—was developed for the needs of this study. Its basic working principle is close to that evoked by Dervieux et al. (2022)—and will be further detailed in Section 3.1.1—while its design is inspired by the early works of Eletr et al. (1978) and Greenspan et al. (1981). The general outline of the sensor and its peripherals can be seen in Figure 1.

The sensor, designed to be placed against the subject's skin by mean of a double-sided adhesive, is connected to three main apparatuses: i) a calibrated, reference thermometer (Testo 735, Testo, Germany) equipped with a type K thermocouple (110-4482, RS Pro, United Kingdom), ii) a Doppler perfusion monitor (Periflux 5000, Perimed, Sweden) equipped with a 407 probe, and iii) a control and supply block, consisting in a thermostat, a power supply unit, and

a Universal Serial Bus (USB) to universal Asynchronous Receiver Transmitter (UART) converter, embedded in a 3D-printed case. For the sake of conciseness though, the control and supply block is only detailed in the Supplementary Material S1, which also contains a thorough analysis of the safety issues that may arise when using this sensor. The sensor itself can be seen in great details in Figure 2. It consists in an aluminium (2017A) body, which serves as a support for the following elements: a CO₂ sensor, a heating resistive wire, a thermistor, a thermocouple, the Doppler probe, a poly-lactic acid (PLA) 3D-printed cover, and an interfacing Printed Circuit Board (PCB). Complete drawings of the sensor's body are provided in the Supplementary Material S1.

The CO₂ sensor is a MinIR (ExplorIR-M5%, CO₂ Meter, United States), an off-the-shelf, compact, Non Dispersive Infra-Red (NDIR) CO₂ sensor, with a full range of 5% and an accuracy of 70 ppm ± 5% of reading at Standard Temperature and Pressure (STP)—see Hodgkinson and Tatam (2012) for further details on the operating principle of such sensors. Its internals—a pair of IR Light Emitting Diode (LED) and photodiode (PD) facing a spherical, gold-plated mirror—may be seen in the cut view of Figure 2. The pCO₂ inside the sensor was recorded with a sampling frequency of 2 Hz.

As the gas-tightness of the measuring chamber is a critical aspect of the sensor's operating principle—see Section 5.1.3—a two-stage sealing was implemented: i) a silicone-grease coated (515520, GEB, France), soft—60 Shore A hardness—silicone O-ring was placed between the aluminium body and the CO₂ sensor itself and ii) liquid epoxy resin (Résine Cristal, Gédéo, France) was cast in the remaining interstice between the latter two elements—illustrated in the cut view of Figure 2 in vivid red.

The thermoregulation of the sensor is performed by mean of a resistive wire for heating, coupled to a thermistor for temperature measurement and regulation. For verification purpose, an additional thermocouple was also added, as stated above. The heating wire consists in 2 × 15 turns of 28 Ω·m⁻¹, 0.15 mm in diameter, enamelled, resistive, constantan wire (Isotan, Thomsen), connected in parallel. The wire delivers a total heating power of 6.1 W under 12 V, and is coiled around the aluminium body, in a dedicated groove. The bottom of the groove is covered with a layer of 0.25 mm thermally conductive double sided tape (8810, 3M, United
<figure><img src="image_5.png" /><figcaption>FIGURE 2</figcaption></figure>

Detailed views of the sensor. (A): exploded view, detailing its different parts. (B): cut view, showing the inner functioning principle of the CO₂ sensor.

Note the epoxy resin sealing, in red. (C): isometric views from above, and below, of the fully assembled sensor, illustrating the grid-shaped sensor's sole.

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

States) prior to coiling the wire, and the latter is finally covered with two nitrile quad rings, as can be seen in **Figure 2**. This covered layout prevents burns caused by direct contact with the heating wires. The thermistor (151-237, RS Pro, United Kingdom) and thermocouple were glued in two dedicated flat-bottom mounting holes which were pre-filled with a thermally conductive, electrically non-conductive, epoxy resin (8329TFM, MG Chemicals, Canada). Care was taken that *i*) the distance between the bottom of the mounting holes and the heating wire and *ii*) that between the sole of the sensor's body and the heating wire were equal, in order to ensure that the temperature measured by the thermistor and thermocouple is close to that of the skin.

The Doppler probe is housed in a dedicated hole, and can slide vertically, in such a way that it can be adjusted to outcrop the sole of the sensor, coming in direct contact with the skin. It holds in place by mean of a cup-pointed, headless, set screw which compresses it radially via an O-ring, so as not to damage the probe. The raw Doppler perfusion signal—originally sampled at 62.5 Hz—was downsampled to 0.625 Hz and low-pass filtered using a tenth-order Butterworth filter prior to further analysis—see **Section 3.1.2**.

The 3D-printed cover and interfacing PCB were added for usability purposes: the 3D-printed PLA cover allows to attach a strap (HTH 833 with H83 hooks, Velcro, United Kingdom) to the sensor in order to maintain it against a subject's skin, as illustrated in **Figure 3B**, while the interfacing PCB gathers the four UART pins from the CO₂ sensor, the two ends of the thermistor, and those of the heating wire into a single eight-pins Microfit connector (0430450812, Molex, United States).

The adhesive itself consists in a disposable laser-cut, double-sided, clinical-grade tape (1567, 3M, United States). For ease of application, a special tooling was developed to accurately align the sensor and the adhesive together—see **Supplementary Material S1**.

## 2.2 Reference tcpCO₂ monitor

In addition to the above-detailed custom-made sensor, a clinical-grade tcpCO₂ monitor (TCM4, Radiometer, Denmark) was also used on the upper deltoid—one of the recommended sites for tcpCO₂ monitoring (SenTec, 2016)—yielding a continuous reference tcpCO₂ reading. The tcpCO₂ sensor itself (tc Sensor 54) was affixed to the skin using an appropriate attachment ring and contact gel, and it was re-membraned and re-calibrated when needed, as per the manufacturer's guidelines. All the accessories used to this end were Radiometer's (Radiometer, 2020).

## 2.3 Sensors positioning

The different sensors and measurement sites chosen in the study are illustrated in **Figure 3A**. All sensors were placed on the subject's left arm: the reference tcpCO₂ monitor was placed on the upper deltoid, as mentioned above, while two custom rate sensors were positioned as follows. The first one was equipped with the Doppler probe, and was attached on the distal side of the upper arm, immediately below the deltoid, at the junction point between the upper part of the biceps, the lower end of the deltoid, and the triceps. The second rate sensor was placed on the dorsal side of the wrist, and did not include a Doppler probe. Both sensors were affixed to the subject's skin by mean of the above-mentioned double-sided adhesive, and secured in place with a Velcro strap. Additionally, the subject's arm laid comfortably onto an arm gutter so that it remains still and relaxed for the whole duration of the experiments.
<figure><img src="image_8.png" /><figcaption>FIGURE 3</figcaption></figure>

(A): outline of the sensors used in this study, and their location. (B): a picture of the sensor with its connection cable assembly and strap, attached on a wrist. (C): simplified model of CO₂ diffusion through the skin inside a closed sensor.

# 3 Methods

## 3.1 Measured metrics

### 3.1.1 Skin CO₂ conductivity and exhalation rate

The rate of diffusion—*a.k.a.* the exhalation rate—of CO₂ through the skin per unit of area—hereafter noted *Q*, of dimension L³·L⁻²·T⁻¹—can be measured by affixing to the skin a cup-like device, which entraps the skin-exhaled CO₂. In this situation, the CO₂ diffusion through the skin can be modelled as presented in Figure 3C. Briefly, the skin is considered as a CO₂-permeable membrane of thickness *w* and diffusivity *D* towards CO₂ (unit of m² s⁻¹). The partial CO₂ pressure inside the sub-cutaneous tissues and inner chamber are tcpCO₂ and *p*<sub>Se</sub>CO₂, respectively, and the sensor area in contact with the skin is *S*<sub>Se</sub>, while its equivalent height and volume are *h*<sub>Se</sub> and *V*<sub>Se</sub>. It can be shown under certain hypotheses (Dervieux et al., 2022) that:

$$ \frac{dp_{Se}CO_2}{dt} = \frac{D}{w \cdot h_{Se}} \cdot (tcpCO_2 - p_{Se}CO_2) \quad (1) $$

Leading to a first-order behaviour for *p*<sub>Se</sub>CO₂, given by:

$$ p_{Se}CO_2(t) = tcpCO_2 \cdot (1 - e^{-t/\tau}) + p_{Se}CO_2(t=0) \cdot e^{-t/\tau}, \quad \text{with } \tau = \frac{h_{Se} \cdot w}{D} \quad (2) $$

And *Q* is then equal to:

$$ Q(t) = \frac{D}{w \cdot P_0} \cdot (tcpCO_2 - p_{Se}CO_2(t=0)) \cdot e^{-t/\tau} \quad (3) $$

wherein *P*₀ is the total atmospheric pressure at measurement site. However, since *Q* depends on both the ambient level of CO₂—via *p*<sub>Se</sub>CO₂(*t* = 0)—and the subject's capnia—via tcpCO₂—we

introduced the skin conductivity—from the thermodynamic or electrical analogy—expressed in m·s⁻¹, and defined as:

$$ K = \frac{D}{w} = \frac{P_0 \cdot Q(t=0)}{tcpCO_2 - p_{Se}CO_2(t=0)} \quad (4) $$

Contrary to *Q*, *K* is an intrinsic property of the skin and is independent of the tcpCO₂/*p*<sub>Se</sub>CO₂ gradient. Additionally, deriving Eq. 2 and evaluating it at *t* = 0 yield:

$$ K = \frac{h_{Se}}{tcpCO_2 - p_{Se}CO_2(t=0)} \cdot \left. \frac{dp_{Se}CO_2}{dt} \right|_{t=0} \quad (5) $$

In practice, *K* was thus measured as follows, choosing arbitrarily *t* = 0 at each temperature change:

- *h*<sub>Se</sub> is known by dividing *V*<sub>Se</sub> by *S*<sub>Se</sub>. The latter is known by construction of the sensor's aluminium body, while *V*<sub>Se</sub> was estimated by filling a clone of the sensor with a low viscosity fluid—pure ethanol—and weighting it.

- the subject's tcpCO₂ was measured using the above-mentioned reference medical grade monitor. The extraction of a single tcpCO₂ reading for each temperature set-point is detailed in Supplementary Material S1.

- *p*<sub>Se</sub>CO₂(*t* = 0) could be measured with a simple reading of the CO₂ sensor.

- $\left. \frac{dp_{Se}CO_2}{dt} \right|_{t=0}$ was estimated by fitting a linear regression on the measured *p*<sub>Se</sub>CO₂. The latter regression was performed on the *p*<sub>Se</sub>CO₂ data starting 3 min after a temperature change—to allow for temperature homogenisation of the different inner parts of the sensor—and up to 18 min after, for a total of 15 min of data. This duration was chosen following a preliminary study performed on ten subjects, which yielded *R*² regression scores above 0.95 for a regression duration above 700 s (about 12 min).

In summary, the diffusion of CO₂ through the skin was quantified by the skin CO₂ conductivity *K*, which was measured
<figure><img src="image_9.png" /><figcaption>FIGURE 4</figcaption></figure>

(A): typical skin perfusion response to local heating, based on data from different sources (Minson et al., 2001; Cracowski et al., 2006; Minson, 2010; Frantz et al., 2012; Roustit and Cracowski, 2012). If we suppose the heat stress to be applied at time origin, three phases are usually observed: ① an onset lag corresponding to i) the heating time required by the sensor to reach its set-point temperature and ii) the time taken by the inner skin layers to reach this temperature and trigger the axon-mediated hyperaemia. ② the axon-mediated hyperaemia which rises quickly and then fades away. ③ the nitric-oxide mediated hyperaemia, whose onset is slower and which slowly fades away if the temperature set-point is not too elevated. (B): perfusion and skin temperature as a function of time as measured at the arm of a test subject.

at five different temperatures, each temperature corresponding to a 18 min measurement window for a total of 90 min of acquisition per subject, as detailed in Section 3.2, below. Additionally, Equation 4 was used to compute the corresponding equivalent Q (t = 0) with $p_{Se}CO_2(t = 0) \approx 0$—i.e. the skin CO₂ exhalation rate in free air as commonly referred to in the literature. Of note, this Q (t = 0) was not observed in practice, since the sensor was left in place—and thus $p_{Se}CO_2(t = 0) \neq 0$ for most temperatures. For the sake of conciseness, in the remainder of this article, the letter Q alone or the mention of “skin CO₂ exhalation rate” without further indications always designate the above-mentioned Q (t = 0). Finally, it should be noted that these Q values were derived mainly as a mean to compare with the existing literature, and that the actual statistical analyses were performed on K—see Section 3.3.

### 3.1.2 Skin blood flow

The skin blood flow—a.k.a. (sub)cutaneous micro-circulation or perfusion—was measured using LDF, and expressed in Perfusion Units (P.U.), a dimensionless arbitrary unit that reflects both the amount and the speed of moving elements—mainly erythrocytes—seen by the Doppler probe (Bonner and Nossal, 1990). When the skin temperature rises, perfusion increases, a phenomenon known as heat-triggered—or thermal—reactive hyperaemia (Minson, 2010), whose dynamics is illustrated in Figure 4A. The respective durations of phases ①–③ were not specified in abscissa since they may vary markedly depending on the heating rate and temperature (Magerl and Treede, 1996; Del Pozzi et al., 2016). To give the reader an order of magnitude, phase ① usually lasts a few minutes, phase ② from 5 up to 10 min, and phase ③ from 30 up to 60 min (Minson et al., 2001; Cracowski et al., 2006; Minson, 2010; Frantz et al., 2012; Roustit and Cracowski, 2012).

Such a behaviour calls for some kind of data processing to yield a single representative perfusion metric for the initial bump, after-bump nadir, and final plateau. In this paper, SkBF₉₀(T) was defined as the 90th percentile of the measured skin blood

flow—SkBF—on an 18 min window at temperature T. This choice was made following preliminary measurements at the arm, an example of which is plotted in Figure 4B. The latter clearly exhibits five perfusion plateaux corresponding to the five temperature set points, and one can also distinguish a small initial bump at the onset of a new temperature, which is especially visible at 35, 38 and 41°C. Due to i) the high variability exhibited by the measured perfusion, especially at high temperature, ii) the fact that the nitric-oxide phase can take up to 30–40 min to establish (Barcroft and Edholm, 1943; Taylor et al., 1984; Minson et al., 2001; Frantz et al., 2012; Del Pozzi et al., 2016), and iii) the fact that each temperature window lasts 18 min, it seemed a good strategy to choose a metric which was more robust than the mean—e.g. the nth percentile—and focused on the very end of the observation window. In this regard, SkBF₉₀(T) seemed to meet the latter requirements and was therefore chosen as perfusion metric. It was additionally normed by the maximal perfusion measured on a given subject—i.e. the SkBF₉₀ value measured at 44°C. Finally, the LDF metric used throughout this study is thus:

$$nSkBF_{90}(T) = \frac{SkBF_{90}(T)}{SkBF_{90}(44^{\circ}\text{C})} \qquad (6)$$

## 3.2 Protocol design

The clinical study was interventional, monocentric, and involved 40 healthy human subjects. Inclusion criteria were an age between 18 and 80, and having given a free and informed consent. Exclusion criteria were the presence of cutaneous lesions at the measurement sites or skin conditions such as dermatitis or psoriasis, and taking vasodilator therapy. The research was approved by the local ethics committee (Comité de protection des personnes Sud-Méditerrannée II, IDRCB ref.: 2020-A02185-38), registered on Clinical Trials (NCT05637138), and it was carried out in accordance with the declaration of Helsinki.
<figure><img src="image_10.png" /><figcaption>**FIGURE 5**

Graphical representation of the data analysis workflow, see the text for further details.</figcaption></figure>

After a preliminary visit during which the subjects were informed of the study modalities and gave consent, all measurements were performed during a single visit, during which the subjects were seated. This visit began with a visual inspection of the measurement sites for detection of cutaneous lesions. The sites were then shaved if needed for a good adhesion of the sensors, using an electric trimmer in order to avoid skin inflammation. The skin was then degreased and cleaned using isopropyl alcohol, and the three above-mentioned sensors were attached to their respective measurement sites. These preliminary steps also allowed for subject acclimation and lasted about 5 min. The measurement itself then began, consisting of five 18 min periods, corresponding to five temperatures for the two rate sensors: NH, 35, 38, 41, and 44°C. At the end of the 90 min measurement period, the sensors were gently peeled off, and the skin was cleaned again. All tcpCO₂ and *p*<sub>Se</sub>CO₂ data were recorded on computers for future analysis, and the room temperature was also recorded using a calibrated thermometer (Testo 735, Testo, Germany).

## 3.3 Data analysis

The data analysis workflow is summarised in Figure 5. Raw data were collected for all 40 subjects at five different temperatures and three metrics were extracted for each *i*th subject/temperature pair: *K* at the arm and wrist, and nSkBF₉₀ at the arm only. For each of those metrics, an ANOVA was performed across all subjects to determine whether their mean values differ significantly between two temperatures. If the ANOVA residuals did not significantly differ from a normal distribution—according to Shapiro-Wilk testing—and if the hypothesis of variance equality between the temperature groups was also verified—according to Bartlett testing—a Tuckey *post hoc* HSD test was then performed. Otherwise, a Kruskal-Wallis test followed by a series of Mann-Whitney U-tests were performed. Additionally, Pearson and Spearman correlation tests were also carried out to study the influence of temperature on the three afore-mentioned metrics (not represented in Figure 5). When applicable, all tests were two-sided and a 5% alpha risk was chosen as significance threshold.

# 4 Results

## 4.1 Demographics and temperatures

The fourty subjects consisted in 24 men and 16 women, aged between 20 and 61 years (mean/median/Standard Deviation (SD): 40/39/13 years). The laboratory temperature was in the 20.1–22.7°C range for the whole duration of the experiments (mean/median/SD: 21.3/21.2/0.7°C). The skin temperature during the non-heated 18 min phase was measured twice, at 10 and 18 min after sensor application, and these two measurements were averaged to yield a single temperature value per subject. The latter was in the 27.1–31.8°C range (mean/median/SD: 29.3/29.2/1.2°C) at the arm, and in the 24.7–33.1°C range (mean/median/SD: 28.6/28.5/1.7°C) at the wrist.

## 4.2 Skin CO₂ conductivity

The linear regressions leading to *K* values—see Section 3.1.1—yielded excellent regression coefficients, with average *R*² values of 0.98 and 0.96 at the arm and wrist, respectively. The resulting skin conductivities are summarised in Figure 6, and show a sharp tendency to increase with an increasing skin temperature. Interestingly, the five upper outlying values at the arm—at all temperatures—and the four upper and below outlying values at the wrist—in the 35–44°C range—belonged each time to a single subject, who exhibited an especially high, or low skin conductivity. However, these two latter subjects were not one and the same person at the arm and at the wrist. The dispersion of skin conductivity values is also glaring with max/min ratios for a given skin temperature/location pair in the 2.9–9.0 range.

A normalisation of the variable *K*—using the change of variable *K* → √*K*—was carried out prior to the ANOVA, resulting in the Shapiro-Wilk and Bartlett tests to be passed (all p-values above 0.05). The ANOVA was significant (p-value below 10⁻¹⁵) and the results of the following Tuckey *post hoc* HSD test are detailed in Table 1. The apparent increase in *K* with an increasing skin temperature seen
<figure><img src="image_11.png" /><figcaption>FIGURE 6</figcaption></figure>

Skin conductivities at the arm and wrist. Each black mark corresponds to a subject/temperature pair, the red circles and texts indicate mean values, and some horizontal jitter was added to the black marks for legibility. The whiskers extend at most to 1.5 times the interquartile range, and descriptive statistics—range, and SD—are provided in Supplementary Material S1—two properties shared by the three box plots of the present paper.

<table><caption>TABLE 1 p-values for the Tuckey HSD post hoc test for differences of the mean K at different temperatures.</caption><thead><tr><th>Arm</th><th>NH</th><th>35°C</th><th>38°C</th><th>41°C</th><th>Wrist</th><th>NH</th><th>35°C</th><th>38°C</th><th>41°C</th></tr></thead><tbody><tr><td>35°C</td><td>&lt; 10<sup>-6</sup>†</td><td>-</td><td>-</td><td>-</td><td>35°C</td><td>&lt; 10<sup>-6</sup>†</td><td>-</td><td>-</td><td>-</td></tr><tr><td>38°C</td><td>&lt; 10<sup>-6</sup>†</td><td>0.72</td><td>-</td><td>-</td><td>38°C</td><td>&lt; 10<sup>-6</sup>†</td><td>0.35</td><td>-</td><td>-</td></tr><tr><td>41°C</td><td>&lt; 10<sup>-6</sup>†</td><td>0.05†</td><td>0.54</td><td>-</td><td>41°C</td><td>&lt; 10<sup>-6</sup>†</td><td>&lt; 10<sup>-3</sup>†</td><td>0.14</td><td>-</td></tr><tr><td>44°C</td><td>&lt; 10<sup>-6</sup>†</td><td>&lt; 10<sup>-3</sup>†</td><td>0.01†</td><td>0.44</td><td>44°C</td><td>&lt; 10<sup>-6</sup>†</td><td>&lt; 10<sup>-6</sup>†</td><td>&lt; 10<sup>-3</sup>†</td><td>0.73</td></tr></tbody></table>

Values with a † are considered significant with a risk $\alpha = 0.05$. (Only the lower-left part of the tables are filled-in for the sake of clarity).

in Figure 6 is hereby confirmed, with significantly different mean K values for most temperature differences except the closest ones—i.e. for the 35/38, 41/38 and 44/41°C pairs.

Pearson and Spearman correlations were both significant (p-values below $10^{-15}$) with correlation coefficients of 0.60 and 0.59 at the arm, respectively, and 0.66 and 0.67 at the wrist. These coefficients indicate a moderate positive influence of the skin temperature on its diffusivity towards CO$_2$.

## 4.3 Skin CO$_2$ exhalation rate

In order to provide a more accessible parameter than K for the reader, as well as to allow direct comparison with existing literature, equivalent initial exhalation rates $Q(t = 0)$ were also computed using Eq. 4. The resulting values are presented in Table 2.

Although our results tend to indicate a higher CO$_2$ exhalation rate at the upper arm than at the wrist—a MANOVA was performed considering the measurement site as independent variable and the five $Q_T$ as dependent variables, and yielded a p-value of 0.01 using Pillai's trace—the size of this effect is moderate, especially in view of the wide Q dispersion.

## 4.4 Laser Doppler Flowmetry

nSkBF$_{90}$ values were computed as described in Section 3.1.2, and are presented in Figure 7. Skin perfusion exhibits a strong increase with temperature, with a tenfold multiplication between the NH basal state and the maximum vasodilation 44°C stage. In particular, a mild heating to a temperature of 38°C already entails a fourfold increase in perfusion. As was the case for skin conductivity, the dispersion of the nSkBF$_{90}$ values is also considerable, with max/min ratios for a given temperature in the 2.2–9.7 range.

Regarding statistical analyses, the normality and variance homogeneity hypotheses could not be verified regardless of the changes of variable performed. A Kruskal-Wallis test was thus performed, followed by a series of Mann-Whitney U-tests, all of which proved significant (all p-values below $10^{-10}$).

Pearson and Spearman correlations were both significant (p-values below $10^{-15}$) with correlation coefficients of 0.90 and 0.96, respectively. These coefficients indicate a strong positive influence of the skin temperature on its perfusion. The fact that Pearson's correlation is below Spearman's is not surprising since the relationship between skin temperature and nSkBF$_{90}$ is strongly non-linear, as emphasised by the sigmoid fit performed in Figure 7.
<table><caption>TABLE 2 Q (t = 0) values, as computed using Eq. 4, expressed in cm³·m⁻²·h⁻¹.</caption><thead><tr><th>Arm</th><th>Mean</th><th>Range</th><th>SD</th></tr></thead><tbody><tr><td>NH</td><td>80.4</td><td>25.1–191.7</td><td>40.3</td></tr><tr><td>35°C</td><td>130.3</td><td>52.4–227.4</td><td>43.5</td></tr><tr><td>38°C</td><td>142.7</td><td>63.0–248.1</td><td>41.4</td></tr><tr><td>41°C</td><td>159.9</td><td>81.9–280.2</td><td>42.8</td></tr><tr><td>44°C</td><td>177.5</td><td>104.6–288.4</td><td>42.5</td></tr></tbody></table><table><thead><tr><th>Wrist</th><th>Mean</th><th>Range</th><th>SD</th></tr></thead><tbody><tr><td>NH</td><td>58.7</td><td>24.2–141.1</td><td>22.7</td></tr><tr><td>35°C</td><td>100.9</td><td>27.7–220.6</td><td>38.3</td></tr><tr><td>38°C</td><td>118.2</td><td>27.1–248.0</td><td>42.8</td></tr><tr><td>41°C</td><td>141.8</td><td>36.0–282.1</td><td>47.2</td></tr><tr><td>44°C</td><td>162.3</td><td>46.5–300.0</td><td>45.9</td></tr></tbody></table>

<figure><img src="image_12.png" /><figcaption>FIGURE 7</figcaption></figure>

Measured nSkBF₉₀ at the arm. Note that since nSkBF₉₀ is normalised with respect to SkBF₉₀ at 44°C, all subjects merge into a single unitary value at this latter temperature. Each black mark corresponds to a subject/temperature pair, the red circles and texts indicate mean values, and the red curve is a least square sigmoid fit. Contrary to Figure 6, no horizontal jitter was added to the data, and the dispersion observed in the 27–32°C range corresponds to the inter-subject variability in non-heated skin temperatures.

# 5 Discussion

The main objectives of this research were to ascertain the influence of skin temperature on *i*) its permeability towards CO₂—through the study of the skin CO₂ conductivity *K*, and exhalation rate *Q*—and *ii*) the skin blood flow—through nSkBF₉₀.

## 5.1 Sensor design

### 5.1.1 Skin contact surface

Although the sensor's aluminium body was precisely machined following the drawing given in Supplementary Material S1, the exact surface area in contact with the skin that participates in gaseous exchange may slightly vary from one subject to another. Indeed, at each hole of the sensor's sole, the skin forms a small dome, whose convexity is essentially function of the mechanical properties of the skin. Yet, those mechanical properties are sex-, moisture-, age-, and temperature-dependant (Salter et al., 1993; Held et al., 2018), thereby introducing small intra- or inter-subject variations in skin contact surface area. Since this area is used to calculate *K* through

$S_{Se}$—see Section 3.1.1—the latter is in turn influenced by these small variations.

While this would likely not change the conclusions of the present study given the order of magnitude of the above-described phenomenon reported in the literature, future research could look into replacing the grid-shaped sole that we used by a thin metal mesh, or a metallic foam. These two latter techniques were for instance implemented by Eletr et al. (1978), McIlroy et al. (1978), and Hansen et al. (1980). However, it must be emphasised that the shape of the sole of a thermally-regulated transcutaneous exhalation rate sensor is essentially a compromise between: *i*) the degree of perforation or porosity of the surface, which should be as high as possible to ensure a large diffusion surface, and *ii*) heat transfer considerations, which call for a plain, dense, surface, with a minimum number of holes, to ensure temperature homogeneity of the skin. Additionally, while the use of a wire mesh, or metallic foam reduces the above-mentioned “dome effect”, it also makes the surface estimation more tedious. Therefore, this aspect of the sensor design should be further investigated to find a satisfactory technical solution which addresses the above concerns.

### 5.1.2 CO₂ sensor choice

The choice of the selected NDIR CO₂ sensor was mainly motivated by its compact form factor and ease of implementation. Additionally, the 5% range was especially adapted for CO₂ diffusion rate measurement. Indeed, tcpCO₂ in healthy subjects is typically in the 35–45 mmHg range (Rithalia et al., 1984), corresponding to 4.6–5.2% of CO₂. Since the CO₂ diffusion rate measurements taking place in the present study were only limited to the first moments of CO₂ diffusion from the skin into the sensor—see Section 3.1.1—measured CO₂ fraction values stayed below 1–2%. The 5% range was thus adapted to our need.

### 5.1.3 Gas tightness

As mentioned in Section 2.1, the gas tightness of the sensor's chamber with respect to ambient air was of paramount importance for the success of the study. Indeed, any leak of inner-chamber CO₂ towards the outer air would subtract from the measured rate of exhalation of CO₂ through the skin, and thus impair the resulting *K* values. During the sensor's design, gas tightness was assessed by sticking the sensor onto a glass plate using the same adhesive as for the human-testing part of the study. The so-obtained glass plate/sensor pair was then put inside a chamber which was successively filled with a 2.5% CO₂/di-nitrogen (N₂) mixture and fresh air. The resulting measurements are presented in Figure 8B, and clearly demonstrate that the greased O-ring alone was not gas tight, while the epoxy sealing was. The Figure 8A also illustrates in
<figure><img src="image_13.png" /><figcaption>FIGURE 8</figcaption></figure>

(A): close up of the section view of Figure 2, showing the two key elements of the gas tight seal: a grease-coated O-ring, and cast epoxy resin. (B): a sealing test comparing the greased O-ring alone, with the greased O-ring and the cast epoxy resin. This test clearly indicates that the greased O-ring alone does not accomplish gas tightness, whereas the cast epoxy resin does. (C): sensor sensitivity towards humidity, showing the onset of condensation onto the gold-plated dome. This condensation drastically reduces the quantity of light reaching the detector, effectively blinding it, which is interpreted as an exceedingly elevated CO₂ concentration.

a cut view how the epoxy seal complements the O-ring. In practice, the grease-coated O-ring acts as a resin-proof sealing that prevents the resin from flowing inside the sensor's chamber during its casting process.

This gas tightness allows CO₂ to accumulate into the sensor chamber until an equilibrium is reached between the subcutaneous tissue and the sensor's chamber—i.e. until $p_{se}CO_2 = tcpCO_2$. While this equilibrium—although it would take several hours given the respective order of magnitudes of Q and $h_{se}$—and the associated CO₂ diffusion process are at the very heart of this study, another undesirable chemical species will also accumulate into the sensor's chamber due to the combined action of diffusion and sweating: water vapour.

While the influence of water vapour on NDIR CO₂ measurements due to the infrared absorbance of water vapour is expected to be negligible given the large gap between CO₂ and water vapour infrared absorption bands (Mranvick, 2023), the onset of condensation onto the reflective part of the sensor—namely the gold-coated reflective dome—can still be an issue. Indeed, the formation of micro-droplets of condensing water onto the latter dome would drastically reduce its reflectance, fooling the sensor into believing that a large amount of CO₂ is present inside the sensor's chamber—a well-known issue in NDIR sensing (Fietzek et al., 2014; Wang et al., 2018). In order to study the influence of condensing humidity levels onto the sensor used in the present research, two experiments were carried out whose outcomes are presented in Figure 8C. The first experiment consisted in placing the sensor on a human thigh at increasing temperatures and waiting for condensation to occur, which happened after 30 min at 44°C. The second experiment consisted in bubbling ambient air (20°C) through pumice stone inside a hot water bath, yielding water-saturated hot air (40°C), which was then flowed onto the un-heated sensor. Even in these unfavourable conditions—i.e. a cold sensor and water-saturated hot air—it took about 20 min to detect the onset of

condensation on CO₂ measurements. Given that the latter onset was particularly sudden and visible on $p_{se}CO_2$ in both experiments, the influence of water vapour condensation on this study was deemed negligible. Indeed, it would be easily detected—were it to happen while measuring a given subject—and the corresponding measurement would be discarded, something which did not happen in practice.

Finally, the reader should bear in mind that the gas tightness of the sensor and the accumulation of humidity underneath it both create a condition called *skin occlusion*. This occlusion, while out of the scope of this paper, has been studied by several authors (Frame et al., 1972; King et al., 1978; Faergemann et al., 1983), who reported much higher CO₂ exhalation rates for long-term—i.e. days—occluded skin, as compared to its basal state. This phenomenon was not investigated in the present study due to the long time scale that it involves, but further research on this topic would be welcome.

### 5.1.4 Sensors positioning

To our knowledge, only three studies compared the influence of the measurement site on the transcutaneous CO₂ diffusion rate in humans: that of Schulze (1943) on twelve subjects (Table 16 *op. cit.*), that of Adamczyk et al. (1966) on one subject and that of Levshankov et al. (1983) on an unspecified number of subjects. The results of the latter two authors are summarised in Figure 9—Schulze indications were difficult to interpret and were thus not illustrated.

The high variability in Adamczyk et al. data—probably caused by the inclusion of only one subject—is glaring, especially when studying left-body/right-body differences. Interesting are the extremely important values reported for the axilla. Those values may be measurement artefacts, or they may be caused by a peculiar behaviour towards CO₂ diffusion of the apocrine glands, which are mainly located in the axilla—see Baker (2019). However, we found no evidence in the literature for or against this hypothesis.
A

<figure><img src="image_14.png" /><figcaption>FIGURE 9</figcaption></figure>

CO₂ diffusion rates through human skin at various sites, expressed in cm³·m⁻²·h⁻¹, based on data from Adamczyk et al. (1966) (A) and Levshankov et al. (1983) (B).

B

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

Alternatively, those elevated values may be caused by the skin temperature, which is much higher at the axilla on resting subjects than at the extremities (Niu et al., 2001; Sund-Levander et al., 2002), since the skin was not heated in their study. At the opposite, the results of Levshenkov et al. are more homogenous concerning left and right body measurements. All in all, and apart from the extreme axilla values, the reported CO₂ diffusion rates exhibit no extreme variations and are of the same order of magnitude, regardless of the measurement site. In this aspect, it thus seems from the limited information at our disposal that no measurement site is far better than the other from the CO₂ diffusion rate perspective.

Consequently, we chose our measurement sites mainly for their ease of access and acceptability, with a view to using these sites for a future wearable tcpCO₂ sensor. In this respect, the dorsal side of the wrist and the upper arm were found to be particularly interesting, as evidenced by the rapidly expanding and widespread use of health-related wristband and armband in the recent years (Al-Eidan et al., 2018; Cosoli et al., 2020; Soon et al., 2020).

## 5.2 Skin CO₂ conductivity

### 5.2.1 Metric choice

It must be emphasised here that in the simplified skin diffusion model introduced in our previous publication (Dervieux et al., 2022) and detailed in Figure 3, the membrane called “skin” does not correspond to an actual physiological membrane. Consequently, its thickness *w* and diffusivity towards CO₂ does not correspond to any physical property that might have been measured on a specific part of the dermis or epidermis. Rather, this “skin” membrane corresponds to a physical modelling of gas transport between the subcutaneous tissue and the outer air. As such, the latter membrane models both the diffusion of CO₂ through the stratum corneum and

the circulation of blood and diffusion of CO₂ in the dermis and subcutaneous tissues.

Moreover, this model also integrates the difference in tcpCO₂ between that measured by the reference Radiometer tcpCO₂ monitor, and that measured at the sensor's location. Indeed, since the reference tcpCO₂ monitor was set to 41°C, it is likely that the tcpCO₂ that we injected in Equation 4 is slightly over-estimated—as per the dilution principle presented in Figure 11—at temperatures below 41°C. Consequently, reported *K*—or *Q*—values below 41°C are likely to be slightly over-estimated. The amplitude of this over-estimation should be in the same order of magnitude as the arterio-venous pCO₂ gradient in resting, healthy subjects—i.e. about 5–15% in the NH-38°C skin temperature range (Kowalchuk et al., 1988; Schneider et al., 2013). However, this state of fact was inevitable since, to the best of our knowledge, no clinical tcpCO₂ monitor working at a temperature below 37°C exists at the time being, and manufacturers recommend using 41–42°C—an injunction that we followed. Future research aiming at extending our work may consider the design of a tcpCO₂ sensor working at low temperature in order to establish the appropriate corrections to the obtained *K* values.

### 5.2.2 Impact on the response time of a future tcpCO₂ sensor

Contrary to perfusion—which increases over elevenfold with skin heating—*K* only doubles from NH to 44°C, and its increase is even smaller between 35–38 and 44°C values. This latter fact is all the more interesting when having in mind the design of a future energy-efficient tcpCO₂ sensor. Indeed, internal studies measuring skin temperature under a wearable device positioned at the upper arm (Bora Band, Biosency, France) on ten healthy subjects revealed that a mean skin temperature of 33.9°C could easily be achieved at the upper arm without additional heating, and that covering
the arm with an additional layer of isolation—i.e. shirt or jumper sleeves—makes it rise even higher to reach 35.1°C.

With such skin temperatures, there is no strong incentive—from a response time point of view—to heat the skin actively any further—i.e. by mean of an external electrical heating system. Indeed, the measured increase of 35% in *K* at the arm from 35 to 44°C—see Figure 6—would result in a decrease in response time of the same magnitude for a given tcpCO₂ sensor, according to the response time model presented in our previous paper (Dervieux et al., 2022). While having a slower sensor may seem like a burning issue for critical care applications, it is not the case for telemonitoring for which long-term tendencies are to be observed over several months (Jang et al., 2021).

Additionally, it should be noted that since the Q values measured in the present study—80–178 cm³·m⁻²·h⁻¹ on average—are in line with that used in our previous publication (Dervieux et al., 2022) for response time calculations—100 cm³·m⁻²·h⁻¹—the afore-proposed sensor thickness of 100 µm for a response time below 10 min remains credible. As a reminder, it was shown in the latter publication that a linear relationship exists between the response time of an equilibrium-based tcpCO₂ sensor, and the volume to surface ratio—i.e. equivalent thickness—of its equilibration medium. Thus, there is essentially a compromise to be made between this thickness, which cannot be infinitely small for technological reasons, and the response time of a so-designed sensor (Dervieux et al., 2022).

Of note, and to the best of our knowledge, there is a lack of clinical guidelines specifying the required response time for tcpCO₂ monitors. Nevertheless, there exists a considerable amount of literature focusing on transcutaneous monitor testing in clinical environments, from which it can be inferred that a typical *in vitro* 90% response time of about 1 min (Bendjelid et al., 2005; Eberhard, 2007) is achievable with current tcpCO₂ monitors. *In vivo* performance reports, on their, part mention an approximatively 10 min initial equilibration time before a first tcpCO₂ reading can be taken (Carter and Banham, 2000; Domingo et al., 2010; Restrepo et al., 2012). Regarding the response time of tcpCO₂ monitors following a sudden change in paCO₂, a lag has been reported in the literature between end-tidal pCO₂—petCO₂—paCO₂, and tcpCO₂, inducing a higher *in vivo* response time than in the ideal *in vitro* case. Reported values for this latter lag fall within the 1–5 min range (Kesten et al., 1991; Carter and Banham, 2000; Cuvelier et al., 2005; Rafl et al., 2018). An overall response time requirement of approximatively 5 min can thus *de facto* be assumed for a tcpCO₂ monitor to meet field expectations. Still, this latter assumption mainly holds for the intensive care of critically ill patients (Mari et al., 2019) and no information exists concerning long-term tcpCO₂ (tele-)monitoring for the obvious reason that the corresponding monitors do not exist yet.

# 5.3 Exhalation Rate

## 5.3.1 An Imperfect Metric

Considering Equation 1, it readily appears that the exhalation rate Q is not constant, and logically depends on the initial $p_{S_e}CO_2$, and of the passing of time. This issue has however been largely ignored by the literature on the topic—see Table 3—and Q has been

considered by most authors as if it had a single constant value. The latter, which has been reported as the CO₂ diffusion rate through the skin, is actually the initial one in free air—i.e. Q (t=0) with $p_{S_e}CO_2(t = 0) \approx 0$—and corresponds to the slopes of the tangents to the $p_{S_e}CO_2$ curves at t = 0 in Figure 10A.

Unfortunately, if measuring Q as illustrated in the latter figure is theoretically feasible at different temperatures, it would also require to remove the sensor at each temperature change, in order to renew the gas inside the inner chamber of the sensor with fresh air. This would in turn require to peel off the sensor from the subject's skin at each temperature change, which would distort the Q measurement, as the skin—and more specifically the *stratum corneum*, its outermost layer—would become thinner and thinner at each sensor replacement—actually, stripping the skin with multiple tape applications is a well-known technique to drastically increase Q (Scheuplein, 1976; Eletr et al., 1978; Greenspan et al., 1981).

Thus, the sensor was left in place in this study, while the temperature was successively changed from NH up to 35, 38, 41, and finally 44°C. This led to measured $p_{S_e}CO_2$ alike that represented in Figure 10B. In that case, using Q as a metric would be impractical, since the $p_{S_e}CO_2$ value at $t_H$ is not null, and Q values would no longer represent initial CO₂ diffusion rates as measured in free air. In practice, the obtained Q values at different temperatures would then not be comparable with each other, each one being measured with a slightly different $p_{S_e}CO_2$ initial value.

## 5.3.2 Comparison with existing literature

Thanks to Equation 4, we can however obtain equivalent Q values in free air from our *K* measurements, and compare them with those of the literature, given in Table 3—at least for the NH case. The values reported in Tables 2 and 3 are of the same magnitude, and the wide amplitudes that we report here—e.g. 25–192 and 24–141 cm³·m⁻²·h⁻¹ at the NH arm and wrist, respectively—are on par with those reported in previous research.

# 5.4 Laser doppler flowmetry

## 5.4.1 Choosing nSkBF₉₀ as a metric

Both inter-subject and inter-site LDF variabilities have often been reported in the literature (Johnson et al., 1984; Cracowski et al., 2006; Minson, 2010; Roustit and Cracowski, 2012; Cracowski and Roustit, 2016; Hodges et al., 2016), and appears to be inherent to this modality of skin blood flow measurement as well as to human physiology in general. Nonetheless, certain guidelines may be followed to obtain the most reproducible results (Cracowski et al., 2006). In particular, when it comes to derive a single explicit LDF metric from a given measurement period—such as a skin-site/sensor-temperature pair, for instance—several techniques have been developed to obtain meaningful results from raw LDF data.

At first, some authors—e.g. Hodges et al. (2016)—prefer to express the skin blood flow as Cutaneous Vascular Conductance (CVC), which is given by the LDF in P.U. or V, divided by the Mean Arterial Pressure (MAP). The CVC is said to be more “physiological” (Cracowski et al., 2006), since an increase in skin blood flow could be caused by an increase in MAP but also by an increase in vascular compliance, for instance. By dividing the LDF-acquired blood flow
<table><caption>TABLE 3 CO₂ exhalation rates through the skin in the literature compared with the present studies.</caption><thead><tr><td>Exhalation rate Q [cm³·m⁻²·h⁻¹]</td><td>Temp. [°C]</td><td>Num. Of Subjects</td><td>Reference</td></tr></thead><tbody><tr><td>25–120</td><td>22–36 (air)</td><td>2</td><td>Shaw and Messer (1930a)</td></tr><tr><td>10–160</td><td>25–37 (air)</td><td>1</td><td>Shaw and Messer (1930b)</td></tr><tr><td><strong>58–169</strong></td><td><strong>26–31 (air)</strong></td><td><strong>38</strong></td><td><strong>Ernstene and Volk (1932), Table 1</strong></td></tr><tr><td>12–143</td><td>23–37 (air)</td><td>1</td><td>Whitehouse et al. (1932)</td></tr><tr><td><strong>32–69</strong></td><td><strong>25–28 (air)</strong></td><td><strong>13</strong></td><td><strong>Schulze (1943), Table 12</strong></td></tr><tr><td>180–2500†</td><td>–</td><td>1</td><td>Adamczyk et al. (1966)</td></tr><tr><td>25–87</td><td>–</td><td>5</td><td>Thiele and Van Kempen (1972)</td></tr><tr><td>11–28</td><td>25–35 (air)</td><td>3</td><td>Frame et al. (1972)</td></tr><tr><td><strong>50</strong></td><td>–</td><td><strong>27</strong></td><td><strong>Levshankov et al. (1983)</strong></td></tr><tr><td><strong>140–221</strong></td><td><strong>36 (skin)</strong></td><td><strong>14</strong></td><td><strong>Eöry (1984)</strong></td></tr><tr><td>25–192</td><td>27–32 (skin)</td><td>40</td><td>This work, NH arm</td></tr><tr><td>24–141</td><td>25–33 (skin)</td><td>40</td><td>This work, NH wrist</td></tr><tr><td>25–288</td><td>27–44 (skin)</td><td>40</td><td>This work, all temp. arm</td></tr><tr><td>24–300</td><td>25–44 (skin)</td><td>40</td><td>This work, all temp. wrist</td></tr></tbody></table>

Past studies with more than ten participants are indicated **in bold** and were used for sample size determination.

<figure><img src="image_17.png" /><figcaption>FIGURE 10</figcaption></figure>

(A): schematic evolution of the pCO₂ inside the sensor's chamber when the skin is heated or not. (B): same as Left, but with a non-heated sensor placed onto the skin for a duration t<sub>fj</sub> before being heated. Note the difference between the two sensing schemes: the left one requires two successive measurements—one heated, the other one non-heated—while the right one consists in a single measurement during which the skin is successively non-heated and then heated. The dashed line on the right represents what would have happened without heating during the whole acquisition, which is equivalent to the solid blue line on the left.

by the MAP, the obtained CVC value is thus in theory more representative of the arteriovenous compliance, a theory supported by several works in haemodynamics (Johnson, 1986; Lautt, 1989; Herring and Paterson, 2018). At the same time, skin blood flow

alone—often abbreviated as SkBF, and either expressed in P.U. or Volts—has been used for several decades (Johnson et al., 1984; Frantz et al., 2012) and remains a good alternative to CVC when MAP is not available.

<footnote>†Axilla measured value, possibly erroneous.</footnote>
Then, once the type of measurement—SkBF or CVC—is chosen, the question that comes next is that of the extraction of a single perfusion metric from a long-lasting acquisition. Indeed, due to the peculiar dynamics of thermal hyperaemia—see Figure 4, above—a simple time-averaging on the whole acquisition duration would make little sense.

To circumvent this issue, several research teams used temporal averaging on manually-set periods of interest in the raw LDF data. The averaging duration that they used depended on the studied phenomenon, with durations of 1–3 min for transient phenomena—i.e. initial bump and after-bump nadir—up to 5–10 min for long-lasting ones—i.e. baseline or maximum perfusion plateau (Minson et al., 2001; Frantz et al., 2012). Other authors, for their parts, chose to average the two to three last minutes of a 10–25 min measurement window at the maximal perfusion value, which obtention is detailed below (Kellogg et al., 2008; Hodges et al., 2016). However, Barcroft and Edholm (1943) and Taylor et al. (1984) mentioned even longer durations for thermal hyperaemia to fully settle following a change in skin temperature—up to 40–60 min. Such a lengthy onset period would result in a total acquisition duration in the 3–5 h range for the five different temperatures involved in the present study. At the opposite, a total experiment duration—including informing the subjects and obtaining their consent—of about 2 h seemed to us to be an acceptable maximum for easily recruiting volunteers. This 2 h duration in turn entails that each temperature window of the present study only lasted 18 min, which may not be enough for the establishment of the nitric-oxide mediated hyperaemia detailed in Figure 4. Thankfully, this 18 min duration is by far long enough for the axon mediated response to take place, and the latter often yields perfusion levels comparable to that reached at the end of the nitric-oxide mediated phase (Kellogg et al., 2008; Minson, 2010; Frantz et al., 2012). Thus, by taking the 90-th percentile of SkBF—see Section 3.1.2—the obtained SkBF₉₀ values are likely to be representative of the SkBF plateau values which would have been observed by increasing the duration of the temperature windows. The latter hypothesis is further confirmed by the similarity between our results and that of the literature, as discussed in the next section.

Finally, it is also common practice to normalise the measured skin blood flow—whether expressed as SkBF or CVC—by its maximum value, often taken after a prolonged (≥15 min) period at an elevated (≥44°C) temperature (Taylor et al., 1984; Vionnet et al., 2014; Hodges et al., 2016) or by direct injection of sodium nitroprusside (Kellogg et al., 2008). Although it has been seldom proposed to normalise the measured values by the baseline blood flow value instead of the maximum one (Magerl and Treede, 1996; Mayrovitz and Leedham, 2001), this is considered bad practice because intra-subject baseline variations can be important even in a temperature controlled room (Bircher et al., 1994; Cracowski et al., 2006).

In this study, CVC was not considered due to the invasiveness of a continuous MAP measurement and SkBF was thus chosen as raw perfusion metric. Then, we proposed to take the 90-th percentile of a given temperature window instead of time averaging. Finally, normalisation by the maximum perfusion value—i.e. the one reached at the end of the 44°C window—was performed, as per literature guidelines.

## 5.4.2 Comparison with the literature

The LDF measurements that were gathered in the present study are consistent with existing literature on the topic. In particular, the sigmoid behaviour observed in Figure 7—revealing a strong onset of hyperaemia in the 35–41°C range—is on par with the observations of Magerl and Treede (1996), Stephens et al. (2001), and Hodges et al. (2016).

## 5.4.3 Impact on the accuracy of a future tcpCO₂ sensor

The fact that the perfusion is doubled at 35°C and quadrupled at 38°C compared to baseline—see Figure 7—is especially encouraging for the development of a future energy-efficient tcpCO₂ sensor, since these temperatures can be easily achieved without—or with minimal—heating, as already discussed in Section 5.2.2 (reaching over 35°C at the upper arm under jumper sleeves). Indeed, since arterialised capillary blood—either obtained by local heating or application of a vasoactive cream—is gaseously close to arterial blood (Zavorsky et al., 2007), it is to be expected that partially arterialised capillary blood obtained by a mild heating—i.e. below 44°C—lies somewhere between venous and arterial blood, from a gaseous content point of view. More specifically, Rooth et al. (1987) hypothesised that the subcutaneous capillary pCO₂—i.e. tcpCO₂—would be a barycentre between venous and arterial pCO₂, as illustrated in Figure 11.

This figure emphasises the fact that—especially for a resting subject—even a mild heating of the skin in the 35–38°C range could be enough to yield a tcpCO₂ only a few mmHg away from the paCO₂. The latter error may be acceptable depending on the clinical application targeted. For example, the Food and Drug Administration (FDA) requires tcpCO₂ monitors to be accurate within 5 mmHg, with an allowed drift of up to 10% of the initial reading over a 1-h period (Food and Drug Administration, 2002).

## 5.5 Sample size

The main objective of the present study was to estimate the mean K value as a function of temperature. The latter mean can be estimated at each temperature T by:

$$ \hat{K}_T = \frac{1}{N} \cdot \sum_{i=1}^{N} K_{S,i} \qquad (7) $$

wherein the index i stands for the ith subject of the study and N stands for its sample size. Contrary to hypothesis testing, for which a sample size may be derived straightforwardly from targeted alpha or beta risks, and some prior knowledge of the data (Ambrosius, 2007, Chap. 19; Chow et al., 2017), sample size determination in the case of an exploratory—or pilot—study is more challenging, with its share of arbitrary decision (Ko and Lim, 2021). Indeed, while a 95% confidence interval can be computed for K_T as:

$$ C.I._{K_T}^{95\%} = [\hat{K}_T - \varepsilon, \hat{K}_T + \varepsilon], \quad \text{and} \quad \varepsilon = -t_{N-1} \left( \frac{0.05}{2} \right) \cdot \frac{s}{\sqrt{N}} \qquad (8) $$

wherein t_{N-1} is the percentile score of a Student distribution with N-1 degrees of freedom, and s is the SD of the sample—i.e. the estimated standard deviation of the population—the value of the latter SD is vastly unknown. In order to estimate an adequate
<figure><img src="image_18.png" /><figcaption>FIGURE 11</figcaption></figure>

Capillary pCO₂ as a function of relative blood flow considering two venous pCO₂ levels: at rest, and while exercising. Relative blood flow values measured in the present study were also added in red with their respective temperature labels. A normal paCO₂ of 40 mmHg (Schneider et al., 2013) was taken for arterial blood, while venous blood levels were set to 46 and 60 mmHg at rest and while exercising, respectively. Of note, while 46 mmHg at rest is generally accepted in the literature (Byrne et al., 2014), the 60 mmHg exercising value was mainly chosen for legibility reasons. Indeed, exercising values may exceed 100 mmHg during heavy exercise, or in case of septic shock (Kowalchuk et al., 1988; Diaztagle Fernández et al., 2017). Modified from Rooth et al. (1987).

sample size for the study at hand—based on an *acceptable* margin of error on $\hat{K}_T$—a prior estimation of *s* is thus needed, and is the object of the upcoming section. Importantly, since *Q*—that is *Q* (*t* = 0) using the above-presented notation—was studied by earlier authors instead of *K*, the following reasoning will be made using the former. This can be done safely since the two values are linked by a proportionality constant—see Equation 4. Of note, Equation 8 holds only if the $K_T$ follow a normal distribution, which had neither been confirmed nor denied in the literature, to the best of our knowledge, but which was verified in the current study—see Section 5.5.3.

### 5.5.1 Literature review

Among the literature studies on the topic of skin CO₂ exhalation rate measurements on human subjects detailed in Table 3, only four of them have been performed on more than ten subjects, and are highlighted in bold in the aforementioned table. Unfortunately, the latter studies are sometimes unclear about *Q* measuring conditions—measurement site and skin temperature, in particular. We did our best not to distort or misinterpret the works of their authors, but what follows is essentially our best interpretation of their writings. These four studies reported—or made it possible to derive from raw data—a $\bar{Q}$ and *s* value for each measurement site investigated. These data can be used to compute 95% confidence intervals on *Q* estimation—as shown in Equation 8—or on *s*, yielding (Ambrosius, 2007, Chap. 4):

$$C.I._{\sigma}^{95\%} = \left[ \frac{s^2 \cdot (N-1)}{\chi_{N-1}^2 \left( \frac{0.05}{2} \right)} ; \frac{s^2 \cdot (N-1)}{\chi_{N-1}^2 \left( 1 - \frac{0.05}{2} \right)} \right] \quad (9)$$

wherein $\chi_{N-1}^2$ is the percentile score of a $\chi^2$ distribution with $N-1$ degrees of freedom. The resulting confidence intervals are reported in Table 4 and tend to indicate a relative uncertainty on *Q* estimation

in the order of 5–30% for relatively small sample sizes—i.e. 13–38 subjects.

### 5.5.2 Chosen sample size

The reported *s* value, as well as the upper and lower bounds of *s* 95% confidence interval were then used to compute the relative uncertainty on *Q* as a function of the number of subjects using Equation 8. While this relative uncertainty decreases when the sample size increases, its reducing rate—as well as the associated uncertainty values—varies wildly depending on the considered data source. Indeed, while a sample size of 20–30 subjects should lead to a relative uncertainty on *Q* in the 5–10% range using Levshankov et al. (1983) data, much larger sample sizes—i.e. 100–150 subjects—would be needed to reach the same level of accuracy using Schulze (1943) or Ernstene and Volk (1932). measurements. In the end, since the works of the latter two authors were much older—1932 and 1943, respectively—than that of Eöry (1984), Levshankov et al. (1983) and, respectively—it was decided to put them aside. The sample size determination was thus grounded only on the works of Eöry (1984), Levshankov et al. (1983), and a sample size of 40 subjects was deemed acceptable, since it should have resulted into a relative uncertainty on *Q* estimation below 10%.

### 5.5.3 Results

Unfortunately, this initial estimation of a 40 subjects cohort proved to be rather optimistic in practice. Indeed, the relative uncertainty on measured *Q* values can be computed using Equation 8, and falls in the 15–32% range, depending on the skin temperature and measurement site. In this aspect, our results are close to those presented by Ernstene and Volk (1932), Schulze (1943), and Eöry (1984) who reported relative uncertainties in the 11–30% range. Yet, the latter authors used less than 40 subjects, and our uncertainty range was thus expected to be narrower than theirs.
<table><caption>TABLE 4 Confidence intervals at the 95% level for σ and Q, computed from s and Q̄ values reported in the literature. Aside from the relative uncertainties—defined as 2·ε/X wherein X is σ or Q—all values are given in cm³·m⁻²·h⁻¹.</caption><thead><tr><th></th><th>Lower bound</th><th>Reported</th><th>Upper bound</th><th>Relative uncertainty (%)</th><th>Measurement site</th><th>Reference</th></tr></thead><tbody><tr><td rowspan="5">σ</td><td>18.7</td><td>22.9</td><td>29.6</td><td>48</td><td>Whole arm</td><td>Ernstene and Volk (1932)</td></tr><tr><td>8.5</td><td>11.9</td><td>19.6</td><td>93</td><td>Abdomen</td><td>Schulze (1943)</td></tr><tr><td>2.65</td><td>3.36</td><td>4.60</td><td rowspan="2">58</td><td>Left forearm</td><td>Levshankov et al. (1983)</td></tr><tr><td>2.41</td><td>3.06</td><td>4.19</td><td>Right forearm</td><td>Levshankov et al. (1983)</td></tr><tr><td>15</td><td>21</td><td>34</td><td rowspan="2">89</td><td>Acupuncture site</td><td>Eöry (1984)</td></tr><tr><td>13</td><td>18</td><td>28</td><td>“adjacent skin area”</td><td>Eöry (1984)</td></tr><tr><td rowspan="6">Q</td><td>113</td><td>120</td><td>128</td><td>13</td><td>Whole arm</td><td>Ernstene and Volk (1932)</td></tr><tr><td>41</td><td>48</td><td>55</td><td>30</td><td>Abdomen</td><td>Schulze (1943)</td></tr><tr><td>48.0</td><td>49.3</td><td>50.7</td><td>5.39</td><td>Left forearm</td><td>Levshankov et al. (1983)</td></tr><tr><td>48.4</td><td>49.6</td><td>50.8</td><td>4.88</td><td>Right forearm</td><td>Levshankov et al. (1983)</td></tr><tr><td>209</td><td>221</td><td>233</td><td>11</td><td>Acupuncture site</td><td>Eöry (1984)</td></tr><tr><td>130</td><td>140</td><td>150</td><td>14</td><td>“adjacent skin area”</td><td>Eöry (1984)</td></tr></tbody></table>

Moreover, the present study also exhibits a higher variability than that of Levshankov et al. (1983)—whose results indicate a relative uncertainty about 5%, see Table 4.

The origin of these discrepancies between literature-driven expectations and the above-presented results is not fully understood at the moment. One possible explanation could be differing measurement sites between the above-mentioned studies and the ones that we chose. Indeed, the data reported by Eöry (1984), Schulze (1943) (Table 16 *op. cit.*.) tend to indicate some degree of variability in the relative uncertainty expected at different sites—in particular when comparing the abdomen and hand in Schulze's data (30 vs. 45%), or the two sites used by Eöry (10 vs. 14%). Thus, it seems plausible that Q variability at the upper arm and wrist is above that reported in earlier studies for differing sites.

Ultimately, we cannot but recommend using larger sample sizes in future studies of a similar nature, considering the significant variability we observed. As a side note related to sample size determination, the normality of the Q distribution was ascertained using a series of Bonferroni-corrected Shapiro-Wilk tests which were non-significant, further justifying the approach presented in Section 5.5. To the best of our knowledge, this is the first report of normality for transcutaneous CO₂ exhalation rates.

## 6 Conclusion

As stated in introduction, the aim of this study was twofold: measuring the influence of skin temperature on the transcutaneous

diffusion of CO₂, and on the skin blood flow. To this end, a custom sensor was designed and used on 40 healthy human subjects at two measurements sites: the upper arm, and the wrist.

Our results indicate comparable behaviours at both sites, with an increasing relationship between temperature on the one hand, and CO₂ exhalation rate, CO₂ conductivity and perfusion on the other hand. These results are encouraging for the development of a future energy-efficient tcpCO₂ sensor for the following reasons:

- Skin conductivity towards CO₂ increases only moderately with an increase in skin temperature, at most doubling from NH to 44°C. Thus, if the response time of the sensor-to-be is not critical—i.e. if a 35% slower response is acceptable compared to the one reachable at maximum skin heating—the latter may not require additional heating. This is especially encouraging in the perspective of building a wearable, battery-operated device.

- Perfusion, for its part, increases strongly with an increase in skin temperature, already doubling from NH to 35°C, and quadrupling from NH to 38°C. This phenomenon is especially interesting since—according to Rooth et al. (1987)—this should bring tcpCO₂ close to paCO₂ even for skin temperatures as low as 35–38°C, which are reachable at the arm without additional heating, given that the latter is covered by warm clothes. However, this latter hypothesis—i.e. the existence of a clinically-satisfying tcpCO₂/paCO₂ correlation in the 35–38°C skin temperature range—is yet to be demonstrated experimentally *in vivo*, which will be the subject of future research.

Additionally, our results highlight the significant variability of transcutaneous CO₂ exhalation rate and conductivity measurements
in human subjects. Hence, we strongly advise future research on the topic to consider large sample sizes—i.e. more than 40 subjects—in order to ensure accurate estimates of the latter metrics. The present study also focuses only on two measurement sites—the upper arm and the wrist—and further investigations at other sites would be welcome. In particular, the remarkably high axilla values reported by some authors is intriguing, and could benefit from a special attention. Of note, the study data—demographics, K, Q, and nSkBF₉₀ values—are provided in **Supplementary Material S2**.

## Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

## Ethics statement

The studies involving humans were approved by the Comité de protection des personnes Sud-Méditerrannée II, IDRCB ref.: 2020-A02185-38. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

## Author contributions

ED: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Visualization, Writing—original draft, Writing—review and editing. FG: Supervision, Validation, Writing—review and editing. WU: Supervision, Validation, Writing—review and editing. M-AG-M: Funding acquisition, Investigation, Supervision, Validation, Writing—review and editing, Project administration. MT: Conceptualization, Funding acquisition, Investigation, Methodology, Supervision, Validation, Writing—review and editing, Project administration.

## Funding

The authors declare financial support was received for the research, authorship, and/or publication of this article. This research was funded by Biosency and the Centre Hospitalier Régional et Universitaire de Brest.

## References

Adamczyk, B., Boerboom, A. J., and Kistemaker, J. (1966). A mass spectrometer for continuous analysis of gaseous compounds excreted by human skin. *J. Appl. Physiology* 21, 1903–1906. doi:10.1152/jappl.1966.21.6.1903

Al-Eidan, R. M., Al-Khalifa, H., and Al-Salman, A. M. (2018). A review of wrist-worn wearable: sensors, models, and challenges. *J. Sensors* 2018, 1–20. doi:10.1155/2018/5853917

Ambrosius, W. T. (2007). *Topics in biostatistics. methods in molecular biology*. Humana Totowa, NJ: Humana Press. Available at: https://link.springer.com/book/10.1007/978-1-59745-530-5

## Acknowledgments

The authors wish to warmly thank Hervé Abitbol and the Tecmatel Engineering team for their design of the adhesive mounting jig, Julien Dolou for his help in the administrative formalities involved in conducting the clinical study, and Alexandre Desvaux, Mikaël Jamet (Cimtech) and Lucien Lecurieux Lafayette (Troisel) for their help in the design and manufacturing of the aluminium sensor body. The authors also acknowledge the following designers/entities for their drawings: Whyakup (computer) and Rflor (power supply) from the Noun Project, Mobius (USB logo) and Mariana Ruiz Villarreal (venous system) from Wikipedia, Pixabay (male anatomy) and the NASA (left arm drawing extracted from Pioneer plaque).

## Conflict of interest

ED works for Biosency, a company which develops wearable cardiopulmonary assessment devices. Although this paper concerns fundamental research in physiology and even if the presented sensor is mainly a research tool—and not the prototype of a commercial sensor—this disclosure may help to avoid the appearance of a conflict of interest.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

## Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2023.1293752/full#supplementary-material

Baker, L. B. (2019). Physiology of sweat gland function: the roles of sweating and sweat composition in human health. *Temperature* 6, 211–259. doi:10.1080/23328940.2019.1632145

Barcroft, H., and Edholm, O. G. (1943). The effect of temperature on blood flow and deep temperature in the human forearm. *J. Physiology* 102, 5–20. doi:10.1113/jphysiol.1943.sp004009

Bendjelid, K., Schütz, N., Stotz, M., Gerard, I., Suter, P., and Romand, J.-A. (2005). Transcutaneous pCO₂ monitoring in critically ill adults: clinical evaluation of a new sensor. *Crit. care Med.* 33, 2203–2206. doi:10.1097/01.CCM.0000181734.26070.26
Bircher, A., de Boer, E. M., Agner, T., Wahlberg, J. E., and Serup, J. (1994). Guidelines for measurement of cutaneous blood flow by laser Doppler flowmetry. A report from the Standardization Group of the European Society of Contact Dermatitis. *Contact Dermat.* 30, 65–72. doi:10.1111/j.1600-0536.1994.tb00565.x

Bonner, R. F., and Nossal, R. (1990). *Principles of laser-Doppler Flowmetry*. Boston, MA: Springer US, 17–45. doi:10.1007/978-1-4757-2083-9_2

Byrne, A. L., Bennett, M., Chatterji, R., Symons, R., Pace, N. L., and Thomas, P. S. (2014). Peripheral venous and arterial blood gas analysis in adults: are they comparable? a systematic review and meta-analysis. *Respirology* 19, 168–175. doi:10.1111/resp.12225

Carter, R., and Banham, S. (2000). Use of transcutaneous oxygen and carbon dioxide tensions for assessing indices of gas exchange during exercise testing. *Respir. Med.* 94, 350–355. doi:10.1053/rmed.1999.0714

Cascales, J. P., Li, X., Roussakis, E., and Evans, C. L. (2022). A patient-ready wearable transcutaneous CO₂ sensor. *Biosensors* 12, 333. doi:10.3390/bios12050333

Chow, S., Shao, J., Wang, H., and Lokhyngina, Y. (2017). *Sample size calculations in clinical research*. London, United Kingdom: Chapman and Hall/CRC Biostatistics Series (Taylor and Francis).

Chung, M., Fortunato, G., and Radacsi, N. (2019). Wearable flexible sweat sensors for healthcare monitoring: a review. *J. R. Soc. Interface* 16, 20190217. doi:10.1098/rsif.2019.0217

Conway, A., Tipton, E., Liu, W.-H., Conway, Z., Soalheira, K., Sutherland, J., et al. (2018). Accuracy and precision of transcutaneous carbon dioxide monitoring: a systematic review and meta-analysis. *Thorax* 74, 157–163. doi:10.1136/thoraxjnl-2017-211466

Cosoli, G., Spinsante, S., and Scalise, L. (2020). Wrist-worn and chest-strap wearable devices: systematic review on accuracy and metricalological characteristics. *Measurement* 159, 107789. doi:10.1016/j.measurement.2020.107789

Cracowski, J.-L., Minson, C. T., Salvat-Melis, M., and Halliwill, J. R. (2006). Methodological issues in the assessment of skin microvascular endothelial function in humans. *Trends Pharmacol. Sci.* 27, 503–508. doi:10.1016/j.tips.2006.07.008

Cracowski, J.-L., and Roustit, M. (2016). Current methods to assess human cutaneous blood flow: an updated focus on laser-based-techniques. *Microcirculation* 23, 337–344. doi:10.1111/micc.12257

Cuvelier, A., Grigoriu, B., Molano, L. C., and Muir, J.-F. (2005). Limitations of transcutaneous carbon dioxide measurements for assessing long-term mechanical ventilation. *Chest* 127, 1744–1748. doi:10.1378/chest.127.5.1744

Dagher, L., Shi, H., Zhao, Y., and Marrouche, N. F. (2020). Wearables in cardiology: here to stay. *Heart rhythm.* 17, 889–895. doi:10.1016/j.hrthm.2020.02.023

Del Pozzi, A. T., Miller, J. T., and Hodges, G. J. (2016). The effect of heating rate on the cutaneous vasomotion responses of forearm and leg skin in humans. *Microvasc. Res.* 105, 77–84. doi:10.1016/j.mvr.2016.01.004

Dervieux, E., Théron, M., and Uhring, W. (2022). Carbon dioxide sensing—biomedical applications to human subjects. *Sensors* 22, 188. doi:10.3390/s22010188

Diaztagle Fernández, J. J., Rodríguez Murcia, J. C., and Sprockel Diaz, J. J. (2017). Venous-to-arterial carbon dioxide difference in the resuscitation of patients with severe sepsis and septic shock: a systematic review. *Med. Intensiva* 41, 401–410. doi:10.1016/j.medin.2017.03.008

Domingo, C., Canturri, E., Moreno, A., Espuelas, H., Vigil, L., and Luján, M. (2010). Optimal clinical time for reliable measurement of transcutaneous co₂ with ear probes: counterbalancing overshoot and the vasodilatation effect. *Sensors* 10, 491–500. doi:10.3390/s100100491

Dunn, J., Runge, R., and Snyder, M. (2018). Wearables and the medical revolution. *Pers. Med.* 15, 429–448. doi:10.2217/pme-2018-0044

Eberhard, P. (2007). The design, use, and results of transcutaneous carbon dioxide analysis: current and future directions. *Anesth. Analgesia* 105, S48–S52. doi:10.1213/01.ane.000027864.16117.f8

Eletr, S., Jimison, H., Ream, A. K., Dolan, W. M., and Rosenthal, M. H. (1978). Cutaneous monitoring of systemic pCO₂ on patients in the respiratory intensive care unit being weaned from the ventilator. *Acta Anaesthesiol. Scand.* 22, 123–127. doi:10.1111/j.1399-6576.1978.tb01406.x

Eöry, A. (1984). *In-vivo skin respiration (CO₂) measurements in the acupuncture loci*. *Acupunct. Electro-Therapeutics Res.* 9, 217–223. doi:10.3727/036012984816714668

Ernstene, A. C., and Volk, M. C. (1932). Cutaneous respiration in man: iv. the rate of carbon dioxide elimination and oxygen absorption in normal subjects. *J. Clin. Investigation* 11, 363–376. doi:10.1172/JCI100418

Faergemann, J., Aly, R., Wilson, D. R., and Maibach, H. I. (1983). Skin occlusion: effect on pityrosporum orbiculare, skin pCO₂, pH, transepidermal water loss, and water content. *Archives Dermatological Res.* 275, 383–387. doi:10.1007/BF000417338

Fietzek, P., Fiedler, B., Steinhoff, T., and Körtzinger, A. (2014). *In situ* quality assessment of a novel underwater pCO₂ sensor based on membrane equilibration and ndir spectrometry. *J. Atmos. Ocean. Technol.* 31, 181–196. doi:10.1175/JTECH-D-1300083.1

Fitzgerald, L. R. (1957). Cutaneous respiration in man. *Physiol. Rev.* 37, 325–336. doi:10.1152/physrev.1957.37.3.325

Food and Drug Administration (2002). *Cutaneous carbon dioxide (PcCO₂) and oxygen (PcO₂) monitors - class II special controls guidance document for industry and FDA. U.S. Department Of Health and Human Services. Tech. rep.*

Frame, G. W., Strauss, W. G., and Maibach, H. I. (1972). Carbon dioxide emission of the human arm and hand. *J. Investigative Dermatology* 59, 155–159. doi:10.1111/1523-1747.ep12625939

Frantz, J., Engelberger, R. P., Liaudet, L., Mazzolai, L., Waeber, B., and François, F. (2012). Desensitization of thermal hyperemia in the skin is reproducible. *Microcirculation* 19, 78–85. doi:10.1111/j.1549-8719.2011.00124.x

Garfan, S., Alamoodi, A. H., Zaidan, B. B., Al-Zobbi, M., Hamid, R. A., Alwan, J. K., et al. (2021). Telehealth utilization during the Covid-19 pandemic: a systematic review. *Comput. Biol. Med.* 138, 104878. doi:10.1016/j.compbiomed.2021.104878

Greenspan, G., Block, A., Haldeman, L., Lindsey, S., and Martin, C. (1981). Transcutaneous noninvasive monitoring of carbon dioxide tension. *Chest* 80, 442–446. doi:10.1378/chest.80.4.442

Hansen, T. N., Sonoda, Y., and McIlroy, M. B. (1980). Transfer of oxygen, nitrogen, and carbon dioxide through normal adult human skin. *J. Appl. Physiology* 49, 438–443. doi:10.1152/jappl.1980.49.3.438

Held, M., Tweer, S., Medved, F., Rothenberger, J., Daigeler, A., and Petersen, W. (2018). Changes in the biomechanical properties of human skin in hyperthermic and hypothermic ranges. *Wounds* 30, 257–262.

Herring, N., and Paterson, D. J. (2018). *Haemodynamics: flow, pressure and resistance*. CRC Press. chap. 8. 121–147. doi:10.1201/9781351107754

Hodges, G. J., McGarr, G. W., Mallette, M. M., Del Pozzi, A. T., and Cheung, S. S. (2016). The contribution of sensory nerves to the onset threshold for cutaneous vasodilatation during gradual local skin heating of the forearm and leg. *Microvasc. Res.* 105, 1–6. doi:10.1016/j.mvr.2015.12.004

Hodgkinson, J., and Tatam, R. P. (2012). Optical gas sensing: a review. *Meas. Sci. Technol.* 24, 012004. doi:10.1088/0957-0233/24/1/012004

Huttmann, S. E., Windisch, W., and Store, J. H. (2014). Techniques for the measurement and monitoring of carbon dioxide in the blood. *Ann. Am. Thorac. Soc.* 11, 645–652. doi:10.1513/AnnalsATS.201311-387FR

Jang, S., Kim, Y., and Cho, W.-K. (2021). A systematic review and meta-analysis of telemonitoring interventions on severe COPD exacerbations. *Int. J. Environ. Res.* 18 (13), 6757. doi:10.3390/ijerph18136757

Johnson, J. M., Taylor, W. F., Shepherd, A. P., and Park, M. K. (1984). Laser-Doppler measurement of skin blood flow: comparison with plethysmography. *J. Appl. Physiology* 56, 798–803. doi:10.1152/jappl.1984.56.3.798

Johnson, P. C. (1986). Autoregulation of blood flow. *Circulation Res.* 59, 483–495. doi:10.1161/01 RES.59.5.483

Kellogg, D. L., Zhao, J. L., and Wu, Y. (2008). Endothelial nitric oxide synthase control mechanisms in the cutaneous vasculature of humans *in vivo*. *Am. J. Physiology-Heart Circulatory Physiology* 295, H123–H129. doi:10.1152/ajpheart.00082.2008

Kesten, S., Chapman, K. R., and Rebuck, A. S. (1991). Response characteristics of a dual transcutaneous oxygen/carbon dioxide monitoring system. *Chest* 99, 1211–1215. doi:10.1378/chest.99.5.1211

King, R., Rl, C., Maibach, H. I., Jh, G., Ml, W., and Jc, J. (1978). The effect of occlusion on carbon dioxide emission from human skin. *Acta dermato-venenologica* 58 (2), 135–138. doi:10.2340/000155558135138

Ko, M. J., and Lim, C.-Y. (2021). General considerations for sample size estimation in animal study. *Korean J. Anesthesiol.* 74, 23–29. doi:10.4097/kja.20662

Koch, G. (1965). Comparison of carbon dioxide tension, pH and standard bicarbonate in capillary blood and in arterial blood with special respect to relations in patients with impaired cardiovascular and pulmonary function and during exercise. *Scand. J. Clin. Laboratory Investigation* 17, 223–229. doi:10.1080/00365516509075339

Kowalchuk, J. M., Heigenhauser, G. J., Lindinger, M. I., Sutton, J. R., and Jones, N. L. (1988). Factors influencing hydrogen ion concentration in muscle after intense exercise. *J. Appl. Physiology* 65, 2080–2089. doi:10.1152/jappl.1988.65.5.2080

Kruse, C., Pesek, B., Anderson, M., Brennan, K., and Comfort, H. (2019). Telemonitoring to manage chronic obstructive pulmonary disease: systematic literature review. *JMIR Med. Inf. 7*, e11496. doi:10.2196/11496

Lautt, W. (1989). Resistance or conductance for expression of arterial vascular tone. *Microvasc. Res.* 37, 230–236. doi:10.1016/0026-2862(89)90040-x

Lermuzeaux, M., Meric, H., Sauneuf, B., Girard, S., Normand, H., Lofaso, F., et al. (2016). Superiority of transcutaneous CO₂ over end-tidal CO₂ measurement for monitoring respiratory failure in nonintubated patients: a pilot study. *J. Crit. Care* 31, 150–156. doi:10.1016/j.jcrc.2015.09.014

Levshankov, A. I., Pushkina, M. A., Slutskia, M. E., and Uvarov, B. S. (1983). Determination of local gas exchange on the body surface by the method of mass spectrometry. *Meditsinskaia tekhnika* 1, 21–26. doi:10.1007/BF00560505
Magerl, W., and Treede, R. D. (1996). Heat-evoked vasodilatation in human hairy skin: axon reflexes due to low-level activity of nociceptive afferents. *J. Physiology* 497, 837–848. doi:10.1113/jphysiol.1996.sp021814

Mari, A., Nougue, H., Mateo, J., Vallet, B., and Vallée, F. (2019). Transcutaneous pCO₂ monitoring in critically ill patients: update and perspectives. *J. Thorac. Dis.* 11, S1558-S1567. doi:10.21037/jtd.2019.04.64

Mayrovitz, H. N., and Leedham, J. A. (2001). Laser-Doppler imaging of forearm skin: perfusion features and dependence of the biological zero on heat-induced hyperemia. *Microvasc. Res.* 62, 74–78. doi:10.1006/mvre.2001.2314

McIlroy, M. B., Simbruner, G., and Sonoda, Y. (1978). Transcutaneous blood gas measurements using a mass spectrometer. *Acta Anaesthesiol. Scand.* 22, 128–130. doi:10.1111/j.1399-6576.1978.tb01407.x

Minson, C. T. (2010). Thermal provocation to evaluate microvascular reactivity in human skin. *J. Appl. Physiology* 109, 1239–1246. doi:10.1152/japplphysiol.00414.2010

Minson, C. T., Berry, L. T., and Joyner, M. J. (2001). Nitric oxide and neurally mediated regulation of skin blood flow during local heating. *J. Appl. Physiology* 91, 1619–1626. doi:10.1152/jappl.2001.91.4.1619

Mranvick, (2023). What is the influence of humidity on NDIR CO₂ measurements? Physics Stack Exchange. Available at: https://physics.stackexchange.com/q/788799 (version: 2023-11-22.

Nanji, A. A., and Whitlow, K. J. (1984). Is it necessary to transport arterial blood samples on ice for pH and gas analysis? *Can. Anaesth. Soc.* J. 31, 568–571. doi:10.1007/BF03009545

Niu, H. H., Lui, P. W., Hu, J. S., Ting, C. K., Yin, Y. C., Lo, Y. L., et al. (2001). Thermal symmetry of skin temperature: normative data of normal subjects in Taiwan. *Chin. Med. J.* 64, 459–468.

Radiometer (2020). TCM CombiM, continuous blood gas monitoring. Available at: https://pdfhost.io/v/DaeEjvgEq_radiometer CombiM.

Rafl, J., Kulhanek, F., Kudrna, P., Ort, V., and Roubik, K. (2018). Response time of indirectly accessed gas exchange depends on measurement method. *Biomed. Eng.* 63, 647–655. doi:10.1515/bmt-2017-0070

Restrepo, R. D., Hirst, K. R., Wittnebel, L., and Wettstein, R. (2012). AARC clinical practice guideline: transcutaneous monitoring of carbon dioxide and oxygen: 2012. *Respir. Care* 57, 1955–1962. doi:10.4187/rescare.02011

Rithalia, S. V. S., Clutton-Brock, T. H., and Tinker, J. (1984). Characteristics of transcutaneous carbon dioxide tension monitors in normal adults and critically ill patients. *Intensive Care Med.* 10, 149–153. doi:10.1007/BF00265805

Rooth, G., Ewald, U. W., and Caligara, F. (1987). Transcutaneous pO₂ and pCO₂ monitoring at 37°C. cutaneous pO₂ and pCO₂. *Adv. Exp. Med. Biol.* 220, 23–32.

Roustit, M., and Cracowski, J.-L. (2012). Non-invasive assessment of skin microvascular function in humans: an insight into methods. *Microcirculation* 19, 47–64. doi:10.1111/j.1549-8719.2011.00129.x

Salter, D., McArthur, H. C., Crosse, J., and Dickens, A. (1993). Skin mechanics measured in vivo using torsion: a new and accurate model more sensitive to age, sex and moisturizing treatment. *Int. J. Cosmet. Sci.* 15, 200–218. doi:10.1111/j.1467-2494.1993.tb00075.x

Scheer, B. V., Perel, A., and Pfeiffer, U. J. (2002). Clinical review: complications and risk factors of peripheral arterial catheters used for haemodynamic monitoring in anaesthesia and intensive care medicine. *Crit. Care* 6, 199–204. doi:10.1186/cc1489

Scheuplein, R. J. (1976). Permeability of the skin: a review of major concepts and some new developments. *J. Investigative Dermatology* 67, 672–676. doi:10.1111/1523-1747.ep12544513

Schneider, A. G., Eastwood, G. M., Bellomo, R., Bailey, M., Lipcsey, M., Pilcher, D., et al. (2013). Arterial carbon dioxide tension and outcome in patients admitted to the intensive care unit after cardiac arrest. *Resuscitation* 84, 927–934. doi:10.1016/j.resuscitation.2013.02.014

Schulze, W. (1943). Untersuchungen über die Alkaliempfindlichkeit, das Alkalineutralisationsvermögen und die Kohlensäureabgabe der Haut. *Arch. für Dermatol. Syph.* 185, 93–161. doi:10.1007/BF02714173

SenTec (2016). *Sentec v-sign brochure, rf-007857-bm*.

Severinghaus, J. W., and Astrup, P. B. (1986). History of blood gas analysis. III. Carbon dioxide tension. *J. Clin. Monit.* 2, 60–73. doi:10.1007/BF01619178

Shaw, L. A., and Messer, A. C. (1930b). Cutaneous respiration in man: II. The effect of temperature and of relative humidity upon the rate of carbon dioxide elimination and oxygen absorption. *Am. J. Physiology-Legacy Content* 95, 13–19. doi:10.1152/ajplegacy.1930.95.1.13

Shaw, L. A., Messer, A. C., and Weiss, S. (1930a). Cutaneous respiration in man: I. factors affecting the rate of carbon dioxide elimination and oxygen absorption. *Am. J. Physiology-Legacy Content* 95, 107–118. doi:10.1152/ajplegacy.1929.90.1.107

Soon, S., Savarsdottir, H., Downey, C., and Jayne, D. G. (2020). Wearable devices for remote vital signs monitoring in the outpatient setting: an overview of the field. *BMJ Innov.* 6, 55–71. doi:10.1136/bmjinnov-2019-000354

Stephens, D. P., Charkoudian, N., Benevento, J. M., Johnson, J. M., and Saumet, J. L. (2001). The influence of topical capsaicin on the local thermal control of skin blood flow in humans. *Am. J. Physiology-Regulatory, Integr. Comp. Physiology* 281, R894–R901. doi:10.1152/ajpregu.2001.281.3.R894

Steventon, A., Bardley, M., Billings, J., Dixon, J., Doll, H., Hirani, S., et al. (2012). Effect of telehealth on use of secondary care and mortality: findings from the whole system demonstrator cluster randomised trial. *BMJ* 344, e3874. doi:10.1136/bmj.e3874

Sund-Levander, M., Forsberg, C., and Wahren, L. K. (2002). Normal oral, rectal, ticpanic and axillary body temperature in adult men and women: a systematic literature review. *Scand. J. Caring Sci.* 16, 122–128. doi:10.1046/j.1471-6712.2002.00069.x

Taylor, W. E., Johnson, J. M., O'Leary, D., and Park, M. K. (1984). Effect of high local temperature on reflex cutaneous vasodilation. *J. Appl. Physiology* 57, 191–196. doi:10.1152/jappl.1984.57.1.191

Thiele, F. A. J., and Van Kempen, L. H. J. (1972). A micro method for measuring the carbon dioxide release by small skin areas. *Br. J. Dermatology* 86, 463–471. doi:10.1111/j.1365-2133.1972.tb16098.x

Tufan, T. B., and Guler, U. (2022). “A miniaturized transcutaneous carbon dioxide monitor based on dual lifetime referencing,” in 2022 IEEE Biomedical Circuits and Systems Conference (BioCAS), 144–148. doi:10.1109/BioCAS54905.2022.9948600

Vionnet, J., Calero-Romero, I., Heim, A., Rotaru, C., Engelberger, R. P., Dischl, B., et al. (2014). No major impact of skin aging on the response of skin blood flow to a submaximal local thermal stimulus. *Microcirculation* 21, 730–737. doi:10.1111/micc.12154

Wagner, P. D. (2015). The physiological basis of pulmonary gas exchange: implications for clinical interpretation of arterial blood gases. *Eur. Respir. J.* 45, 227–243. doi:10.1183/09031936.00039214

Wang, J. N., Xue, Q. S., Lin, G. Y., and Ma, Q. J. (2018). Mid-infrared carbon dioxide sensor with wireless and anti-condensation capability for use in greenhouses. *Spectrosc. Lett.* 51, 266–273. doi:10.1080/00387010.2018.1468785

Whitehouse, A. G. R., Hancock, W., and Haldane, J. S. (1932). The osmotic passage of water and gases through the human skin. *Proc. R. Soc. Lond. Ser. B, Contain. Pap. A Biol. Character* 111, 412–429. doi:10.1098/rspb.1932.0065

Wimberley, P. D., Grønlund Pedersen, K., Olsson, J., and Siggaard-Andersen, O. (1985). Transcutaneous carbon dioxide and oxygen tension measured at different temperatures in healthy adults. *Clin. Chem.* 31, 1611–1615. doi:10.1093/clinchem/31.10.1611

Yetisen, A. K., Martinez-Hurtado, J. L., Únal, B., Khademhosseini, A., and Butt, H. (2018). Wearables in medicine. *Adv. Mater.* 30, 1706910. doi:10.1002/adma.201706910

Yun, J. E., Park, J.-E., Park, H.-Y., Lee, H.-Y., and Park, D.-A. (2018). Comparative effectiveness of telemonitoring versus usual care for heart failure: a systematic review and meta-analysis. *J. Cardiac Fail.* 24, 19–28. doi:10.1016/j.cardfail.2017.09.006

Zavorsky, G. S., Cao, J., Mayo, N. E., Gabbay, R., and Murias, J. M. (2007). Arterial versus capillary blood gases: a meta-analysis. *Respir. Physiology Neurobiol.* 155, 268–279. doi:10.1016/j.resp.2006.07.002