Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Med., 25 September 2025

Sec. Intensive Care Medicine and Anesthesiology

Volume 12 - 2025 | https://doi.org/10.3389/fmed.2025.1630871

This article is part of the Research TopicClinical Advances in Acid-base Physiology and Hydroelectrolyte BalanceView all 4 articles

Mathematical correction of the effects of storage time and gas contamination on blood sample measurements

Bahareh Nevirian
Bahareh Nevirian1*Lars Pilegaard ThomsenLars Pilegaard Thomsen1Steen Kre Fagerberg,Steen Kåre Fagerberg2,3Jette Nybo,Jette Nybo2,4Mette Krogh PedersenMette Krogh Pedersen3Kjeld Asbjrn Jensen DamgaardKjeld Asbjørn Jensen Damgaard3Lisha ShastriLisha Shastri1Sren Risom Kristensen,Søren Risom Kristensen2,4Stephen Edward Rees
Stephen Edward Rees1*
  • 1Respiratory and Critical Care (Rcare) Group, Aalborg University, Aalborg, Denmark
  • 2Department of Clinical Medicine, Aalborg University, Aalborg, Denmark
  • 3Anesthesia and Intensive Care Unit, Aalborg University Hospital, Aalborg, Denmark
  • 4Department of Clinical Biochemistry, Aalborg University Hospital, Aalborg, Denmark

Objective: Analysis delay and gas contamination can affect the accuracy of blood measurements. This study uses a mathematical model of blood acid-base chemistry and gas in the sample tubes to calculate values of pH, partial pressures of carbon dioxide (pCO2), partial pressure of oxygen (pO2), oxygen saturation in whole blood (SO2), glucose, and lactate at sample time from measurements with delayed analysis and gas contamination.

Methods: Data were analyzed from two published studies. Study 1: Samples were obtained from 30 critically ill patients in standard blood gas syringes and analyzed after 0, 36, 54, 72, 90, 108, 126, 144, 162, and 180 min. Study 2: Samples were taken from 20 healthy participants in standard blood gas syringes and vacuum tubes (2 mL and 4 mL) and analyzed after 0, 20, and 90 min. Calculated values from the mathematical model were compared to measured values at sample time.

Results: For delays of up to 90 min, the accuracy (mean) and precision (standard deviation (SD)) values calculated at the sample time using syringes and 4-mL vacuum tubes remained within clinically acceptable limits when compared to measured values, with the exception of SO₂ in vacuum tubes. Values represent the mean difference ± standard deviation between calculated and measured values. For syringes, the results were as follows: pH = −0.004 ± 0.011, pCO₂ = 0.08 ± 0.18 kPa, pO₂ = 0.05 ± 0.34 kPa, SO₂ = 0.39 ± 2.21%, glucose = 0.07 ± 0.35 mmol/L, and lactate = 0.13 ± 0.22 mmol/L. For 4-mL vacuum tubes, the results were as follows: pH = 0.006 ± 0.007, pCO₂ = −0.07 ± 0.11 kPa, pO₂ = −0.37 ± 0.34 kPa, SO₂ = −7.79 ± 4.95%, glucose = 0.01 ± 0.11 mmol/L, and lactate = −0.00 ± 0.20 mmol/L. In addition, 2-mL vacuum tubes had poorer accuracy and precision values than syringes and 4-mL vacuum tubes in a subset of cases.

Conclusion: This study has shown that a mathematical model can accurately and precisely calculate blood values at sample time, even following delayed analysis, using both standard blood gas syringes and selected vacuum tubes. This method may have clinical applications in improving the logistics of blood sampling and analysis.

Introduction

Anaerobic metabolism in erythrocytes causes biochemical changes that, over time, affect measured values of acid-base, oxygenation, and metabolism in blood samples. Consequently, the time passed between blood collection and analysis is an important factor influencing the validity of these measurements (19). In addition to the timing of analysis, the type of tube or syringe used for blood collection can influence the validity of measurements. Sampling with vacuum tubes that contain residual air after blood collection results in altered partial pressures of carbon dioxide (pCO₂) and oxygen (pO₂) in the blood due to the diffusion of CO₂ and O₂ between the blood and the residual air in the tube (10).

