Science review: Quantitative acid–base physiology using the Stewart model
Critical Care volume 8, Article number: 448 (2004)
There has been renewed interest in quantifying acid–base disorders in the intensive care unit. One of the methods that has become increasingly used to calculate acid–base balance is the Stewart model. This model is briefly discussed in terms of its origin, its relationship to other methods such as the base excess approach, and the information it provides for the assessment and treatment of acid–base disorders in critically ill patients.
Acid–base derangements are commonly encountered in the critical care unit , and there is renewed interest in the precise description of these disorders in critically ill patients [2–5]. This new interest has led to a renovation of the quantitative assessment of physiological acid–base balance, with increasing use of the Stewart model (strong ion difference [SID] theory) to calculate acid–base balance in the critically ill [2, 3, 6, 7]. This method is discussed, particularly as it pertains to the metabolic component of acid–base derangements, as one of several approaches that may be used in the intensive care unit for quantitative evaluation. As with any mathematical model, a basic understanding of its principles is useful for proper application and interpretation.
All equilibrium models of acid–base balance utilize the same basic concept. Under the assumption of equilibrium or a steady-state approximation to equilibrium, some property of the system (e.g. proton number, proton binding sites, or charge, among other possible properties) is enumerated from the distribution of that property over the various species comprising the system, according to the energetics of the system manifested through the relevant equilibrium constants of the various species under a given set of conditions [5, 8–12]. This function is calculated at the normal values and then the abnormal values; from these the degree of change is obtained to give information about the clinical acid–base status of the patient. All of the apparently 'different' methods for assessing acid–base balance arise from this common framework [5, 12].
In the Stewart method, charge is taken as the property of interest [7, 11, 13]. Using this property, acid–base status may be expressed for a single physiologic compartment, such as separated plasma, as follows [7, 10, 11, 13]:
Strong ions are those that do not participate in proton transfer reactions, and the SID is defined as the difference between the sum of positive charge concentrations and the sum of negative charge concentrations for those ions that do not participate in proton transfer reactions. Cn are the analytical concentrations of the various buffer species also in the compartment (e.g. of the buffer amino acid groups on albumin), and are the average charges of those various species. The can be expressed as functions of pH and equilibrium constants [11, 12], and it is therefore convenient to calculate SID using Eqn 1 from the pH and the concentrations of relatively few buffer species, as opposed to a direct calculation from a measurement of all of the various strong ion species. In many implementations of the Stewart method, contributions from the water equilibrium and from carbonate species other than bicarbonate are neglected, because these are small under physiologic conditions [11, 14, 15]. The first term in Eqn 1 may then be equated with the bicarbonate concentration, with the remaining terms referring to other buffer species [11, 14].
Plasma physiologic pH is then determined by the simultaneous solution of Eqn 1 and the Henderson-Hasselbalch Equation:
Where for human plasma pK' = 6.103. S = 0.0306 is the equilibrium constant between aqueous and gas phase CO2[16, 17]. [HCO3-] is the concentration of plasma bicarbonate in mmol/l, and PCO2 is the partial CO2 tension in Torr.
SID = [HCO3-] (3)
In this approach the metabolic component of an acid–base disorder is quantified as the change in plasma bicarbonate concentration (Δ[HCO3-]) , which by Eqn 3 is also equal to ΔSID. This method is often sufficient and has been used successfully to diagnose and treat countless patients, but it has also been criticized as not strictly quantitative [19, 20]. [HCO3-] depends upon the PCO2 and does not provide complete enumeration of all species, because albumin and phosphate also participate in plasma acid–base reactions [15, 17, 20, 21].
A more complete calculation may be undertaken for better approximation by including more terms in the series in Eqn 1. In addition, although
is a nonlinear function of pH, it can be approximated over the physiologic range by a more computationally convenient linear form, such that for plasma the following explicit expression is obtained [11, 12, 15]:
SID = [HCO3-] + CAlb (8.0pH - 41) + CPhos (0.30pH - 0.4) (4)
Where CAlb and CPhos are plasma albumin and phosphate concentrations, respectively. All concentrations are in mmol/l. One may multiply albumin in g/dl by 0.15 to obtain albumin in mmol/l, and phosphate in mg/dl by 0.322 to get phosphate in mmol/l. The factors 8.0 and 0.30 are the molar buffer values of albumin and phosphate, respectively. The buffer value is the change in of a species for a one unit change in pH [5, 11, 17]. Note that the ability of a system to resist pH change also increases with CAlb and CPhos.
Equation 4 was obtained via a term by term summation over all of the buffer groups in albumin and of phosphoric acid, as performed by Figge and coworkers [15, 21]. The theoretical basis for the validity of this approach is well established , and Eqn 4 has been shown to reproduce experimental data well [11, 12, 15, 21, 22]. Some authors have argued that the effects of plasma globulins should also be considered for better approximation [17, 20, 23, 24], although other calculations suggest that the consideration of globulins would be of little clinical significance in humans .
Consideration of the change in SID using Eqn 4 between normal and abnormal states at constant albumin and phosphate concentrations gives the following:
ΔSID = Δ[HCO3-] + (8.0CAlb + 0.30CPhos)ΔpH (5)
Which is recognized to be of the same form and numerically equivalent to the familiar Van Slyke equation for plasma, yielding the plasma base excess (BE) [5, 11, 17, 25]. Furthermore, Eqn 4 is of the same form as the CO2 equilibration curve of the BE theory presented by Siggaard-Andersen [11, 17, 20, 25]. The BE approach and the Stewart method are equivalent at the same level of approximation [11, 12, 26].
Strong ion gap
A widely used concept arising from the Stewart approach is the strong ion gap (SIG), which was popularized by Kellum  and Constable . This relies upon a direct calculation of the SID as, for example, the following:
Where SIDm is the measured SID . This direct measurement is then compared with that generated via Eqn 4:
SIG = SIDm - SID (7)
This gives a higher level version of the familiar plasma anion gap [1, 18]. Some publications have used the notation SIDa (for SID apparent) to refer to the variable SIDm calculated using Eq. 6, and SIDe (SID effective) to refer to that calculated using Eqn 4 [2, 3, 15, 27]. SIG has been shown to predict the presence of unmeasured ions better than the conventional anion gap , as might be expected, given that more variables are taken into account. Some unmeasured ions that are expected to contribute to the SIG are β-hydroxybutyrate, acetoacetate, sulfates, and anions associated with uremia .
Changes in noncarbonate buffer concentration
ΔSID expressed through the relationship of Eqn 5 unambiguously quantifies the nonrespiratory component of an acid–base disturbance in separated plasma [11, 17], with the total concentrations of amphoteric species such as albumin and phosphate remaining constant [11, 12, 17]. An amphoteric substance is one that can act as both an acid and a base. Stewart and other investigators [4, 7, 29–33], though, have emphasized the role played by changes in the noncarbonate buffer concentrations in acid–base disorders. When the noncarbonate buffer concentrations change, the situation becomes more complex, and in general a single parameter such as ΔSID no longer necessarily quantifies the metabolic component of an acid–base disorder, and enough variables must be examined to characterize the disorder unambiguously. Examples below demonstrate this point when the concentrations of noncarbonate buffers change, through a pathologic process or through resuscitation.
Table 1 gives several examples for separated human plasma, including the normal values of case 1. Case 2 demonstrates a metabolic acidosis with constant noncarbonate buffer concentrations, in which the ΔSID of -10 mmol/l quantifies the metabolic component of the acid–base disorder , which has been described as a strong ion acidosis . Case 3 gives values for the fairly common occurrence of isolated hypoproteinemia. This too gives a ΔSID of -10 mmol/l, although the total weak acid and weak base concentrations have both decreased . The physiological interpretation of this condition in terms of acid–base pathology is the subject of debate [3, 6, 12, 20, 31, 34]. Considering this to be an acid–base disorder, some authors would classify this case as hypoproteinemic alkalosis with a compensating SID acidosis [4, 6, 30–32]. More generally, this has been termed a buffer ion alkalosis with compensating strong ion acidosis . If the mechanism of hypoalbuminemia is en bloc loss of charged albumin with counterions in tow, for example in nephrotic syndrome, then it seems dubious to describe this process as compensation in the usual physiologic sense. Also, note that both cases 2 and 3 have the same decrease in SID, but the individual in case 2 is expected to be quite sick with acidemia whereas the patient in case 3 is probably not acutely ill, except for the effects of low oncotic pressure.
Although it has been suggested that alkalosis can result from hypoproteinemia, with patients without adequate compensation becoming alkalemic [29, 32], the idea of alterations in protein concentration as acid–base disorders per se has been questioned [3, 20]. The concept of the normal SID changing as a function of protein concentration has been suggested [3, 11, 12]. In such an instance, ΔSID again quantifies the metabolic component of an acid–base disturbance, essentially renormalizing the noncarbonate buffer concentrations to the abnormal values [11, 12]. This is basically what has been advocated in the past for BE [20, 34], in which Eqn 5 uses the abnormal protein and phosphate concentrations for CAlb and CPhos. Thus, the SID of 29 mmol/l in case 3 is said to be normal for the decreased albumin concentration , giving a ΔSID of 0 mmol/l. This individual will, however, be more susceptible to acidemia or alkalemia for a given derangement, as expressed through the molar buffer values and noncarbonate buffer concentrations, than would a normal individual . If SID is not renormalized as described above, then BE and ΔSID differ by an added constant [11, 12].
Another interesting issue is raised in the treatment of patients with intravenous albumin or other amphoteric species. Kellum previously pointed out that, based on the SID, one might think that albumin solutions with a SID of 40–50 mmol/l would be alkalinizing to the blood, even though their pH is close to 6.0 . This apparent paradox is resolved by again realizing that, for amphoteric substances, one is not only changing the SID but also increasing both the total weak acid and weak base concentrations by increasing the total protein concentration [9, 11]. This highlights the point made by Stewart concerning the necessity of considering all variables in assessing acid–base balance [7, 13]. A complete calculation yields what is intuitively predicted – that such a solution is in fact acidifying to blood (unpublished data). One might further speculate that the administration of 'unbuffered' albumin to patients may contribute to the reason why this treatment has not been more successful in the critically ill . Extensive quantitative discussions regarding the acid–base balance of administered fluids have typically not been given in publications on resuscitation with amphoteric colloids [36–39], although this is an issue that should be examined. Constable  recently gave a brief quantitative discussion of acid–base effects of giving various crystalloids.
Model for whole blood
Several points arise in the comparison of SID with BE, as has been performed in a number of studies [33, 38, 41–44]. This is in some respects a misplaced comparison, because BE represents a difference whereas SID does not [11, 26]. The corresponding variable to SID in the BE formalism is the concentration of total proton binding sites, while the BE represents the change in this quantity from the normal value, and corresponds to ΔSID [11, 12, 17, 26]. More significant, clinical studies using Stewart theory have calculated the separated plasma SID, while making comparison with the BE for whole blood or the standard base excess (SBE) [33, 38, 41, 42], rather than the corresponding plasma BE. Furthermore, consideration of only the plasma compartment creates a potential source of error, because separated plasma versions of the Stewart method quantify only a portion of the acid–base disorder [12, 17, 45]. An equation for the SID of whole blood has recently been derived, partly to address this issue .
Where φ(E) is the hematocrit, CHgb(B) is the hemoglobin concentration of whole blood, and CDPG(E) is the 2, 3-diphosphoglycerate concentration in the erythrocyte. Again, concentrations are in mmol/l, and one may multiply hemoglobin in g/dl by 0.155 to obtain hemoglobin in mmol/l. The normal 2, 3-diphosphoglycerate concentration in the erythrocyte is 6.0 mmol/l . The 'P', 'B', and 'E' designations stand for plasma, whole blood, and erythrocyte fluid, respectively. The corresponding Van Slyke form has also been obtained, and is numerically identical to BE for whole blood .
The SBE, as mentioned above, is also widely used [3, 17, 20, 25]. This parameter reflects the extracellular acid–base status and approximates the in vivo BE for the organism [17, 20, 25]. The Van Slyke equation for SBE approximates this situation via a 2:1 dilution of whole blood in its own plasma [17, 20, 25]. It should be borne in mind, therefore, that Eqn 4 may prove more concordant with clinical data than Eqn 8, since the plasma expression may produce values closer to the in vivo condition because of the distribution functions of various species across the whole organism .
Stewart theory and mechanism
Finally, the Stewart model is taken by some to be a mechanistic description of acid–base chemistry in which changes only occur by alteration in PCO2, SID, or noncarbonate buffer concentrations because these are the only true independent variables; changes never occur by addition or removal of H+to the system or by changes in [HCO3-] because these are dependent variables [7, 13]. It is said that because the Stewart theory provides mechanistic information, it is superior to the BE approach [3, 35, 46, 47]. Support for this point of view is offered in the form of philosophic arguments regarding the nature of independence [7, 13], as well as studies showing that the Stewart model accurately predicts what is observed experimentally [30, 42, 44, 48]. However, like the BE approach and like any other method derived from considerations involving the calculation of interval change via the assessment of initial and final equilibrium states, the Stewart method does not produce mechanistic information [8, 35]. These are basically bookkeeping methods. To believe otherwise risks falling prey to the computo, ergo est (I calculate it, therefore it is) fallacy. What is thus required for mechanistic understanding is the collection of actual mechanistic data, perhaps obtainable through isotopic labeling and kinetics experiments.
Both experimental and theoretical data have shown that the Stewart method is accurate for describing physiological acid–base status, and the use of the SIG potentially offers an improvement over the traditional anion gap, but because the Stewart method proceeds from the same common framework as the BE approach, it theoretically offers no quantitative advantage over BE at corresponding levels of approximation [11, 12, 26, 35, 49]. As such, it remains to be seen whether the renovation of acid–base assessment afforded by the Stewart approach constitutes a radical new architecture for understanding acid–base physiology, or whether it is simply a new façade.
partial CO2 tension
standard base excess
strong ion difference
strong ion gap.
Irwin RS, Rippe JM, editors: Intensive Care Medicine. 2003, Philadelphia, PA: Lipincott, Williams & Wilkins
Kellum JA: Metabolic acidosis in the critically ill: lessons from physical chemistry. Kidney Int. 1998, Suppl 66: S81-S86.
Kellum JA: Determinants of blood pH in health and disease. Crit Care. 2000, 4: 6-14. 10.1186/cc644.
Constable PD: Clinical assessment of acid-base status: comparison of the Henderson-Hasselbalch and strong ion approaches. Vet Clin Path. 2000, 29: 115-128.
Corey HE: Stewart and beyond: new models of acid–base balance. Kidney Int. 2003, 64: 777-787. 10.1046/j.1523-1755.2003.00177.x.
Constable PD: Clinical assessment of acid-base status: strong ion difference theory. Vet Clin N Am. 1999, 15: 447-471.
Stewart PA: Modern quantitative acid-base chemistry. Can J Physiol Pharmacol. 1983, 61: 1441-1461.
Tanford C: Physical Chemistry of Macromolecules. 1961, New York: John Wiley & Sons, Inc
Guenther WB: Unified Equilibrium Calculations. 1991, New York: John Wiley & Sons, Inc
Butler JN, Cogley DR: Ionic Equilibrium: Solubility and pH Calculations. 1998, New York: John Wiley & Sons, Inc
Wooten EW: Analytic calculation of physiological acid-base parameters in plasma. J Appl Physiol. 1999, 86: 326-334.
Wooten EW: Calculation of physiological acid-base parameters in multicompartment systems with application to human blood. J Appl Physiol. 2003, 95: 2333-2344.
Stewart PA: Independent and dependent variables of acid-base control. Respir Physiol. 1978, 33: 9-26. 10.1016/0034-5687(78)90079-8.
Constable PD: A simplified strong ion model for acid-base equilibria: application to horse plasma. J Appl Physiol. 1997, 83: 297-311.
Figge J, Mydosh T, Fencl V: Serum proteins and acid-base equilibria: a follow-up. J Lab Clin Med. 1992, 120: 713-719.
Burus CA, Ashwood ER, Editors: Tietz Textbook of Clinical Chemistry. 1994, Philadelphia, PA: Saunders, 2
Siggaard-Andersen O: The Acid–base Status of the Blood. 1974, Baltimore, MD: Williams and Wilkins, 4
Nairns RG, Emmett M: Simple and mixed acid-base disorders: a practical approach. Medicine (Baltimore). 1980, 59: 161-187.
Severinghaus JW: Siggaard-Andersen and the 'Great TransAtlantic Acid-Base Debate'. Scand J Lab Invest. 1993, Suppl 214: 99-104.
Siggaard-Andersen O, Fogh-Andersen N: Base excess or buffer base (strong ion difference) as a measure of a non-respiratory acid–base disturbance. Acta Anaesthesiol Scand. 1995, Suppl 107: 123-128.
Figge J, Rossing TH, Fencl V: The role of serum proteins in acid–base equilibria. J Lab Clin Med. 1991, 117: 453-467.
Watson PD: Modeling the effects of proteins on pH in plasma. J Appl Physiol. 1999, 86: 1421-1427.
Siggaard-Andersen O, Rorth M, Strickland DAP: The buffer value of plasma, erythrocyte fluid and whole blood. In Blood pH, Gases and Electrolytes. Workshop on pH and Blood Gases, National Bureau of Standards, 1975. National Bureau of Standards (US), Special Publications. 1977, 450: 11-19.
Staempfli HR, Constable PD: Experimental determination of net protein charge and Atot and Ka of nonvolatile buffers in human plasma. J Appl Physiol. 2003, 95: 620-630.
Siggaard-Andersen O: The Van Slyke equation. Scand J Clin Lab Invest. 1977, Suppl 146: 15-20.
Schlichtig R: [Base excess] vs. [strong ion difference]. Which is more helpful?. Adv Exp Med Biol. 1997, 411: 91-95.
Kellum JA, Kramer DJ, Pinsky MR: Strong ion gap: a methodology for exploring unexplained anions. J Crit Care. 1995, 10: 51-55. 10.1016/0883-9441(95)90016-0.
Constable PD, Hinchcliff KW, Muir WW: Comparison of anion gap and strong ion gap as predictors of unmeasured strong ion concentration in plasma and serum from horses. Am J Vet Res. 1998, 59: 881-887.
McAuliffe JJ, Lind LJ, Lieth DE, Fencl V: Hypoproteinemic alkalosis. Am J Med. 1986, 81: 86-90. 10.1016/0002-9343(86)90187-7.
Rossing TH, Maffeo N, Fencl V: Acid–base effects of altering plasma protein concentration in human blood in vitro. J Appl Physiol. 1986, 61: 2260-2265.
Jabor A, Kazda A: Modelling of acid–base equilibria. Acta Anaesthesiol Scand. 1995, Suppl 107: 119-122.
Wilkes P: Hypoproteinemia, strong ion difference, and acid-base status in critically ill patients. J Appl Physiol. 1998, 84: 1740-1748.
Rocktaeschel J, Morimatsu H, Uchino S, Goldsmith D, Poustie S, Story D, Gutteridge G, Bellomo R: Acid–base status of critically ill patients with acute renal failure: analysis based on Stewart-Figge methodology. Crit Care. 2003, 7: R60-R66. 10.1186/cc2333.
Moon JB: Abnormal base excess curves. Pediatr Res. 1967, 1: 333-340.
Wooten EW: Strong ion difference theory: more lessons from physical chemistry. Kidney Int. 1998, 54: 1769-1770. 10.1046/j.1523-1755.1998.00170.x.
Cochrane Injuries Group Reviewers: Human albumin administration in critically ill patients: systematic review of randomized controlled trials. BMJ. 1998, 317: 235-240.
Hayhoe M, Bellomo R, Liu G, McNicol L, Buxton B: The aetiology and pathogenesis of cardiopulmonary bypass-associated metabolic acidosis using polygeline pump prime. Int Care Med. 1999, 25: 680-685. 10.1007/s001340050930.
Rehm M, Orth V, Scheingraber S, Kreimeier U, Brechtelsbauer H, Finsterer U: Acid–base changes caused by 5% albumin versus 6% hydroxyethyl starch solution in patients undergoing acute normovolemic hemodilution. Anesthesiol. 2000, 93: 1174-1183. 10.1097/00000542-200011000-00007.
Rehm M, Finsterer U: Treating intraoperative hyperchloremic acidosis with sodium bicarbonate or tris-hydroxylmethyl aminomethane: a randomized prospective study. Anesth Analg. 2003, 96: 1201-1208. 10.1213/01.ANE.0000048824.85279.41.
Constable P: Stewart approach is not always a practical tool. Anesth Analg. 2004, 98: 271-272.
Kellum JA: Fluid resuscitation and hyperchloremic acidosis in experimental sepsis: improved short-term survival and acid-base balance with Hextend compared with saline. Crit Care Med. 2002, 30: 300-305. 10.1097/00003246-200202000-00006.
Morgan TJ, Venkatesh B, Hall J: Crystalloid strong ion difference determines metabolic acid-base change during in vitro hemodilution. Crit Care Med. 2002, 30: 157-160. 10.1097/00003246-200201000-00022.
Scheingraber S, Rehm M, Sehmisch C, Finsterer U: Rapid saline infusion produces hyperchloremic acidosis in patients undergoing gynecologic surgery. Anesthesiology. 1999, 90: 1265-1270. 10.1097/00000542-199905000-00007.
Waters JH, Bernstein CA: Dilutional acidosis following hetastarch or albumin in healthy volunteers. Anesthesiol. 2000, 93: 1184-1187. 10.1097/00000542-200011000-00008.
Davenport HW: The ABC of Acid-base Chemistry. 1974, Chicago, IL: University of Chicago Press, 6
Waters JH, Scanlon TS, Howard RS, Leivers D: Role of minor electrolytes when applied to Stewart's acid–base approach in an acidotic rat model. Anesth Analg. 1995, 81: 1043-1051. 10.1097/00000539-199511000-00026.
Constable PD: Hyperchloremic acidosis: the classic example of strong ion acidosis. Anesth Analg. 2003, 96: 919-922. 10.1213/01.ANE.0000053256.77500.9D.
Alfaro V, Torras R, Ibanez J, Palacios L: A physical–chemical analysis of the acid-base response to chronic obstructive pulmonary disease. Can J Physiol Pharmacol. 1996, 74: 1229-1235. 10.1139/cjpp-74-11-1229.
Cameron JN: Acid–base homeostasis: past and present perspectives. Phys Zool. 1989, 62: 845-865.
About this article
Cite this article
Wooten, E.W. Science review: Quantitative acid–base physiology using the Stewart model. Crit Care 8, 448 (2004). https://doi.org/10.1186/cc2910
- base excess
- Stewart model