Original Research ARTICLE
Gas analyzer's drift leads to systematic error in maximal oxygen uptake and maximal respiratory exchange ratio determination
- 1Studies, Research and Sports Medicine Center, Government of Navarre, Pamplona, Spain
- 2Laboratory of Performance, Sport-Occupational Activities-Biology-Association, Lyon-Chassieu, France
The aim was to examine the drift in the measurements of fractional concentration of oxygen (FO2) and carbon dioxide (FCO2) of a Nafion-using metabolic cart during incremental maximal exercise in 18 young and 12 elderly males, and to propose a way in which the drift can be corrected. The drift was verified by comparing the pre-test calibration values with the immediate post-test verification values of the calibration gases. The system demonstrated an average downscale drift (P < 0.001) in FO2 and FCO2 of −0.18% and −0.05%, respectively. Compared with measured values, corrected average maximal oxygen uptakevalues were 5–6% lower (P < 0.001) whereas corrected maximal respiratory exchange ratio values were 8–9% higher (P < 0.001). The drift was not due to an electronic instability in the analyzers because it was reverted after 20 min of recovery from the end of the exercise. The drift may be related to an incomplete removal of water vapor from the expired gas during transit through the Nafion conducting tube. These data demonstrate the importance of checking FO2 and FCO2 values by regular pre-test calibrations and post-test verifications, and also the importance of correcting a possible shift immediately after exercise.
Maximal oxygen uptake () is defined as the highest rate at which oxygen can be taken up and utilized by the body during exercise. In laboratory settings, is commonly measured during incremental exercise to exhaustion, during which expired air is analyzed. The key variables needed to calculate are the ventilator flow and the inspired and the expired fractional concentrations of oxygen (FIO2 and FEO2, respectively) and carbon dioxide (FICO2 and FECO2, respectively) (Hodges et al., 2005; Gore et al., 2013).
One of the main potential sources of error in the calculation of using automated systems is related to the stability of FEO2 and FECO2 measurements, because the electronic oxygen (O2) and carbon dioxide (CO2) analyzers are prone to drift over time (Winter, 2012; Gore et al., 2013). To our knowledge, there is surprisingly relatively little information available on the stability of O2 and CO2 analyzing systems over time during incremental exercise (Hodges et al., 2005; Salier Eriksson et al., 2012). In virtually all the publications that have measured , the authors have mentioned performing a pre-test calibration. As it has been pointed out in a recent Editorial (Winter, 2012), in the majority of these studies it is rare to see, however, equivalent post-test verifications. For instance, after reviewing more than 50 studies measuring published between 1973 and 2012, we have found only 8 studies (~16%) in which the authors mentioned that the analyzers' drift at the completion of exercise was assessed (Wilmore et al., 1976; Armstrong and Costill, 1985; Prieur et al., 1998; McLaughlin et al., 2001; Rietjens et al., 2001; Day et al., 2003; Gore et al., 2003; Bowen et al., 2012). Only 4 of these 8 studies reported the average numerical drift values in O2% and CO2% (Wilmore et al., 1976; Armstrong and Costill, 1985; Prieur et al., 1998; Rietjens et al., 2001), which ranged from 0.02 to 0.22%. These reported drift values, according to the equations governing gas concentrations (Beaver et al., 1973; Wasserman et al., 1994a), would have caused an error in up to 8–9% in standard laboratory conditions (~20°C of temperature, ~40% of relative humidity and ~720 mmHg of barometric pressure). Furthermore, none of these 4 studies gave any criterion for the maximum drift error that can be accepted. It is still unknown whether the drift magnitude is related to some physical or physiological exercise variables and how long any particular drift remains after the end of exercise. It is also unclear how the drift readings should be adjusted or corrected to overcome the inaccuracy due to the drift (Winter, 2012).
Clearly, it seems that the process of post-test verification tends to be overlooked and there is insufficient data available on how stable specific gas analysis systems are during exercise conditions (Atkinson et al., 2005; Salier Eriksson et al., 2012). This issue may be particularly relevant in several modern analyzers, in which the exhaled gas is not dried but is equilibrated with the laboratory environment by the use of a length of semi-permeable Nafion tubing (Medbø et al., 2002; Larsson et al., 2004). The purpose of the present study was, therefore, to examine the drift over time of a Nafion-using O2 and CO2 analyzing system during maximal incremental exercise in experienced athletes and elderly sedentary males. By including sedentary elderly and young athletic subjects, as well as short and long-duration exercise protocols, a large range of metabolic responses and exercise durations were examined and the influence of the drift on oxygen uptake (), CO2 output () and respiratory exchange ratio (RER) assessment was determined. This study also proposed a way in which the error might be reduced.
Materials and Methods
Eighteen male amateur athletes (young group) and twelve older men (elderly group) volunteered to participate in the study. Athletes were recruited from various regional Sports Federations (Swimming, Athletics, Basketball, Basque-Ball, Paddle Tennis, Mountaineering and Climbing, Karate, Taekwondo, Judo, and Boxing). Athletes' mean (± SD) age, height, body mass,and percentage of body fat were 22 ± 6 years, 182 ± 7 cm, 79.3 ± 8.3 kg and 10.4 ± 3.1%, respectively. Participants in the elderly group were recruited from a Physical Activity Program for persons over 55. Mean (± SD) age, height, body mass, and percentage of body fat of the participants constituting the elderly group were 69 ± 6 years, 167 ± 7 cm, 85.9 ± 13.3 kg and 27.3 ± 4.3%, respectively. A detailed medical history was taken on the day of the study. No subject reported a history of abnormal dyspnea on exertion or of angina.
Written informed consent was obtained from all volunteers prior to their participation. The study was approved by the Institutional Review Committee of the Instituto Navarro del Deporte y Jueventud (Government of Navarre, Spain), according to the requirements of the Declaration of Helsinki.
Two different maximal incremental exercise protocols, with different exercise stage duration, were used for each population to examine whether the drift is influenced by the duration of the test. All testing sessions within each group were performed at the same time of the day in an air-controlled and well ventilated laboratory with a volume of 1121 m3. Young and elderly individuals reported to the laboratory at least 2 h after their last meal and having refrained from caffeine, alcohol, and strenuous or non-habitual exercise for 24 h before testing. Participants were habituated to the exercise testing equipment and procedures, as they were previously tested in the same laboratory using similar testing procedures.
Young Exercise Trials
Participants were habituated to the exercise testing equipment and procedures, as they were previously tested in the same laboratory using similar testing procedures. was determined by a continuous maximal graded exercise test while sitting on a mechanically braked cycle-ergometer (Monark, Ergomedic 839-E, Varberg, Sweden). The exercise started at 20 W and the load was increased by 20 W every 2 min until volitional exhaustion. This exercise protocol was designed to reach volitional exhaustion within 23–33 min. It has been shown that relatively short (8–12 min) or long (~30 min) protocols do not affect attainment of in highly motivated athletes (Gore et al., 2013). Participants maintained a constant cycling pedaling cadence of 60 rpm. Exhaustion was defined as the subject not being able to maintain the required pedaling cadence, despite vigorous verbal encouragement during the last min of exercise.
Elderly Exercise Trials
was determined by a continuous incremental maximal exercise test on a treadmill ergometer (Kuntaväline, Hyper Treadmill 2040, Finland). The exercise test started at 5.5 km· h−1, after one min the speed was increased to 6.1 km· h−1 for another min, and thereafter grade was increased 1.1% every min until volitional exertion. Exhaustion was defined as the subjects not being able to maintain the required exercise intensity or they wished to stop.
At least two of the following criteria had to be met to determine in both groups (American College of Sports Medicine, 2009): (1) no increase in despite increased workload, defined as a increment of less than 120 ml· min−1 per stage in the young group or a increment of less than 1.75 ml· kg−1· min−1 per stage in the elderly group. This criterion implies that any increment lower than 50% of the metabolic demand of these protocols' stages was accepted as a plateau (Taylor et al., 1955). (2) A maximal respiratory exchange ratio (RERmax) greater than 1.10 (Robergs et al., 2010); (3) peak blood lactate concentration greater than 8 mmol·L−1, and (4) peak heart rate exceeding 90% of age predicted maximum (220-age). Heart rate (Polar Electro Oy, RS800CX, Kempele, Finland) was monitored throughout the exercise in both groups. Capillary blood samples from hyperemic earlobe were obtained at rest, on completion of the trial and at the 1st and 3rd min of recovery. After cleaning and puncturing, the single-use enzyme-coated electrode test strip was directly filled by a 5 μl whole-blood sample and blood lactate concentration was amperometrically determined (Arkray KDK Corporation, Lactate Pro LT-1710, Shiga, Japan).
Collection of Respiratory Gases
Participants were fitted with an appropriately sized mouth and nasal breathing mask (Series 7930, Hans Rudolph, Kansas City, MO, USA) adjusted with a headgear (Vacu-Med, Ventura, CA, USA). Metabolic data was continuously collected using a Vista Mini-CPX (Vacu-Med, Silver Edition 17670, Ventura, CA, USA) computer-integrated metabolic system. The Vista Mini-CPX is a high precision mass flowmeter instrument composed of a turbine flow sensor and O2 and CO2 analyzers designed to measure the flow of the exhaled gases and the concentrations in the O2 and CO2 gases on-line. At the start of each test, room temperature (RT), barometric pressure (PB), and relative room humidity (RH) were measured (Precision Barometer, Lufft, Fellbach, Germany) and these data were entered manually into the computer. The environmental laboratory conditions were kept within the recommended values (18–23°C with a relative humidity lower than 70%) (Gore et al., 2013) by means of a heating system.
Minute expired ventilation () is calculated by a signal generated by the volume transducer of the turbine flow sensor. FEO2 is measured at RH through a disposable galvanic fuel cell (Teledyne Analytical Instruments, R-22MED Oxygen Sensor, Industry, CA, USA). FECO2 is measured at RH through a nondispersive infrared system (Servomex, Ir1507 CO2 infrared transducer, Crowborough, UK). According to the manufacturer, the CO2 and O2 analyzers have zero drift (< 1.5 Torr in 1 h for the CO2 analyzer and 0.3% a week at constant temperature for the O2 analyzer) and their response times are 90 to 130 ms (CO2 analyzer) and 5 s (O2 analyzer). This time delay is automatically assessed and the length of the airline is taken into account according to the manufacturer's specifications. From these measurements the metabolic cart's computer calculates the mass flow of (in liters per minute), (in liters per minute), and the ratio of to (RER) with an accuracy (according to the manufacturer) of ±1% in measures of FEO2 and FECO2, of ±2% in measures of , and of ±3% in measures of and .
This metabolic system uses a proportional sampling approach in the process of mixing the exhaled gases. Thus, the flow rate of this sampling is closely related to the flow of exhalation at ~0.5% of its rate, and directs the exhaled gases in three steps into the O2 and CO2 gas analyzers connected in parallel: (1) through a capillary tube, into a miniature mixing chamber, (2) through a built-in Nafion gas dryer humidifier conducting 180 tube (29 cm long × 1 mm inner diameter), and (3) through a capillary tube system with the same configuration of sampling tube length, diameter and pump flow rate for both analyzers. The Nafion tube is a semi-permeable membrane to water vapor made of copolymer of tetrafluoroethylene (Teflon®) and perfluoro-3,6-dioxa-4-methyl-7-octene-sulfonic acid, highly selective in the removal of water from the vapor phase. The Nafion tube allows water vapor to pass in and out of the tube by absorption and conveys the exhaled gases to the gas analyzers once an equilibrium is reached with the ambient humidity (Macfarlane, 2001). According to the Nafion manufacturer, during exercise the water vapor tension of the aspirated gas sample (relative humidity ~100%) (Bageant, 1976; Macfarlane, 2001; Atkinson et al., 2005) is reduced in milliseconds to the level of RH of the laboratory environment (~27%, Table 1) by moving the water through the Nafion membrane wall and evaporating it very quickly into the surrounding air. Conversely, the typically dry calibration gas is humidified by the Nafion tubing to the level of RH. This system provides a constant value of water vapor tension of the exhaled and calibration gases just prior to the entry of the samples into the gas analyzers. The Nafion tube was replaced at least every 3 years according to the manufacturer. All tests were carried out within the 18 months following the last Nafion tube replacement.
The metabolic measurement software supplied with the analyzer (Vacu-Med, TurboFit 5, Ventura, CA, USA) was set to report mean metabolic data over a 30 s time period and to adjust the volume of the expired air to standard conditions (STPD) for temperature (0°C), pressure (760 mmHg), and dry (absence of water vapor). was defined as the highest 30-s value averaged over two consecutive readings, and its time-corresponding values of , and RER were considered as , , and RERmax, respectively.
Pre-test Calibration and Post-test Verification Processes
The instrument was warmed up for at least 2 h prior to every exercise test to minimize any possible electrical drift. Calibration of the O2 and CO2 analyzers was performed immediately prior to every test using two-point calibration with two precision-analyzed gas mixtures. One calibration point was room air (O2: 20.93%; CO2: 0.00%) Non-hygroscopic soda lime CO2 absorbent (Vacu-Med, Ventura, CA, USA) was used for maximum precision of ambient CO2 measurement. Thus, fractional concentrations of room air were assumed to be 20.93% O2 and 0.00% CO2. The second point was a high-precision certified calibration tank gas containing 15.05% O2, 5.99% CO2 and balanced nitrogen. This high-precision gas was determined gravimetrically, was obtained from a reliable gas supplier (Praxair, Madrid, Spain) and had a claimed accuracy of ±0.02%. Turbine flow calibration was determined using a high-precision 3-L calibration syringe (Vacu-Med, Calibringe 1092, Ventura, CA, USA), in a five-pump series. A series of complete pumps of the syringe and of gas calibrations were repeated until the difference between the current and the previous calibration was less than 0.05 L for volume and less than 0.02% for O2 and CO2. When the calibration process was finished, the gas sample line was connected to the subject's mask.
Within 15 s of the completion of each exercise trial the sample line was removed from the connection to the face mask/turbine and the after trial verification of FO2, FCO2 and turbine flow measurements was performed. Both calibration gases (room air and tank gas) were run through the metabolic system to check for the drift of the analyzer over the course of the measurement period. Verification readings of the calibration gases and the flow sensor were noted down and compared with the calibration references.
Correction of Metabolic Data
Post-test verifications readings were used to correct the metabolic data measured by the Vista Mini-CPX. Corrected () in STPD condition was calculated as follows:
where “3” was the volume (L) of the syringe used to calibrate the flow sensor, “” was the minute ventilation (L· min−1) in STPD condition measured by the metabolic cart, “Cal” was the calibration readout (L) recorded before the exercise and “Ver” was the verification readout (L) recorded after the exercise.
Correction of FO2 is illustrated in Figure 1. During the pre-test calibration process we adjusted the gain settings of the span potentiometers to the corresponding voltage outputs, so that readings of O2% (tank gas: y1 = 15.05%; room air: y2 = 20.93%) equaled real O2% (x1 = y1; x2 = y2). The equation of the pre-test calibration regression line is therefore:
Figure 1. Correction of fractional concentrations of oxygen. x1 and x2, true tank (x1 = 15.05%) and room air (x2 = 20.93%) fractions of oxygen; y1 and y2, fractions of tank (y1 = 15.05%) and room oxygen (y2 = 20.93%) read by the analyzer during the pre-test calibration process when the true tank (x1) and room air (x2) gases were aspirated by the analyzers; y′1 and y′2, fractions of tank (y′1), and room oxygen (y′2) read by the oxygen analyzer during the post-test verification process when true tank (x1) and room air (x2) gases were aspirated by the analyzers.
During the post-exercise verification process we used the same pre-test calibration gases (x1 = 15.05% O2; x2 = 20.93% O2), but the %O2 values read (y′1 and y′2) were different from the O2% read during the pre-test calibration process. In this case, the equation of the post-test verification regression line is:
For a given value of (y′n) measured at during exercise, we can calculate the corresponding value of x (xn) from the equation of the post-test verification line (Y = A′X + B′) as follows:
Therefore, the corrected FEO2 value at (y′n) in the pre-test calibration line (Y = X) is:
FECO2 was corrected using this same procedure. Once the corrected FEO2, FECO2, and were obtained, formulas provided by the manufacturer [see Beaver et al. (1973) or Wasserman et al. (1994a) for further detail] were employed to correct and as follows:
where FiCO2, FiO2, and FiN2 are fractions of inspired carbon dioxide, oxygen and nitrogen respectively, FECO2, FEO2, and FeN2 are fractions of corrected expired carbon dioxide, oxygen and nitrogen respectively, is the corrected minute ventilation (L·min−1) in STPD condition, and HF is the humidity factor defined as:
where PB is the barometric pressure (mmHg) and PH2O is the pressure of water (mmHg) at room temperature (RT) and humidity (RH). Standard tables provided by the manufacturer, also presented by Wasserman et al. (1994b), were used to determine PH2O.
The metabolic system calculates FiN2 and FeN2 using the next two formulas:
where it is assumed that FiN2 is constant and FeN2 is the remaining fractional gas of HF, FEO2 and FECO2.
All corrections were performed off-line using specific routines developed in a commercial software package (The MathWorks Inc., MATLAB R2008a, Natick, MA, USA).
To check the stability of the analyzers, each exercise test was pair-matched on duration, time of the day and number of pre-test calibrations and post-test verifications assessed, with a non-exercise trial, accounting for a total of 30 non-exercise trials (one per subject). These non-exercise trials consisted of performing the identical calibration and verification processes of the gas analyzers over the same time interval to that used during each exercise trial. Between the calibrations and verifications, the metabolic system worked throughout but no subject was connected to the metabolic cart. No flow or volume measures were recorded.
The pattern of change in FO2 and FCO2 during the first 30 min of recovery after the completion of the exercise trials, and after disconnecting the gas sample line from the mask, was investigated immediately after 9 exercise trials. These recovery trials consisted of performing the post-test verifications of the gas analyzers within 15 s of the completion of each exercise trial, but also at 3, 5, 10, 15, 20, and 30 min of recovery from each exercise trial.
Standard statistical methods were used for the calculation of means, standard deviations (SD), standard errors of the estimates (SEE), and confidence intervals (CI). Data were analyzed using parametric statistics following confirmation of normality, homoscedasticity, and when appropriate sphericity. Gas measure readings after the trials (verification readings) were compared with the concentrations of the standard calibration gases (calibration readings) using two-tailed one-sample Student's t-tests. Two-tailed Student's paired t-tests were used to analyze differences between verification readings of the exercise trials and their paired non-exercise trials, as well as between the non-corrected (measured) and corrected values of the respiratory parameters. Respiratory values of the elderly and young groups were compared by two-tailed independent samples t-tests, with Levene's tests used to assess equality of variances. Relationships between variables of interest were assessed by linear regression analyses. Pearson product-moment correlation coefficients (r) were used to indicate the magnitude and direction of each linear relationship. The slopes of the regression lines in elderly and young groups were compared using analysis of covariance (ANCOVA). Differences between pre- and post-test values in FO2 and FCO2 during the recovery period were analyzed using one factor ANOVA with repeated measures. When significance was found, Student's t-test with Bonferroni correction for multiple comparisons was used to locate the significance. Significance was set at P < 0.05. Statistical analyses were performed using SPSS 17.0 (SPSS Inc., Chicago, USA). Data in the text, tables and figures are reported as mean ± SD.
As designed, the duration of the young cycling exercise trials (26:53 ± 3 min) was higher (P < 0.001) than the duration of the elderly treadmill exercise trials (9:29 ± 3 min). Maximal power output reached by young athletes was 294 ± 34 W (3.74 ± 0.54 W· kg−1). Maximal grade attained by elderly individuals at 6.1 km· h−1 was 8.5 ± 3.9%. Young athletes attained significantly higher (P < 0.001) peak heart rate and peak blood lactate concentration values (195 ± 11 b· min−1 and 10.3 ± 2.2 mmol·L−1) compared to elderly individuals (144 ± 24 b· min−1 and 6.6 ± 1.7 mmol·L−1, respectively).
The pre-test calibration and post-test verification values of FO2 and FCO2 of the room air and tank gases assessed within 15 s of the completion of each exercise trial in the whole group of subjects are presented in Table 2. The system showed a downscale drift (P < 0.001) in FO2 and FCO2 from pre- to post-test values in the exercise trials. Mean absolute differences between pre- and post-test values were −0.18% (room air) and −0.14% (tank gas) in O2 and 0.00% (room air) and −0.05% (tank gas) in CO2. Expressed as a percentage of the average pre-test calibration values, the magnitude of the downscale drift was similar (~0.9%) in both analyzers. There was no statistical difference (P = 0.08; 95% CI: −0.00 to 0.01 L) in the registered air volumes between post-test verification (2.99 ± 0.01 L) and pre-test calibration values (2.99 ± 0.01 L). This means that the calibration factor for ventilation volume was essentially constant throughout the test period.
Table 2. Calibration (pre-test) and verification (post-test) readings of the exercise and non-exercise trials.
Figure 2 presents the relationships in the total sample between the individual values of and the individual post-test verification values of FO2 and FCO2 of both calibration gases (room air and tank gas). Regression analyses indicated significant negative correlations between and post-test verification values of room air FO2 in the total sample (r = −0.48; P = 0.007; SEE = 0.056%; 95% CI: 20.77–20.91%) and in the young group (r = −0.49; P = 0.03; SEE = 0.057%; 95% CI: 20.75–21.04%). The gradients of the rest of the relationships presented in Figure 2 were not different from zero (P < 0.05). According to the ANCOVA results, the slopes of the regression lines were not different among groups (P > 0.05). No significant relationships were observed between test duration and post-test verification values of FO2 and FCO2 (P > 0.05). No other relevant significance was found between respiratory parameters and post-trial readings.
Figure 2. Relationships between the individual values of maximal minute ventilation and the individual post-test verification values of fractional concentrations of oxygen (FO2; A,B) and carbon dioxide (FCO2; C,D) when both calibration gases (room air and tank gas) were run through the metabolic system after maximal exercise. Open circles: elderly sedentary subjects. Filled circles: young athletes.
During the non-exercise trials, the drift over time in the electronic gas analysis system was minimal because FO2 and FCO2 remained very stable throughout the time (Table 2). The highest individual difference in the post-test verification during the non-exercise trials was only of 0.03% in FO2 and of 0.02% in FCO2.
Measured and Corrected Respiratory Values
Measured FEO2 values reached at during exercise by the young and elderly groups were 17.39 ± 0.29% and 17.18 ± 0.52%, respectively. When these values were corrected with the proposed correction equation, the corresponding FEO2 values at were 17.54 ± 0.32% and 17.32 ± 0.53% for the young and elderly groups, respectively. Measured FECO2 values reached at by the young and elderly groups were 3.77 ± 0.28% and 4.00 ± 0.51%, respectively. When these values were corrected, the corresponding FECO2 values at were 3.80 ± 0.29% and 4.03 ± 0.50% for the young and elderly groups, respectively. Inasmuch as no drift was observed in the calibration factor for ventilation volume during exercise, there were no differences between corrected and measured values of in any of the groups (P > 0.05). Average was 88% higher (P < 0.001; 95% CI 49 to 80 L·min−1) in the young group compared with the elderly group (137 vs. 73 L·min−1).
The corrected FEO2 and FECO2 values resulted in systematic significant changes in and values. Measured values reached by the young and elderly groups were 4.93 ± 0.57 L·min−1 and 2.82 ± 0.61% L·min−1, respectively. When these values were corrected, the average values (5.07 ± 0.59 L·min−1 and 2.98 ± 0.62 L·min−1 for the young and elderly groups respectively) were 3–5% higher (P < 0.001) than the corresponding measured values. The measured average values in the young and elderly groups were 4.64 ± 0.56 L·min−1 and 2.62 ± 0.50 L·min−1, respectively. Corrected average values (4.35 ± 0.46 and 2.50 ± 0.47 L·min−1 for the young and elderly groups, respectively) were 5–6% lower (P < 0.001) than the corresponding measured values. The individual overestimation of the measured values ranged from 0.3 to 11%. Figure 3A shows the average and the individual measured and corrected values, expressed relative to kilogram of body mass, in the young and elderly subjects. Average corrected values were 3.6 ml·kg−1·min−1 (young) and 1.4 ml·kg−1·min−1 (elderly) lower (P < 0.001) than the average measured values. In every subject, the corrected value was lower than the measured value.
Figure 3. Individual measured (M) and corrected (C) values of maximal oxygen uptake (A) and maximal respiratory exchange ratio (B) in the young and elderly groups. Maximal oxygen uptake is expressed in ml·kg·min−1. The bars indicate mean values. ***Significant difference between the corrected and the corresponding measured values (P < 0.001). †Significant difference between groups (P < 0.05).
Figure 3B shows the average and the individual measured and corrected RERmax values in the young and elderly subjects. The average measured RERmax values were 1.06 ± 0.05 in the young group and 1.07 ± 0.05 in the elderly group. When these values were corrected, the average RERmax values (1.16 ± 0.06 and 1.15 ± 0.06 for the young and elderly groups, respectively) were 8–9% higher (P < 0.001) than the corresponding measured values. In every subject, the corrected RERmax value was higher than the measured value.
When the measured values were taken into account, 14 out of the 18 young subjects (78%) and 9 out of the 12 old subjects (75%) satisfied at least two of the criteria established to verify attainment of . When the RERmax and the values were corrected, the ratio of the subjects who met these criteria increased to 89 and 83% in the young and elderly groups respectively.
Figure 4 shows the average and individual FO2 changes observed in 9 subjects when the post-test verification process was repeated several times during the first 30 min of recovery after the completion of the exercise trials, and after disconnecting the gas sample line from the subjects' mask. During the first 5 min of recovery the average FO2 remained similar to the significantly diminished values (P < 0.001) read immediately after the end of the exercise trials. From that time on, the FO2 reading values increased progressively and linearly over the time. The disappearance of the drift was completed after 20 min of recovery, although at this time the average FO2 readings still tended to be slightly lower than the pre-test calibration values (P = 0.20). Similar patterns were observed for the time course of FCO2 changes (data not shown).
Figure 4. Average and individual time course of the fractional concentrations of oxygen (FO2) during recovery after maximal exercise. The post-test verification values read by the gas analyzers using the room air (A) and the tank gas (B) were measured at 20 s, 3, 5, 10, 15, 20, and 30 min of recovery. The number of observations made at each time-point was 9. ***Significantly different from pre-test (P < 0.001); **Significantly different from pre-test (P < 0.01); *Significantly different form pre-test (P < 0.05).
The main finding of this study is that the pre-test calibration and the post-test verification values of O2 and CO2 demonstrated a downscale drift in the O2 and CO2 readings. The drift was observed in all the exercise tests and was higher than the absolute accuracy of at least ±0.03% (Gore et al., 2013) and ±0.05% (Jones, 1988) that laboratories should strive to attain for electronic O2 and CO2 analyzers. This indicates that the present metabolic system systematically underestimates FEO2 and FECO2 values during maximal exercise.
Several potential sources of error, working separately or together, could explain the FEO2 and FECO2 downscale drifts during maximal exercise (Robergs et al., 2010). One potential source of error may be due to an electrical instability in the analyzers over time (Kannagi et al., 1983). Evidence of this mechanism, however, has not been provided. When a series of calibrations were assessed during the non-exercise pair-matched trials without any subject being connected to the metabolic system, the O2 and CO2 readings remained unchanged over the course of the period (Table 2). This suggests that no base-line drift of the analyzers occurred due to an electronic error, indicating that the analyzers were electrically stable for a long period of time.
The most likely factor explaining the reduction in O2 and CO2 percentages may be associated with how water vapor is handled in the aspirated gas by the analyser. The metabolic system used in the present study sends the exhaled gas to the O2 and CO2 gas analyzers through a built-in Nafion gas dryer humidifier conducting tube. This tube provides a constant value of water vapor tension of the exhaled and the calibration gases just prior to the entry of the samples into the gas analyzers. It is possible that the observed downscale drift could be partly explained by an incomplete removal of the water vapor tension of the aspirated gas by the analyser to equilibrate the partial water vapor pressure (PH2O) into and out of the Nafion tube wall. Since the O2 and CO2 analyzers are partial pressure sensors that measure gas fractions of the total gas volume including water vapor, and they are sensitive to the presence of water vapor molecules, the passage of excessive water vapor to the gas analyzers could raise the PH2O of the sample. A rise in PH2O would reduce O2 and CO2 fractions by the factor [(PB − PH2O excess)·(PB)−1] or [(1 − FH20)] (Gore et al., 2013) and the analyzer would read lower concentration values (Auchincloss et al., 1970). The observation that the O2 and CO2 drifts were almost completely reversed in a few min after exercise by simply disconnecting the sampling line from the flow-meter and the subject's mask, and by flushing the system with room air (Figure 4), supports the notion that some failure in the drying process occurred during exercise.
Under the assumption of an incomplete removal of the water vapor, it is possible to estimate the average extra amount of PH2O at a given temperature that was not removed by the Nafion tube to equilibrate the aspirated gas by the analyzers to the level of ambient humidity during exercise. This can be calculated from the average drift values observed in the O2 (from 20.93 to 20.75% and from 15.05 to 14.91%) and CO2 analyzers (from 5.99 to 5.94%) (Table 2) using the following formula (Bageant, 1976; Gore et al., 2013):
where read O2% is the oxygen percentage read during the post-test verification, true O2% is the oxygen percentage read during the pre-test calibration, and PB is ambient barometric pressure (in our case: ~724 mmHg).
In that case, the estimated average PH2O that could not be removed was 6.2 mmHg (range 0.7–11 mmHg) for the O2 calibration with room air, 6.7 mmHg (range: 0–11 mmHg) for the O2 calibration with the tank, and 6.0 mmHg (range: 2.4–12.1 mmHg) for the CO2 calibration with the tank. Inasmuch as the PH2O of the exhaled gas leaving the body is ~47 mmHg (on the basis of ~100% of relative humidity, at body temperature) (Bageant, 1976), an incomplete average removal of around 6.3 mmHg of water vapor corresponds to ~13% of excess in relative humidity (6.3 · 100 · 47−1) that cannot be cleared from the circuit, with individual values ranging from 2 to 24%.
The reason why the Nafion tube could not fully equilibrate the gas being conveyed to the analyzers with the ambient humidity is unknown. However it can be related to:
(1) A saturation process that reduces active surface area in the Nafion tubing. It is known that some saturation process occurs in the Nafion tubing since the wall of the tubing always retains some residual water, because the sulphonic acid groups within the Nafion polymer will never give up all their water (Mauritz and Moore, 2004). When the dryer becomes progressively physically wet over time, a failure to dry occurs. This failure to dry may be more relevant when the exhaled air flow is high and, therefore, when the aspirated gas sample's flow rate (0.5% of the exhaled flow rate) and its water vapor content are high. For example, in the young exercise trials the amount of water vapor content to be removed out of the Nafion tube can be 16 times higher at maximal exercise (exhaled flow gas: 190 L·min−1; aspired gas: 950 ml·min−1) than at rest (exhaled flow gas: 12 L·min−1; aspired gas: 60 ml·min−1). This is in agreement with the significant linear negative correlation observed in this study between and the magnitude of the drift in FO2 (Figure 2). This strongly suggests that the higher the and the amount of water vapor to be removed, the higher the absolute magnitude of the drift.
(2) The inability of the system to maintain a very low water pressure outside, in the air surrounding the Nafion tube wall. An excess of condensate water vapor may be surrounding the Nafion tube as a consequence of the release of the excess of moisture out of the tube. This process may be more pronounced when the Nafion tube is located inside the metabolic measurement cart, such as in the metabolic system used in this study. In such a case, the fan of the metabolic cart cannot remove this excess water vapor condensed inside the metabolic cart.
(3) Factors like the accumulation of sweat, saliva, foreign bodies and condensation generated by the subject can enter the internal lumen of the sampling line; a portion of exhaled air is drawn and, therefore, a change in the resistance of the delivery tubing or in the gas sampling rate can occur. This could contribute to a decrease in the gas flow rate and pressure in the sampling tube, leading to irregular results (Atkinson et al., 2005; Gore et al., 2013).
The present results support the above theoretical possibilities that cause an incomplete removal of water vapor of the aspired gas transported from the mouth to the analyzers before the gas enters the analyzers. This would explain the significant downscale drift in the O2 and CO2 analyzers that occurs during continuous measurement during human exercise.
The practical question to consider is the influence of the analyzers' drifts on . The correction used for the difference between the pre- and post-test conditions indicated that the corrected values were on average 3.6 ml·kg−1·min−1 (young subjects) and 1.4 ml·kg−1·min−1 (older subjects) lower than those of the measured values. When expressed relative to the individual values, the average difference between the measured and the corrected values was similar (5–6%) in the young and the elderly subjects. This suggests that, in relative terms, there is a systematic and considerable overestimation in the measurement of that is uniform over a full range of values regardless of exercise duration. The average technological error of 5–6% may be considered unacceptable because it is larger than the ±0.5 to ±3% (technological error) or the ±2.2 to ±4% (technological plus biological variation) accuracy standards accepted for the precision of measurement by most certifying organizations that supervise the accreditation process of the metabolic systems (American Thoracic Society, 1987; Gore et al., 2013). The present results may explain, at least partly, the reason why a measurement error of 5% in between laboratories and metabolic systems is nowadays a difficult goal to achieve, owing to the combined technical error and the biological variation (Hodges et al., 2005).
A question raised is the comparison of the corrected and RERmax data with published values. When compared to the measured values, the average corrected values in the elderly group (29.6 ml·kg−1·min−1) and the average corrected -to cycling work rate values in the young group (14.7 ml O2·W−1) are lower than the measured values (31.0 ml·kg−1·min−1 and 15.7 ml O2·W−1), and compare favorably with those estimated for the elderly group (29.5 ml·kg−1·min−1) using the formula of the American College of Sports Medicine (2009) and with the average ratio (14.1 ml O2·W−1; range: 12.1–18.6) reported by other investigators using different metabolic systems during long duration (15–27 min) incremental maximal cycling tests (Pollock et al., 1982; Armstrong and Costill, 1985; Storer et al., 1990; American College of Sports Medicine, 2009; Bowen et al., 2012; Petot et al., 2012; Adami et al., 2013). The average corrected RERmax was 9% higher than the measured RERmax in the young group (1.16 vs. 1.06) and 8% higher in the elderly group (1.15 vs. 1.07). When RERmax values were not corrected, only 17% of the young and 25% of the elderly subjects reached a RERmax greater than 1.10, the most widely used secondary criterion to verify attainment of (Howley et al., 1995; American College of Sports Medicine, 2009). A major effect of correcting the RERmax values was that the ratio of the subjects reaching a RERmax greater than 1.10 was increased to 72% in the young group and to 75% in the elderly group. The difference between the corrected and measured RERmax values suggests that some inconsistencies and failures found in several studies to satisfy RERmax criterion for achievement of may be largely due to an artifact related to technological error (Bowen et al., 2012). This indicates that correction of and RERmax values, on the basis of the FEO2 and FECO2 drifts observed, produced more reasonable and satisfactory values than the measured ones.
This study has several limitations. The major drawback comes from the fact that we did not corroborate the validity of the correction method suggested. There is also a lack of consensus on which method is the most appropriate to assess the reliability and validity of measures (Salier Eriksson et al., 2012). The conventional Douglas bag procedure has been regarded as the gold standard method to validate metabolic measurement systems (McLaughlin et al., 2001; Rietjens et al., 2001). This method remains, however, very limited (Salier Eriksson et al., 2012). In any case, in close agreement with our corrected values, Medbø et al. (2002) and Larsson et al. (2004) found that a commercial metabolic system (Metamax II), utilizing a built-in Nafion conducting tube, significantly overestimated by 4–13% and underestimated RER by 6% compared to the Douglas bag method. However, other validation studies have produced more varied results (Versteeg and Kippersluis, 1989; Bassett et al., 2001; McLaughlin et al., 2001). An alternative method to validate and RER measures is to use a metabolic calibrator system. However, the external validity of such a test is limited since it often uses dry gases and does not involve challenging factors such as humidified gases and irregular breathing patterns (Macfarlane, 2001). In the absence of a reliable gold standard method, the rationale for the analyzer's drift correction method used in this study is that the time point at , which was reached close to the end of the test, is close to the time point at which the post-test verification was undertaken (within 15 sec of the end of each test). It seems, therefore, justifiable to remove and correct the variations observed in FEO2 and FECO2 at by adjusting the analyzer pre-exercise base-line values to the post-exercise verification values.
Another limitation of this study is that we used a single metabolic system. Therefore, the generalizability of our findings is constrained to the metabolic cart and the analyzers used. However, four studies reporting the average numerical downscale drifts in O2 immediately after exercise using other metabolic systems have found values ranging from −0.02 to −0.22% (Wilmore et al., 1976; Armstrong and Costill, 1985; Prieur et al., 1998; Rietjens et al., 2001). This indicates that an absolute downscale drift also occurs in other metabolic systems. If the main source of the error is related to the built-in Nafion gas dryer humidifier conducting tube, a lower error (or none) should occur when gas fractions are measured as fractions of dry gas, when ambient relative humidity is higher than in the present study (e.g., 60%) or when the condensate water vapor surrounding the Nafion tube is more efficiently removed. A wider study is needed to extend the present findings to the wide metabolic systems' population.
Conclusion and Perspectives
In conclusion, the present experiment indicates that, under controlled laboratory conditions, a physiologically significant downscale drift in FO2 and FCO2 was observed over time at the end of maximal exercise in elderly sedentary and young athletes using a metabolic cart equipped with a built-in Nafion conducting tube. The most likely explanation for the drift is an accumulation of excess water vapor in the sample line which could not be completely removed during transit through the Nafion conducting tube. The correction method proposed indicates that ignoring the effects of the drift would induce an average overestimation of 5–6% and a RERmax underestimation of 8–9%, with errors ranging up to 11–12% () and up to 15–16% (RERmax). Therefore, ignoring the drift can have an important influence on the accurate calculation of these variables. The disagreement between the measured and the corrected and RERmax values observed in this particular metabolic system is not acceptable to test athletes, to prescribe exercise intensities, to calculate the fat oxidation rate from RER values, or to use the respiratory values for some other clinical purposes, such as to guide treatment in patients with chronic heart failure (Bowen et al., 2012), to enter in cardiac transplantation listing, to indicate the health status or to predict prognosis and mortality (Myers et al., 2002; Mehra et al., 2006). The implications of the present study point to the necessity to check FO2 and FCO2 values by carefully calibrating the pre-test calibration gases and verifying a possible shift immediately after exercise, as well as to correct the respiratory data in situations where the drift in O2 and CO2 analyzers occurs. Special care must be taken in studies where a Nafion conducting tube is used. Further research in this area is certainly warranted to establish valid correction factors for each device.
EG, IG, and JE conceived and designed the experiments; EG, IG, and JA contributed to the acquisition and analysis of the data; EG, IG, JE, and JA interpreted the data; EG wrote the first draft; EG, IG, JE, and JA critically reviewed and edited the drafts; all authors approved the final version of the manuscript.
Conflict of Interest Statement
The 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.
We gratefully thank Irene Madariaga for her professional assistance in exercise testing.
Adami, A., Sivieri, A., Moia, C., Perini, R., and Ferretti, G. (2013). Effects of step duration in incremental ramp protocols on peak power and maximal oxygen consumption. Eur. J. Appl. Physiol. 113, 2647–2653. doi: 10.1007/s00421-013-2705-9
Armstrong, L. E., and Costill, D. L. (1985). Variability of respiration and metabolism: responses to submaximal cycling and running. Res. Q. Exerc. Sport 56, 93–96. doi: 10.1080/02701367.1985.10608441
Atkinson, G., Davison, R. C., and Nevill, A. M. (2005). Performance characteristics of gas analysis systems: what we know and what we need to know. Int. J. Sports Med. 26(Suppl. 1), S2–S10. doi: 10.1055/s-2004-830505
Bassett, D. R. Jr., Howley, E. T., Thompson, D. L., King, G. A., Strath, S. J., McLaughlin, J. E., et al. (2001). Validity of inspiratory and expiratory methods of measuring gas exchange with a computerized system. J. Appl Physiol. (1985). 91, 218–224. Available online at: http://jap.physiology.org/content/91/1/218.full-text.pdf+html
Bowen, T. S., Cannon, D. T., Begg, G., Baliga, V., Witte, K. K., and Rossiter, H. B. (2012). A novel cardiopulmonary exercise test protocol and criterion to determine maximal oxygen uptake in chronic heart failure. J. Appl. Physiol. (1985). 113, 451–458. doi: 10.1152/japplphysiol.01416.2011
Day, J. R., Rossiter, H. B., Coats, E. M., Skasick, A., and Whipp, B. J. (2003). The maximally attainable VO2 during exercise in humans: the peak vs. maximum issue. J. Appl. Physiol. (1985). 95, 1901–1907. doi: 10.1152/japplphysiol.00024.2003
Gore, C. J., Clark, R. J., Shipp, N. J., Van der Ploeg, G. E., and Withers, R. T. (2003). CPX/D underestimates VO(2) in athletes compared with an automated Douglas bag system. Med. Sci. Sports Exerc. 35, 1341–1347. doi: 10.1249/01.MSS.0000079045.86512.C5
Gore, C. J., Tanner, R. K., Fuller, K. L., and Stanef, F. (2013). “Determination of maximal oxygen consumption (VO2max),” in Physiological Tests for Elite Athletes, eds R. K. Tanner and C. J. Gore (Champaign, IL: Human Kinetics), 103–122.
Hodges, L. D., Brodie, D. A., and Bromley, P. D. (2005). Validity and reliability of selected commercially available metabolic analyzer systems. Scand. J. Med. Sci. Sports 15, 271–279. doi: 10.1111/j.1600-0838.2005.00477.x
Kannagi, T., Bruce, R. A., Hossack, K. F., Chang, K., Kusumi, F., and Trimble, S. (1983). An evaluation of the Beckman Metabolic Cart for measuring ventilation and aerobic requirements during exercise. J. Cardiac. Rehabil. 3, 38–53.
Larsson, P. U., Wadell, K. M., Jakobsson, E. J., Burlin, L. U., and Henriksson-Larsén, K. B. (2004). Validation of the MetaMax II portable metabolic measurement system. Int. J. Sports Med. 25, 115–123. doi: 10.1055/s-2004-819953
McLaughlin, J. E., King, G. A., Howley, E. T., Bassett, D. R. Jr., and Ainsworth, B. E. (2001). Validation of the COSMED K4 b2 portable metabolic system. Int. J. Sports Med. 22, 280–284. doi: 10.1055/s-2001-13816
Medbø, J. I., Mamen, A., Welde, B., von Heimburg, E., and Stokke, R. (2002). Examination of the Metamax I and II oxygen analysers during exercise studies in the laboratory. Scand. J. Clin. Lab. Invest. 62, 585–598. doi: 10.1080/003655102764654321
Mehra, M. R., Kobashigawa, J., Starling, R., Russell, S., Uber, P. A., Parameshwar, J., et al. (2006). Listing criteria for heart transplantation: international Society for Heart and Lung Transplantation guidelines for the care of cardiac transplant candidates–2006. J. Heart. Lung. Transplant. 25, 1024–1042. doi: 10.1016/j.healun.2006.06.008
Myers, J., Prakash, M., Froelicher, V., Do, D., Partington, S., and Atwood, J. E. (2002). Exercise capacity and mortality among men referred for exercise testing. N. Engl. J. Med. 346, 793–801. doi: 10.1056/NEJMoa011858
Petot, H., Meilland, R., Le Moyec, L., Mille-Hamard, L., and Billat, V. L. (2012). A new incremental test for VO(2)max accurate measurement by increasing VO(2)max plateau duration, allowing the investigation of its limiting factors. Eur. J. Appl. Physiol. 112, 2267–2276. doi: 10.1007/s00421-011-2196-5
Pollock, M. L., Foster, C., Schmidt, D., Hellman, C., Linnerud, A. C., and Ward, A. (1982). Comparative analysis of physiologic responses to three different maximal graded exercise test protocols in healthy women. Am. Heart J. 103, 363–373. doi: 10.1016/0002-8703(82)90275-7
Prieur, F., Busso, T., Castells, J., Bonnefoy, R., Benoit, H., Geyssant, A., et al. (1998). Validity of oxygen uptake measurements during exercise under moderate hyperoxia. Med. Sci.Sports Exerc. 30, 958–962. doi: 10.1097/00005768-199806000-00028
Rietjens, G. J., Kuipers, H., Kester, A. D., and Keizer, H. A. (2001). Validation of a computerized metabolic measurement system (Oxycon-Pro) during low and high intensity exercise. Int. J. Sports Med. 22, 291–294. doi: 10.1055/s-2001-14342
Robergs, R. A., Dwyer, D., and Astorino, T. (2010). Recommendations for improved data processing from expired gas analysis indirect calorimetry. Sports Med. 40, 95–111. doi: 10.2165/11319670-000000000-00000
Salier Eriksson, J., Rosdahl, H., and Schantz, P. (2012). Validity of the Oxycon Mobile metabolic system under field measuring conditions. Eur. J. Appl. Physiol. 112, 345–355. doi: 10.1007/s00421-011-1985-1
Wasserman, K., Hansen, J. E., Sue, D. Y., Whipp, B. J., and Casaburi, R. (1994a). “Appendix D. Calculations, formulas, and examples,” in Principles of Exercise Testing and Interpretation, eds K. Wasserman, J. E. Hansen, D. Y. Sue, B. J. Whipp, and R. Casaburi (Philadelphia, PA: Williams & Wilkins), 454–464.
Wasserman, K., Hansen, J. E., Sue, D. Y., Whipp, B. J., and Casaburi, R. (1994b). “Appendix F. Tables and nomogram,” in Principles of Exercise Testing and Interpretation, eds K. Wasserman, J. E. Hansen, D. Y. Sue, B. J. Whipp, and R. Casaburi (Philadelphia, PA: Williams & Wilkins), 468–470.
Keywords: exercise testing, maximal oxygen consumption, gas exchange, calibration, verification
Citation: Garcia-Tabar I, Eclache JP, Aramendi JF and Gorostiaga EM (2015) Gas analyzer's drift leads to systematic error in maximal oxygen uptake and maximal respiratory exchange ratio determination. Front. Physiol. 6:308. doi: 10.3389/fphys.2015.00308
Received: 06 August 2015; Accepted: 13 October 2015;
Published: 30 October 2015.
Edited by:Kimberly Huey, Drake University, USA
Reviewed by:Kenneth Harrington McKeever, Rutgers, The State University of New Jersey, USA
Jordan A. Guenette, University of British Columbia, Canada
Copyright © 2015 Garcia-Tabar, Eclache, Aramendi and Gorostiaga. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor 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.
*Correspondence: Esteban M. Gorostiaga, firstname.lastname@example.org