In previous studies (11, 12), we used a physicochemical mathematical model to simulate changes in acid-base status, oxygenation, electrolytes, and metabolic status under conditions of delayed analysis and gas contamination. In those studies, the model was used in the forward direction—from the sample time to the analysis time—to describe how blood values change with delay and gas exposure. However, the main challenge in clinical practice is not predicting forward, but rather being able to back-calculate the values at the sample time from blood sample measurements affected by delayed analysis and/or gas contamination. In the present study, we therefore applied an inverse use of the model, calculating values at the sample time from the values obtained at the analysis time. This mathematical inversion has not been implemented in the previous study (11, 12) and represents the essential novelty of our approach. Moreover, by calculating the accuracy and precision of these back-calculated values for different delay durations, we provide a practical framework that allows clinicians and researchers to determine the maximum acceptable delay for reliable blood gas analysis. This practical application extends the utility of the physicochemical mathematical model beyond the description of changes toward the correction of delayed and/or gas-contaminated blood samples stored in either standard blood gas syringes or vacuum tubes.

Method

The applied mathematical model describes three compartments that may exist in a blood sample: (1) the gas phase, i.e., a significant gas or air bubble occupying a fraction of the sample tube, as occurs in vacuum tubes; (2) the plasma fraction of blood; and (3) the erythrocyte fraction of blood. Equations were formulated to describe mass balance, mass action, and the physicochemical properties of these compartments, as well as the interfaces between gas and blood and between plasma and erythrocytes, as illustrated in Figure 1. Figure 1 provides a schematic overview of the mathematical model, including key equations and compartmental interactions.

Figure 1
Chemical equilibrium equations for mass balance and gas properties are presented alongside a diagram illustrating gas exchange at the blood-gas interface. The diagram shows the interactions between gas phases (CO₂, O₂), blood phases (Electrolytes, Glucose, Lactate), and erythrocytes, highlighting processes like glycolysis and gas exchanges with reactions involving bicarbonate and hemoglobin.

Figure 1. Equations and physicochemical processes describing the blood and gas phases (see equation numbers in the figure). Mass balance equations (1–18), mass action equations (19–30), gas properties (31–37), physicochemical properties (38–44), and the interface between gas and blood and between the plasma and erythrocyte (45–49).

Detailed explanations of the equations and the corresponding table of abbreviations have been moved to the Supplementary material in accordance with the journal’s guidelines. Readers are encouraged to consult this material to fully understand the methodology before proceeding.

In essence, the model represents the dynamic exchange and equilibration of gases (O₂ and CO₂), ions, and metabolites between the gas, plasma, and erythrocyte compartments of a blood sample. By incorporating diffusion processes, buffering of hydrogen ions, and hemoglobin (Hb) binding reactions, it simulates the alterations in measured values that occur during storage or in the presence of residual gas in the blood tube. In this study, the model is applied in an inverted manner to back-calculate the values at the sample time from those obtained at the analysis time.

Model simulation

The mathematical model in Figure 1 was applied to simulate the plasma values of pH (pHp), partial pressures of carbon dioxide (pCO2p) and oxygen (pO2p), glucose (Glu), lactate (Lac), and oxygen saturation in whole blood (SO2b) at the sample time from measured values at the analysis time. This was performed for the data of two previously conducted studies.

Study 1: Venous blood samples were obtained from 30 critically ill patients admitted to an intensive care unit (ICU), collected in standard blood gas syringes, and analyzed after 0, 36, 54, 72, 90, 108, 126, 144, 162, and 180 min.

Study 2: Peripheral venous blood samples were taken from 20 healthy participants, collected into both standard blood gas syringes and vacuum tubes (2 mL and 4 mL), and analyzed after 0, 20, and 90 min. However, for the simulations presented here and the analysis results shown in the “Results” section, only the vacuum tube data at 20 and 90 min were included. Data from standard blood gas syringes at 90 min in Study 2 were not included in this analysis, as our focus in this study was specifically on data from vacuum tubes. The vacuum tubes used in this study were identical in total volume but designed to draw either 2 mL or 4 mL of blood, leaving a defined volume of the remaining air above the fill line. This remaining air, determined by the manufacturer’s design, was included in the model to reflect typical clinical conditions and potential gas exchange during storage.

In both studies, blood samples were stored at room temperature prior to analysis, and the applied model accounted for the metabolic effects of erythrocyte changes occurring at room temperature.

