Potential pitfalls of using a correction to normative values for the assessment of acid–base compensation during early ascent to high altitude
Bibliographic record
Abstract
We read with interest the article by Skalla and colleagues recently published in Experimental Physiology (Skalla et al., 2025). The addition of data highlighting the early time course of respiratory and renal responses with ascent to a moderately high altitude contributes to our understanding of the early phases of integrated responses to acute exposure and acclimatisation in healthy humans with ascent. Many organ systems respond to the acute hypoxic stimulus with high altitude ascent, and indeed, the response of one organ system (e.g., the hypoxic ventilatory response and resulting hypocapnia and respiratory alkalosis) is the stimulus for another (e.g., renal compensation via bicarbonate diuresis). Assessing the early time course of exposure to high altitude is of interest to applied human physiologists and clinicians alike, due to the dynamic nature of integrated acute responses and acclimatisation, and potential implications for acute mountain sickness. Skalla et al. (2025) measured arterial blood gases and acid–base variables before and 24 and 44 h after arrival at 3100 m. We noted that despite an incremental reduction in [HCO3−]a, pH remained alkalaemic, even after 44 h at 3100 m. We were particularly pleased to see our ‘renal reactivity’ metric applied to early ascent, where each individual's renal response (Δ[HCO3−]a) was indexed against the renal stimulus, taken as Δ P aC O 2 ${P_{{\mathrm{aC}}{{\mathrm{O}}_2}}}$ , following ascent (Zouboules et al., 2018). We noted that Skalla et al. (2025) also assessed early blood acid–base acclimatisation using altitude-corrected extracellular fluid base excess (ecf-BE), mathematically equivalent to titratable hydrogen ion difference (THID), a correction initially introduced and intended for application in acclimatised and chronically dwelling high-altitude residents (e.g., Zubieta-Calleja et al., 2011). With ascent to high altitude and renal compensation, the Van Slyke equation results in a negative BE (or, in other words, a base deficit), congruent with reduced arterial [bicarbonate], as expected. However, the THID equation takes into account the altitude, normalising BE toward assumed baseline values. We have two main critiques regarding the use of this derivation of BE with early ascent, namely: (a) THID was not developed to assess early-acclimatisation responses, but rather developed as a steady-state tool to assess acid–base status with sustained high-altitude residence, after which time acid–base homeostasis has stabilized, and (b) the use of assumed normative values in these correction factors is potentially misleading, particularly when researchers have each participant's actual measured values available. THID (ecf-BE) is a steady-state normalisation technique originally intended to assess acid–base status in acclimatised high-altitude residents. It was introduced to prevent erroneous diagnosis of acid–base disturbances in high altitude residents whose ventilatory and renal responses or compensation established a new steady-state setpoint, in other words, a stable homeostatic target for P aC O 2 ${P_{{\mathrm{aC}}{{\mathrm{O}}_2}}}$ , [HCO3−]a and pHa under chronic exposure and residence in the context of hypobaric hypoxia (Zubieta-Calleja et al., 2011). Accordingly, the original Van Slyke equation was modified so that individuals could be assessed at this new altitude-acclimatised setpoint, allowing clinicians to distinguish a newly-imposed (i.e., clinical) acid–base disturbance from the resident's ‘normal’ altitude-compensated baseline (Zubieta-Calleja et al., 2011). Thus, THID was developed in high altitude resident populations (e.g., indigenous highlanders in Bolivia) and assumes population-level reference values for P aC O 2 ${P_{{\mathrm{aC}}{{\mathrm{O}}_2}}}$ , [HCO3−]a and pHa consistent with full (or near-full) renal compensation at that particular altitude (Zubieta-Calleja et al., 2011). This specific clinical application does not translate to lowlanders during the first 24–44 h of ascent, a non-steady-state period, during which lowlanders typically remain alkalaemic for days (e.g., Zouboules et al., 2018), if not indefinitely. Applying correction factors to visitors can (a) normalise away the very perturbation we seek to quantify with acute ascent, namely, acute respiratory alkalosis with early but incomplete renal compensation, and (b) invert or obscure the direction and magnitude of these responses. In this non-steady-state time domain, marked by substantial within- and between-individual variability (e.g., Zouboules et al., 2018), using a correction anchored to assumed ‘textbook’ values is unnecessary when each participant's P aC O 2 ${P_{{\mathrm{aC}}{{\mathrm{O}}_2}}}$ , [HCO3−]a and pHa are measured directly, and it risks distorting the interpretation of inter-individual variability in response magnitudes. When describing physiological responses in this early time domain of high-altitude ascent, the critical information is the within-individual trajectory (Δ P aC O 2 ${P_{{\mathrm{aC}}{{\mathrm{O}}_2}}}$ , Δ[HCO3−]a and ΔpHa), as Skalla et al. (2025) also address. THID/ecf-BE was designed for the clinical problem of misclassification of acid–base disorders in high-altitude residents who are persistently alkalaemic, despite reduced bicarbonate (with associated negative BE values), not for quantifying acute physiology in high-altitude sojourners, and accordingly, in our view, should not be used to infer rapid renal compensation over 24–44 h, as in Skalla et al. (2025). THID/ecf-BE and related corrections, like the original BE calculation using the Van Slyke equation, rely on population-level reference values for P aC O 2 ${P_{{\mathrm{aC}}{{\mathrm{O}}_2}}}$ , [HCO3−]a and pHa. In the first 24–44 h after ascent to high altitude, where both baseline (pre-ascent) and post-ascent response values are heterogeneous, applying a single assumed reference value can distort the signal when each participant's P aC O 2 ${P_{{\mathrm{aC}}{{\mathrm{O}}_2}}}$ , [HCO3−]a and pHa have been measured directly. In this early time domain, the variables of interest are the within-individual change from baseline (Δ P aC O 2 ${P_{{\mathrm{aC}}{{\mathrm{O}}_2}}}$ , Δ[HCO3−]a and ΔpHa), not proximity to any external average population reference set. Reference value choice also introduces biological bias. For example, baseline P aC O 2 ${P_{{\mathrm{aC}}{{\mathrm{O}}_2}}}$ and [HCO3−]a differ by sex and menopausal status, as premenopausal females typically exhibit reduced P aC O 2 ${P_{{\mathrm{aC}}{{\mathrm{O}}_2}}}$ and [HCO3−]a compared with postmenopausal females and males (e.g., Loeppky et al., 2001), so using a single assumed reference value to calculate a correction index will systematically bias the reporting of female responses, and can inflate apparent ‘compensation’. Population heterogeneity further limits the utility of generalised correction factors. Indigenous highlanders are not interchangeable; Andean and Sherpa groups demonstrate distinct acid–base profiles across comparable altitudes (e.g., Tymko et al., 2022), and recent data demonstrate lowlanders differ from Tibetan highlanders with ascent and acclimatisation (e.g., Johnson et al., 2025). Consequently, applying a population-specific, Bolivian-derived correction factor to non-Bolivian highland populations or lowlanders with early ascent assumes a uniformity that does not exist. In addition, we recently showed that Tibetan highlanders have differential P C O 2 ${P_{{\mathrm{C}}{{\mathrm{O}}_2}}}$ and [HCO3−] compared to lowlanders even at 1400 m (Bushfield et al., 2025), so there is appreciable variability within and between groups, even prior to ascent. We appreciate Skalla and colleagues’ careful staging and early time course measurements with ascent to ∼3100 m. The dataset is valuable for adding to our emerging understanding of integrated respiratory–renal responses and elucidating interactions between organ systems over the first 2 days of acclimatisation to moderately high altitude. We suggest centring the presentation of data on participant-level trajectories of the biochemical variables of interest, namely, P aC O 2 ${P_{{\mathrm{aC}}{{\mathrm{O}}_2}}}$ , [HCO3−]a and pHa, where the conclusions do not depend on a correction to a single normative reference set. We offer these suggestions on the potential pitfalls of correction factors and hope that this discussion sharpens physiological inferences drawn for this and other contributions to the literature. Both authors wrote, reviewed, edited and approved the final manuscript and agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. All persons designated as authors qualify for authorship, and all those who qualify for authorship are listed. We are grateful to Mr Jordan Bird for helpful and constructive conversations on this topic. None declared.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame distilled prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. Learned from the 10,348 direct Codex labels and 10,348 direct Gemma labels. Candidate is the union of thresholded teacher heads; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels or direct frontier model labels.
Codex and Gemma teacher scores by category
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.000 | 0.000 |
| Meta-epidemiology (narrow) | 0.000 | 0.000 |
| Meta-epidemiology (broad) | 0.000 | 0.000 |
| Bibliometrics | 0.000 | 0.000 |
| Science and technology studies | 0.000 | 0.000 |
| Scholarly communication | 0.000 | 0.000 |
| Open science | 0.000 | 0.000 |
| Research integrity | 0.000 | 0.000 |
| Insufficient payload (model declined to judge) | 0.000 | 0.000 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one teacher head, not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".