Study 1 and Study 2 differed. In Study 1, blood samples were taken in standard blood gas syringes (PICO safe heparin syringes, Radiometer, Denmark) with no gas phase, and observable changes in samples were attributed primarily to the metabolic effect of erythrocytes. In Study 2, blood samples were collected in vacuum tubes (VACUETTE® 454001, 454088, Greiner Bio-One, Kremsmünster, Austria), allowing observable changes in the samples to reflect a mix of gas diffusion between the gas phase and blood, as well as erythrocyte metabolism. Accordingly, the simulations performed for each of these two studies differ but use the same computational model illustrated in Figure 1. These simulations are described here and illustrated by the flowchart in Figure 2.

Figure 2
Flowchart comparing Study 1 (standard syringes, no gas phase) and Study 2 (vacuum tubes, gas phase present). Both studies start by measuring values. Study 1 adjusts for anaerobic metabolism and recalculates values, while Study 2 involves blood cooling, gas-blood equilibration, accounting for metabolism, and recalculating values. Steps include model state variable calculations and adjustments for different conditions, like temperature and metabolism.

Figure 2. Model simulation strategy.

In Study 1, the values of pHp, pCO2p, pO2p, SO2b, Hb, Glu, and Lac measured using a blood gas analyzer (ABL 835 FLEX, Radiometer) at analysis time were used to calculate the state variables from the model: blood buffer base (BBb), total concentration of CO₂ (tCO2b) and O₂ (tO2b) in blood, SO2b, and total concentration of Glu and Lac (tGluc and tLac). For Study 1, values of concentrations of O2 and CO2 in the gas compartment were set to zero, reflecting blood sampling in standard blood gas syringes without air contamination. This modeling assumption implies that no significant gas phase was present, and thus, gas exchange between blood and a gas compartment was considered negligible. Consequently, all equations numbered 1, 6, 31–37, 42, and 45–46 related to the gas phase were not applied.

Simulating the conditions at the sample time based on values measured at the analysis time requires adjusting for the acid produced by anaerobic erythrocyte metabolism. To do so, the acid added as part of erythrocyte metabolism due to sampling delay was removed. This removal was performed by reducing the base excess (BE) and buffer base (BB) of whole blood by an acid production rate of 0.008 mmol/L/min, taking into account the duration of the delay, removing an equivalent amount of lactate, and adding half the amount of glucose to account for usual anaerobic metabolism. The value of 0.008 mmol/L/min was the value obtained in Nevirian et al. (11). This modifies the values of the state variables for the mathematical model (Figure 1), which can then be re-solved to calculate all model variables, including pHp, pCO2p, pO2p, SO2b, Glu, and Lac.

In Study 2, calculation from analysis time to sample time requires accounting for both anaerobic metabolism and CO₂ and O₂ distribution between blood and gas compartments. Values measured in blood were reported at 37 °C. However, during the delay before analysis, the blood in the vacuum tube cooled to approximately room temperature. Therefore, the first step was to mathematically convert the measured values to those corresponding to the lower temperature. This conversion was performed at 22 °C using equations from the mathematical model in Figure 1, mainly by adjusting the oxygen dissociation curve for temperature effects (13), while assuming a constant total O₂ concentration in the blood (Equation 28).

Following this cooling, the resulting simulated values of pCO2 and pO2 in blood, i.e., those calculated at 22 °C, were assumed to be in partial equilibrium with those in the gas phase according to equations 45 and 46 of the model. Calculating pCO2 and pO2 in the gas phase following this method allows for calculating the total concentration of CO2 and O2 in the whole system, which includes both the blood and gas phase (tCO₂s and tO₂s), according to equations 1 and 6. Due to mass conservation in this closed tube, these totals remain the same at both sample time and analysis time and can therefore be used—along with two other pieces of information—to enable a complete solution of the mathematical model to calculate all values at sample time. The required information is the composition of the gas phase and the rate of anaerobic metabolism. We assumed that the gas in the vacuum tube at sample time was air, with pO₂ = 21 kPa and pCO₂ = 0.04 kPa, and that the effects of anaerobic metabolism could be calculated, as in the previous study, by modifying blood BE/BB for an acid production rate of 0.008 mmol/L/min.

Some data from Study 2 were excluded from the analysis. After correcting the measured values of pCO₂ and pO₂ at analysis time to 22 °C, occasionally the calculated partial pressure of oxygen in blood remained higher than in air, i.e., >21 kPa. Since all subjects breathed ambient air, pO₂ values above this threshold were not physiologically possible and indicate numerical error rather than true measurement. In these cases, the model attempted to reverse the gas exchange that occurred during storage; however, the result indicated that oxygen remained in the blood at a higher level than that attained under atmospheric breathing conditions. As this issue can be identified within the method, it reflects a limitation in applicability rather than introducing errors. Three data points were excluded due to this, representing 4% of the total calculations performed.

Presentation of results

Results from both Study 1 and Study 2 are presented in two complementary formats: tables and plots.

Tables 1, 2 summarize the mean differences and standard deviations (SD) of the differences for six variables—pH, pCO₂, pO₂, SO₂, Glu, and Lac—for both pre-correction (before applying the model) and post-correction (after applying the model) analyses. For pre-correction, the differences were calculated as the measured baseline values (time 0) minus the measured values from samples analyzed at specific time points: 36, 54, 72, 90, 108, 126, 144, 162, and 180 min for Study 1 and 20 and 90 min for 4-mL and 2-mL vacuum tubes in Study 2. For post-correction, the differences were calculated as the baseline values minus the values back-calculated by the mathematical model at the sample time. The model used measurements from the later analysis times to correct for delayed analysis in Study 1 and both delayed analysis and gas contamination in Study 2, thereby calculating the original blood values at the time the sample was taken.

Table 1
www.frontiersin.org

Table 1. Mean (±SD) of the differences for both pre-correction and post-correction analyses in standard blood gas syringes from Study 1.

Table 2
www.frontiersin.org

Table 2. Mean (±SD) of the differences for both pre-correction and post-correction analyses in 4-mL and 2-mL vacuum tubes from Study 2.

Corresponding plots visualize the post-correction results at these same time points. Additionally, Bland–Altman plots were constructed for each variable across all patients/subjects and are provided in the Supplementary material to evaluate agreement between measured and model-calculated values at sample time (14). The data summarized in Tables 1, 2 are also visualized as error bars showing mean differences and SDs for both pre- and post-correction in the Supplementary material.

Results

As shown in Tables 1, 2, applying the model resulted in a reduction of mean differences (bias) across all time points, except for SO₂ in Study 1, for which the reason will be discussed in the “Discussion” section. The bias between measured values at sample time and measured values at different analysis time points decreased after correction (post-correction), resulting in values closer to the actual measured values at sample time compared to pre-correction. This indicates that the model effectively reduces systematic errors caused by delayed analysis and/or gas contamination in blood samples in both Study 1 and Study 2.

Figures 35 illustrate the mean differences (accuracy) and SD of the differences (precision) obtained by comparing measured and calculated values at the sample time from different analysis time points in Studies 1 and 2. Common to these figures are measurements taken 90 min after sampling, a duration that may be important in relation to blood sample transport to the laboratory. The mean and SD values at 90 min for Study 1 are as follows: pH = −0.004 ± 0.011, pCO₂ = 0.08 ± 0.18 kPa, pO₂ = 0.05 ± 0.34 kPa, SO₂ = 0.39 ± 2.21%, Glu = 0.07 ± 0.35 mmol/L, and Lac = 0.13 ± 0.22 mmol/L. For Study 2, the results at 90 min for 4-mL vacuum tubes are as follows: pH = 0.006 ± 0.007, pCO₂ = −0.07 ± 0.11 kPa, pO₂ = −0.37 ± 0.34 kPa, SO₂ = −7.79 ± 4.95%, Glu = 0.01 ± 0.11 mmol/L, and Lac = −0.00 ± 0.20 mmol/L. For 2-mL vacuum tubes, the results are as follows: pH = 0.021 ± 0.030, pCO₂ = −0.28 ± 0.52 kPa, pO₂ = −1.33 ± 3.96 kPa, SO₂ = −12.64 ± 18.27%, Glu = 0.00 ± 0.14 mmol/L, and Lac = −0.07 ± 0.24 mmol/L.

Figure 3
Six line graphs showing changes in various blood parameters over time: pH, pCO2 (kPa), pO2 (kPa), SO2 (%), Glucose (mmol/L), and Lactate (mmol/L). Each graph plots values at intervals from 0 to 180 minutes, with error bars indicating variability.

Figure 3. Mean differences and standard deviations (SD) of the differences between the measured and model-calculated values at sample times in standard blood gas syringes, calculated from time points of 36, 54, 72, 90, 108, 126, 144, 162, and 180 min in Study 1.

Figure 4
Six line graphs display changes over time for various parameters: pH, pCO2, pO2, SO2%, glucose, and lactate. Time intervals are shown at 0, 20, 50, and 90 minutes. Each graph features mean values with error bars, indicating little variation. The y-axes represent respective units, showing minor fluctuations across all parameters.

Figure 4. Mean differences and standard deviations (SD) of the differences between the measured and model-calculated values at the sample time, calculated from analysis time points of 20 and 90 min in 4-mL vacuum tubes in Study 2.

Figure 5
Six line graphs display changes over time (0, 20, 50, and 90 minutes) for pH, pCO2, pO2, SO2, glucose, and lactate. Each graph shows data points with error bars indicating variability.

Figure 5. Mean differences and standard deviations (SD) of the differences between the measured and model-calculated values at the sample time, calculated from analysis time points of 20 and 90 min in 2-mL vacuum tubes in Study 2.

As illustrated in Figure 3, accuracy changed only slightly over 180 min for all variables except SO₂, while precision remained relatively stable up to 90 min. For Figures 4, 5, accuracy and precision changed only slightly between separate analyses at 20 and 90 min.

Comparison of standard blood gas syringes (Figure 3) and 4-mL vacuum tubes (Figure 4) with 2-mL vacuum tubes (Figure 5) shows that 2-mL tubes exhibit higher mean differences (lower accuracy) and higher SDs (lower precision) in all variables at both 20 and 90 min compared to the other two groups. The reason for this difference is addressed in the “Discussion” section.

Discussion

In this study, we used a mathematical model to simulate acid-base, oxygenation, and metabolism in venous blood samples at the sample time, using measured values of blood samples analyzed at different times of analysis.

The results in standard blood gas syringes and 4-mL vacuum tubes illustrate that, for delays in the analysis up to 90 min, values describing acid-base, oxygenation, and metabolism at sample time can be accurately and precisely calculated. The differences between model-calculated and measured values of pH and pCO₂ at 90 min remain within the range previously reported for repeatability, where consecutive blood gas measurements from the same patient showed standard deviations of 0.01 for pH, 0.16 kPa for pCO₂, 0.15 kPa for pO₂, and 1.49% for SO₂ (15). Although the SD of differences between model-calculated and measured pO₂ values in standard blood gas syringes and 4-mL vacuum tubes was higher than Mallat’s SD (15), this is likely due to oxygen degradation occurring during storage. Our model substantially corrects the bias and variability introduced by storage, bringing the values closer to the true values measured at the sample time, as shown in Supplementary Figures 7, 10. However, the results from 4-mL vacuum tubes in Study 2 are better than those from standard blood gas syringes in Study 1, which can be attributed to greater Hb variability observed in Study 1. Thus, the residual variability reflects physiological and technical limits rather than model inadequacy, supporting the model’s validity for calculating pO₂ in stored samples from both standard blood gas syringes and 4-mL vacuum tubes. For vacuum tubes, pO₂ could be calculated accurately and precisely. However, because the pO₂ values were low (2–6 kPa), even small calibration errors in the oxygen dissociation curve led to large errors in the calculated SO₂, resulting in relatively poor accuracy for SO₂. This would not be expected at higher oxygenation levels. For SO₂ in Study 1, the bias did not decrease substantially, which could be attributed to variability in Hb present in this study. To accurately model the oxygen saturation curve, the model assumes stable Hb levels. Otherwise, shifts in the oxygen dissociation curve can make SO₂ estimation unreliable—even if pO₂ was well-modeled.

In addition, errors in pH and pCO2 at 90 min for standard blood gas syringes and 4-mL vacuum tubes were within those considered useful for clinical interpretation in clinical emergencies, i.e., 0.05 (95% CI 0.04–0.06) for pH and 0.88 (95% CI: 0.75, 1.01 kPa) for pCO2 reported by Rang et al. (16). Errors in calculated glucose and lactate were always within 0.5 mmol/L, even in the 2-mL vacuum tube data. These values were adequate to identify glucose abnormalities or hyperlactatemia.

The method presented in this study may therefore be useful in clinical practice and change current opinion as to the use of blood with analysis delay for acid-base status and oxygenation (19). Using blood following a delay of 90 min might improve the blood sampling logistics in both out-clinic and in-hospital situations. Samples could be taken at the start of an ambulance transport prior to the start of therapy, analyzed in the hospital, and compared to a newly drawn blood sample in order to evaluate prehospital treatment. Similarly, samples could be taken in the ward and analyzed centrally, reducing the need for decentralized point-of-care technology. Such an approach might be integrated into the clinical workflow, with the only additional data required being the delay duration and the type of tube used. In addition, the use of vacuum tubes may reduce the number of venous punctures, as samples could be taken together with those tubes for other purposes.

The same mathematical model has been applied here for the analysis of standard blood gas syringes, vacuum tubes, different delays, and both healthy subjects and critically ill patients. The model seems to be robust over a wide range of applications and patients. This, however, was not the case for calculations made in vacuum tubes with only 2 mL of blood. For vacuum tubes with a small blood volume—that is, a large gas volume—errors in calculating pO₂ and pCO₂ in blood during cooling result in large errors in the calculation of the mass of O₂ and CO₂ in the gas phase. This finding indicates that the method presented here is sensitive to vacuum tubes with large gas volumes. However, as shown in Supplementary Figure 9, these errors appear only in a small subset of subjects and are associated with poor Bland–Altman agreement, largely due to a few outliers. This reflects a numerical limitation of the current modeling approach, rather than a physiological inconsistency, and occurs only when corrected pO₂ values exceed ambient air levels. Further studies are required to understand the conditions resulting in these outliers, the identification of which may allow the use of the method in vacuum tubes with large gas volumes. However, it is encouraging that 4-mL vacuum tubes are the most widely used in clinical practice. In addition, the study is limited to only normal subjects in vacuum tubes and severely ill patients in standard blood gas syringes and limited to specific brands of syringes and vacuum tubes. Further analysis will be needed to assess the generalizability and limitations of the method.

The estimation of acid production (0.008 mmol/L/min) was derived from the dataset in Study 1, as published previously (11), and verified in the prior publication of Study 2 (12). However, this value can change under varying physiological conditions, including different hematocrit, glucose levels, oxygenation, temperature, and other patient-specific factors. While the fixed rate applied here provided good agreement across both critically ill patients in Study 1 and healthy subjects in Study 2, further research is required to assess model outcomes under broader physiological conditions.

We acknowledge that the cited values from Mallat and Rang (15, 16) reflect practical clinical thresholds for interpretation rather than statistical measures of variability, such as the smallest detectable difference (SDD). While these values provide useful context for clinical relevance, statistical benchmarks such as the SDD may offer additional value in future methodological assessments focused on test–retest reliability.

The applied mathematical model uses the exact time of sample analysis relative to collection. However, in routine clinical practice, sampling times were not always recorded with minute-level precision and may be approximate. Further studies are required to assess how this uncertainty could influence model outcomes.

Conclusion

This study evaluated a method for calculating values at the sample time following delayed analysis in standard blood gas syringes and vacuum tubes. The results demonstrate that accurate and precise values can be obtained for syringes and certain vacuum tubes. This method may have clinical applications in improving the logistics of blood sampling and analysis and increasing the usability of blood samples acquired in out-of-hospital or prehospital settings, where analysis of drawn blood may take several hours.

Data availability statement

The data analyzed in this study is subject to the following licenses/restrictions: the dataset consists of anonymized patient data and is not publicly available due to ethical and privacy restrictions. Access to the data may be granted upon reasonable request and with appropriate institutional and ethical approvals. Requests to access these datasets should be directed to BN email: YmFoYXJlaG5AaHN0LmFhdS5kaw==.

Ethics statement

This study used two previously approved datasets: (1) blood samples from critically ill patients approved as a quality control evaluation by the local hospital administration at Region Hospital, North Jutland, Denmark and (2) blood samples from healthy volunteers approved by the Ethics Committee for the Region of North Jutland (Reference: N-20230009). 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

BN: Writing – review & editing, Validation, Investigation, Formal analysis, Software, Resources, Data curation, Writing – original draft, Visualization, Conceptualization. LT: Writing – review & editing. SF: Investigation, Writing – review & editing, Resources. JN: Resources, Writing – review & editing. MP: Resources, Investigation, Writing – review & editing. KD: Investigation, Writing – review & editing, Resources. LS: Writing – review & editing. SK: Writing – review & editing, Resources. SR: Conceptualization, Funding acquisition, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded by Roche Medical, who played no role in the design of the study or the evaluation of the data.

Conflict of interest

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.

Generative AI statement

The authors declare that Gen AI was used in the creation of this manuscript. Generative AI was used to assist with proofreading, grammar checking, and language polishing.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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/fmed.2025.1630871/full#supplementary-material

References

1. Zavorsky, GS, Van Wijk, XMR, Gasparyan, S, Stollenwerk, NS, and Brooks, RA. Stability of whole blood electrolyte specimens at room temperature vs. slushed ice conditions. J Appl Lab Med. (2022) 7:541–54. doi: 10.1093/jalm/jfab089

PubMed Abstract | Crossref Full Text | Google Scholar

2. Zavorsky, GS, Gasparyan, S, Stollenwerk, NS, and Brooks, RA. Stability of whole blood lactate specimens at room temperature versus slushed ice conditions. Respir Care. (2021) 66:494–500. doi: 10.4187/respcare.08023

PubMed Abstract | Crossref Full Text | Google Scholar

3. Zavorsky, GS, and Van Wijk, XMR. The stability of blood gases and CO-oximetry under slushed ice and room temperature conditions. Clin Chem Lab Med. (2023) 61:1750–9. doi: 10.1515/cclm-2022-1085

PubMed Abstract | Crossref Full Text | Google Scholar

4. Dent, RG, Boniface, DR, Fyffe, J, and Yousef, Z. The effects of time delay and temperature on capillary blood gas measurements. Respir Med. (1999) 93:794–7. doi: 10.1016/S0954-6111(99)90264-0

PubMed Abstract | Crossref Full Text | Google Scholar

5. Knowles, TP, Mullin, RA, Hunter, JA, and Douce, FH. Effects of syringe material, sample storage time, and temperature on blood gases and oxygen saturation in arterialized human blood samples. Respir Care. (2006) 51:732–6.

Google Scholar

6. Liss, HP, and Piyne, CP. Stability of blood gases in ice and at room temperature. Chest. (1993) 103:1120–2. doi: 10.1378/chest.103.4.1120

PubMed Abstract | Crossref Full Text | Google Scholar

7. Srisan, P, Udomsri, T, Jetanachai, P, Lochindarat, S, and Kanjanapattanakul, W. Effects of temperature and time delay on arterial blood gas and electrolyte measurements. J Med Assoc Thail. (2011) 94:S9–S14.

Google Scholar

8. Woolley, A, and Hickling, K. Errors in measuring blood gases in the intensive care unit: effect of delay in estimation. J Crit Care. (2003) 18:31–7. doi: 10.1053/jcrc.2003.YJCRC7

PubMed Abstract | Crossref Full Text | Google Scholar

9. Biswas, CK, Ramos, JM, Agroyannis, B, and Kerr, DNS. Blood gas analysis: effect of air bubbles in syringe and delay in estimation. Br Med J. (1982) 284:923–7. doi: 10.1136/bmj.284.6320.923

PubMed Abstract | Crossref Full Text | Google Scholar

10. Tang, NY, Leanse, JH, Tesic, V, and van Wijk, XMR. A comparison of venous blood gas analysis with paired specimens collected in syringes and evacuated blood collection tubes. Clin Chim Acta. (2020) 510:671–4. doi: 10.1016/j.cca.2020.08.023

PubMed Abstract | Crossref Full Text | Google Scholar

11. Nevirian, B, Fagerberg, SK, Pedersen, MK, Kristensen, SR, Damgaard, KAJ, Rees, SE, et al. A physio-chemical mathematical model of the effects of blood analysis delay on acid-base, metabolite and electrolyte status: evaluation in blood from critical care patients. Clin Chem Lab Med. (2024) 63:1139–52. doi: 10.1515/cclm-2024-1350

PubMed Abstract | Crossref Full Text | Google Scholar

12. Nevirian, B, Thomsen, LP, Nybo, J, Shastri, L, Kristensen, SR, and Rees, SE. A physio-chemical mathematical model of the effects of O2 and CO2 diffusion between the blood and gas phases of venous blood sampled in vacuum tubes. Clin Chim Acta. (2025) 572:120270. doi: 10.1016/j.cca.2025.120270

PubMed Abstract | Crossref Full Text | Google Scholar

13. Schmidt, C, and Müller-Plathe, O. Stability of pO2, pCO2 and pH in heparinized whole blood samples: influence of storage temperature with regard to leukocyte count and syringe material. Eur J Clin Chem Clin Biochem. (1992) 30:767–73. doi: 10.1515/cclm.1992.30.11.767

PubMed Abstract | Crossref Full Text | Google Scholar

14. Bland, JM, and Altman, DG. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. (1986) 1:307–10.

PubMed Abstract | Google Scholar

15. Mallat, J, Lazkani, A, Lemyze, M, Pepy, F, Meddour, M, Gasan, G, et al. Repeatability of blood gas parameters, PCO2 gap, and PCO2 gap to arterial-to-venous oxygen content difference in critically ill adult patients. Medicine. (2015) 94:e415. doi: 10.1097/MD.0000000000000415

PubMed Abstract | Crossref Full Text | Google Scholar

16. Rang, LCF, Murray, HE, Wells, GA, and MacGougan, CK. Can peripheral venous blood gases replace arterial blood gases in emergency department patients? CJEM. (2002) 4:7–15. doi: 10.1017/s1481803500006011

Crossref Full Text | Google Scholar

17. Rees, SE, and Andreassen, S. Mathematical models of oxygen and carbon dioxide storage and transport: the Acid-Base chemistry of blood. Crit Rev Biomed Eng. (2005) 33:209–64. doi: 10.1615/CritRevBiomedEng.v33.i3.10

PubMed Abstract | Crossref Full Text | Google Scholar

18. Rees, SE, Klaestrup, E, Handy, J, Andreassen, S, and Kristensen, SR. Mathematical modelling of the acid-base chemistry and oxygenation of blood: a mass balance, mass action approach including plasma and red blood cells. Eur J Appl Physiol. (2010) 108:483–94. doi: 10.1007/s00421-009-1244-x

PubMed Abstract | Crossref Full Text | Google Scholar

19. Funder, J, and Wieth, JO. Chloride and hydrogen ion distribution between human red cells and plasma. Acta Physiol Scand. (1966) 68:234–45. doi: 10.1111/j.1748-1716.1966.tb03423.x

Crossref Full Text | Google Scholar

20. Holtkamp, HC, Verhoef, NJ, and Leijnse, B. The difference between the glucose concentrations in plasma and whole blood. Clin Chim Acta. (1975) 59:41–9. doi: 10.1016/0009-8981(75)90216-8

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: acid-base, mathematical modeling, analysis delay, gas contamination, venous blood sampling

Citation: Nevirian B, Thomsen LP, Fagerberg SK, Nybo J, Pedersen MK, Damgaard KAJ, Shastri L, Kristensen SR and Rees SE (2025) Mathematical correction of the effects of storage time and gas contamination on blood sample measurements. Front. Med. 12:1630871. doi: 10.3389/fmed.2025.1630871

Received: 18 May 2025; Accepted: 05 September 2025;
Published: 25 September 2025.

Edited by:

Zhongheng Zhang, Sir Run Run Shaw Hospital, China

Reviewed by:

Francesco Zadek, University of Milano Bicocca, Italy
Claudio Pellegrino, Agostino Gemelli University Polyclinic (IRCCS), Italy
Martin Krbec, Charles University, Czechia
Micah Liam Arthur Heldeweg, Amsterdam University Medical Center, Netherlands

Copyright © 2025 Nevirian, Thomsen, Fagerberg, Nybo, Pedersen, Damgaard, Shastri, Kristensen and Rees. 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) 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.

*Correspondence: Bahareh Nevirian, YmFoYXJlaG5AaHN0LmFhdS5kaw==; Stephen Edward Rees, c3JAaHN0LmFhdS5kaw==

Disclaimer: 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.