Document Detail

On the kinetics of anaerobic power.
Jump to Full Text
MedLine Citation:
PMID:  22830586     Owner:  NLM     Status:  MEDLINE    
BACKGROUND: This study investigated two different mathematical models for the kinetics of anaerobic power. Model 1 assumes that the work power is linear with the work rate, while Model 2 assumes a linear relationship between the alactic anaerobic power and the rate of change of the aerobic power. In order to test these models, a cross country skier ran with poles on a treadmill at different exercise intensities. The aerobic power, based on the measured oxygen uptake, was used as input to the models, whereas the simulated blood lactate concentration was compared with experimental results. Thereafter, the metabolic rate from phosphocreatine break down was calculated theoretically. Finally, the models were used to compare phosphocreatine break down during continuous and interval exercises.
RESULTS: Good similarity was found between experimental and simulated blood lactate concentration during steady state exercise intensities. The measured blood lactate concentrations were lower than simulated for intensities above the lactate threshold, but higher than simulated during recovery after high intensity exercise when the simulated lactate concentration was averaged over the whole lactate space. This fit was improved when the simulated lactate concentration was separated into two compartments; muscles + internal organs and blood. Model 2 gave a better behavior of alactic energy than Model 1 when compared against invasive measurements presented in the literature. During continuous exercise, Model 2 showed that the alactic energy storage decreased with time, whereas Model 1 showed a minimum value when steady state aerobic conditions were achieved. During interval exercise the two models showed similar patterns of alactic energy.
CONCLUSIONS: The current study provides useful insight on the kinetics of anaerobic power. Overall, our data indicate that blood lactate levels can be accurately modeled during steady state, and suggests a linear relationship between the alactic anaerobic power and the rate of change of the aerobic power.
John F Moxnes; Kjell Hausken; Øyvind Sandbakk
Related Documents :
7882126 - Adherence to pelvic floor muscle exercise and long-term effect on stress urinary incont...
2938306 - Noninvasive screening procedure for evaluation of bladder emptying function.
10491346 - Non-steroidal anti-inflammatory drugs and renal response to exercise: a comparison of i...
16217326 - Factors that distinguish continent from incontinent older adults with detrusor overacti...
20930016 - Comparison of 6-min walk test distance and cardiopulmonary exercise test performance in...
21318926 - Changing exercise through targeting affective or cognitive attitudes.
Publication Detail:
Type:  Journal Article     Date:  2012-07-25
Journal Detail:
Title:  Theoretical biology & medical modelling     Volume:  9     ISSN:  1742-4682     ISO Abbreviation:  Theor Biol Med Model     Publication Date:  2012  
Date Detail:
Created Date:  2013-01-21     Completed Date:  2013-06-24     Revised Date:  2013-07-12    
Medline Journal Info:
Nlm Unique ID:  101224383     Medline TA:  Theor Biol Med Model     Country:  England    
Other Details:
Languages:  eng     Pagination:  29     Citation Subset:  IM    
Department of Human Movement Science, Norwegian University of Science and Technology, 7491, Trondheim, Norway.
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Computer Simulation
Models, Biological*
Muscle, Skeletal / physiology*

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Theor Biol Med Model
Journal ID (iso-abbrev): Theor Biol Med Model
ISSN: 1742-4682
Publisher: BioMed Central
Article Information
Download PDF
Copyright ©2012 Moxnes et al.; licensee BioMed Central Ltd.
Received Day: 13 Month: 3 Year: 2012
Accepted Day: 20 Month: 6 Year: 2012
collection publication date: Year: 2012
Electronic publication date: Day: 25 Month: 7 Year: 2012
Volume: 9First Page: 29 Last Page: 29
PubMed Id: 22830586
ID: 3548746
Publisher Id: 1742-4682-9-29
DOI: 10.1186/1742-4682-9-29

On the kinetics of anaerobic power
John F Moxnes1 Email:
Kjell Hausken2 Email:
Øyvind Sandbakk3 Email:
1Department for Protection, Norwegian Defence Research Establishment, P.O. Box 25, 2007, Kjeller, Norway
2Faculty of Social Sciences, University of Stavanger, 4036, Stavanger, Norway
3Department of Human Movement Science, Norwegian University of Science and Technology, 7491, Trondheim, Norway


During human exercise, it is well known that intracellular adenosine triphosphate (ATP) can be produced aerobically in the mitochondria by oxidative phosphorylation, anaerobically due to glycolysis or glycogenolysis, or by the breakdown of phosphocreatine (PCr) into Creatine (Cr) in the Creatine Kinease (CK) reaction. With the aerobic energy as input, the current paper develops mathematical models that simulate the kinetics of anaerobic ATP production and power.

The rate of oxygen (O2) consumption can be set to the sum of a constant rate (resting rate of O2 consumption), a rate due to unloaded body movements and at a rate that is used to perform work [1]. For moderate constant work rates, the O2 consumption increases to a steady state condition. However, at the onset of exercise or due to a change in work rate, there is a rate of change of VO2 (or aerobic power) that is named as the VO2 kinetics. Pulmonary rate of O2 consumption (VO2p) has been used as a proxy for VO2. For moderate intensity exercise at constant work rate below the lactate threshold three distinct phases have been observed for VO2p. Phase I is the cardio dynamic phase, which represents the circulatory transit delay from muscles to lungs. Phase II is the mono exponential increase of VO2p. This phase reflects the adjustment of VO2 due to the use of active skeletal muscles. Phase III is the steady state phase of VO2p and VO2 during moderate exercise intensities [2,3]. For work rates associated with sustained acidosis, the mono-exponential component is slowed compared with lower intensities below the lactate threshold. In addition, a delayed slow component is superimposed. The slow component begins around 100-200 s into the exercise and can result in either a delayed sub-maximal steady state or a steady state equal to the maximal oxygen uptake (VO2 max) [4]. However, the mechanism of this slow component has not been resolved [5]. In this paper, we model phase II and III because only these two phases are considered relevant for the anaerobic alactic power model.

When exercise intensity increases and the rate of ATP production by oxidative sources becomes insufficient, anaerobic ATP production is required. When ATP is produced by glycolysis or glycogenolysis the endpoint is pyruvate, which can be reduced to lactate or oxidized to CO2 or H2O. The blood lactate concentration in the lactate pool is the result of the appearance of lactate from working muscles and various tissues and the disappearance of lactate in the skeletal muscles, the heart, the liver and the kidney cortex [6-9]. During steady state, lactate production (influx) is equal to lactate removal (outflux). As a result, the lactate concentration in the lactate pool stays constant, and the rate of oxygen consumption is the measure of the whole body energy expenditure regardless of the magnitude of lactate production and removal or the absolute blood lactate concentration [10]. The concept of a maximal lactate steady state can be defined as the highest level of intensity where a steady state condition of lactate can be obtained, which is also referred to as the lactate threshold. At exercise intensities above the lactate threshold the rise in the lactate concentration could be attributed to an increase in the rate of lactate appearance or the result of a decrease in the rate of lactate disappearance [10].

The maximum anaerobic energy that can be utilized is proportional to the sum of Cr and lactate that can accumulate in the body. PCr is an energy buffer that supports the transient failure of other metabolic pathways to support ATP production. The equilibrium constant of the CK reaction is around 20 and the slightest drop in ATP allows the reaction to proceed to ATP production [11]. Thus, the ATP concentration stays nearly constant until almost all the PCr is utilized. Rossiter et al. [12] found that the PCr levels follow an exponential time course after changes in work rate before approaching a steady state condition at moderate exercise intensities. In such cases, a strong similarity has been reported for the time constants of the VO2 kinetics and the PCr consumption [12]. Margaria [13] was the first to propose a hydraulic model for examining the whole body energy process during exercise. Despite this breakthrough, Margaria’s model was not quantified. Morton [14] presented a generalization of this model that was solved mathematically and compared with experimental data. Here Morton [14] modeled aerobic power and alactic anaerobic power with an exponential time devolvement without any anaerobic lactic power during steady state exercises below the lactate threshold. However, for exercise intensities that are above the lactate threshold, the anaerobic glycolytic energy supply is significant. The association between PCr and VO2 rate constants for exercises at such intensities has not yet been systematically reported. Furthermore, during recovery after high intensity exercise, the level of PCr must be restored, the pH must be re-established and ADP removed. While the PCr recovery is mainly due to oxidative ATP synthesis, the PCr stores may be rebuilt by anaerobic glycolysis [15-17]. Altogether, these insights provide a theoretical background for the models developed in this paper.

The O2-deficit formula of Medbø et al. [18] is an alternative model for the anaerobic power that accounts for lactic and alactic anaerobic power in which the chemical coupling efficiencies are assumed to be similar. Medbø et al. [18] suggested that the O2-deficit can be calculated by assuming that the metabolic power at intensities above VO2max can be estimated by extrapolating the steady state linear relationship between work rate and VO2 at submaximal intensities. The validity of the O2-deficit method has been widely debated [e.g., [19-21]. However, a rationale for the O2-deficit model is the assumption that the chemical coupling efficiencies of the three sources of ATP synthesis are similar.

The current study investigated two different mathematical models for the kinetics of anaerobic power during whole body exercise at different intensities. In order to test the models during exercise, oxygen uptake and blood lactate concentration were measured while a cross country skier ran with poles on a treadmill. Aerobic power, based on oxygen uptake measurements, was used as input. The lactic anaerobic power was calculated with a model presented by Moxnes and Hausken [10] using lactate concentration averaged over the whole lactate space, and with a model of Moxnes and Sandbakk [22] where the lactate concentration was separated into two compartments: muscles + internal organs and blood. The current study has input from these two previous studies, and focuses on alactic power as a novel contribution. Therefore, the power due to PCr break down was calculated theoretically and compared against the results of Jeneson et al. [23]. Finally, the models were used to compare PCr break down during continuous and interval exercise.

Overall design

The current study simulated the kinetics of anaerobic powers during whole body exercise at different exercise intensities through the use of Mathematica 8 (Wolfram Research Inc., Champaign, IL, USA). In a first model (Model 1), we suggested that for all exercise intensities, the work power would be linear to the work rate. As an alternative model (Model 2) where the lactic power and energy were the same as for Model 1, we hypothesized that the alactic anaerobic power would be proportional to the rate of change of the aerobic power [24,25]. For both models we used a first order differential equation for the time development of aerobic power as a function of work rate. As most parameters in the models needed to be fitted to each individual athlete, we fitted all parameters to a male Norwegian national level cross-country skier, with a body mass of m = 77.5 kg, body height of 181 cm and 600 h of endurance training per year. Thereafter, the simulations were compared with experimental data where this skier was running with poles on a treadmill (see details below).

The aerobic and anaerobic powers

ATP is produced by three different sources:

[Formula ID: bmcM1a]

[Formula ID: bmcM1b]

In these equations, “model” means model assumption. Qr is the resting metabolic power, set to 80 J/s based on oxygen uptake measurements of this skier during rest in our laboratory. QG is ATP production by glycolysis/glycogenolysis and CK by phosphocreatine break down to creatine. Qaw is aerobic power due to internal and external work. The aerobic power is Qa=Qr+Qaw.

For each ATP used (produced and consumed) heat and work are outputs. Heat when producing ATP differs between aerobic and anaerobic sources. However, heat when consuming ATP is the same for the aerobic and anaerobic sources. We defined work power and heat power due to oxidative phosphorylation (a), glycolysis/glycogenolysis (G) and phosphocreatine break down to creatine (CK) as:

[Formula ID: bmcM2a]

[Formula ID: bmcM2b]

[Formula ID: bmcM2c]

[Formula ID: bmcM2d]

[Formula ID: bmcM2e]

where “def” means definition. ηa, ηG and ηCK are the chemical coupling efficiencies when producing ATP aerobically, anaerobically by glycolysis/glycogenolysis (G), and anaerobically by CK, respectively. The two chemical coupling efficiencies ηG and ηa are similar, whereas ηCK is larger [26]. η is the mechanical efficiency, which means the work per unit use of energy of 1 ATP. PT includes both internal and external work.

At exercise intensities that exceed the maximal aerobic power Qmax, the virtual steady state is the steady state that would be attained if it was possible to carry out the exercise under purely aerobic conditions [27]. This virtual steady state is never reached as the aerobic power ends when Qmax is achieved. Figure 1 shows the steady state aerobic power Q¯a for an elite skier running with poles at the angle of inclination of 0.105 radians on a large treadmill. A linear relationship applies for Q¯a versus work rate up to approximately 95% of Qmax. When repeating the experiment with the same skier, the difference did not exceed 25 W for this skier which is considered an upper bound error margin.

We defined a steady state virtual power Q¯vir which is a straight line that crosses Qmax. Thus we set:

[Formula ID: bmcM3a]

[Formula ID: bmcM3b]

where c is a parameter. P is the work rate. Qun is the metabolic power due to unloaded running with poles. In general, Qun depends on the angle of inclination and the cycle frequency of the activity. For a given incline the frequency depends of the work rate. Thus Qun becomes a function of the work rate even for a fixed incline. P¯vir includes both internal and external work power. Q¯vir is input to the further modeling. We set that Q¯vir is a linear interpolation function through our data points. For values above or below the maximum and minimum measured Q¯vir we used linear extrapolation. During steady state, Qaw=Q¯a−Qr=MinQmax−Qr,Q¯vir−Qr.

The aerobic power is delayed with a time lag during steady state work rate. Qa(t) is the aerobic power and Qvir is the virtual aerobic power, Qat=modelMinQmax,Qvirt. To account for the delay mathematically, we used a first order differential equation of the virtual power, to read:

[Formula ID: bmcM4a]

[Formula ID: bmcM4b]

The “dot” above a variable means time derivative. τa is the so-called e-folding time, which is the time interval when an exponentially growing quantity increases by a factor of e. Thus τa is a time parameter that scales from the initiation of the activity until the aerobic power reaches a steady state asymptotically when assuming a constant work rate. In practice, reaching steady state means having less than 1% change in aerobic power per second. The cardio dynamic phase is not necessary to account for since only the true aerobic power is used as input for our anaerobic models. Typical values for τa are shown to be 10-36 s for moderate intensity exercise [22,27-33]. For example, di Prampero [27] suggested 10-24 s, Ceretelli et al. [32] found that τa increases linearly with the concentration of lactate up to 36 s, and Bizoni et al. [33] found τa = 23 s for all work rates. We set that τa = 30 s as a compromise. Equation (4b) ensures that the aerobic power is less than the maximal aerobic power.

If the virtual metabolic power is equal to or below Qmax we have that Qa=MinQmax,Qvir=Qvir. This gives from equation (4a and b) that:

[Formula ID: bmcM5]

For steady work rate equation (5) gives the analytical solution:

[Formula ID: bmcM6]

If the work rate is sufficiently low, Q¯a=Q¯vir. Notable is that solution (6) will only apply for a restricted time period unless Q¯a=Q¯vir.

The anaerobic power due to anaerobic glycolysis or glycogenolysis can be calculated from the total lactate concentration by using the equation from di Prampero and Ferretti [27], to read:

[Formula ID: bmcM7]

C is the lactate concentration defined as the amount of lactate per unit volume of lactate space (including muscles and blood). m is the mass of the body. λ converts the lactate per unit body mass to oxygen consumption. Note that during dynamic situations the blood lactate concentration differs significantly from the muscle lactate concentration. In this article, we measured the blood lactate concentration, which is different from the muscle lactate concentration. However, earlier research has demonstrated that the lactate concentration can be calculated as summarized in Appendix A [10].

The lactic energy used from time t0 to time t is:

[Formula ID: bmcM8]

Experimentally, EG(t0, t) can be found by measuring changes in the lactate concentration before and after exercise.

From (2) we have that:

[Formula ID: bmcM9]

The metabolic power is then given by:

[Formula ID: bmcM10]

Steady state is achieved at exercise intensities below the lactate threshold. We defined steady state with the lactate concentration regarded as steady (i.e., no lactic power) and the aerobic power was steady, to read Thus Qaw is a constant and Q−Qr=Qaw during steady state. Hence, it follows that QCK = 0. Thus the alactic power is zero. To achieve QCK = 0 during steady state in (9) we must have PT=ηηaQaw. However, during steady state (below LT) Qaw=Qun+cP. Thus, below the lactate threshold PT=Qun+cPηηa=P¯vir during steady state.

To close the equations in (9) and (10), which produce an indeterminate solution by themselves unless we have steady state, we needed a model for the work power PT in general. As Model 1 we forecasted that the work power PT was similar to the steady state virtual aerobic metabolic power. Thus, Model 1:

[Formula ID: bmcM11]

Inserting (11) into (9)-(10) gives the alactic power as:

[Formula ID: bmcM12a]

and the metabolic rate as:

[Formula ID: bmcM12b]

To calculate QCK1t for a work rate P(t) in (12a) we inserted Q¯vir from (3), Qa(t) from (4) and QG(t) from (7). It should be noted that the mechanical efficiency η is absent in (12a and b). We conceived that the work power PT was proportional to the rate of consumption of ATP. Thus, an equivalent model to (11) would be that the rate of ATP consumption is linear with the work rate.

The alactic energy is for Model 1:

[Formula ID: bmcM13]

The anerobic power and the anaerobic energy used from time t0 to time t for Model 1 are:

[Formula ID: bmcM14a]


[Formula ID: bmcM14b]

where EDF(t0, t) is in the literature named the oxygen deficit of the exercise. Medbø et al. [18] and Losnegard et al. [29] calculated the anaerobic energy during exercise as EDF(t0, t) (assuming that the rate of change of O2 consumption is proportional to the aerobic power). Equation (14b) shows that EDFt0,t=Eant0,t only if ηG=ηCK=ηa. Thus, we deduced that the O2-deficit model for the anaerobic energy was equal to (14b) if the chemical coupling efficiencies were alike.

Model 1 is a general version of the O2-deficit model. As the O2-deficit model has been widely debated, we proposed an alternative in Model 2 where the lactic power and energy were the same as for Model 1. However, we hypothesized that the alactic power was proportional to the rate of change of the aerobic power [24,25]. Thus, for Model 2:

[Formula ID: bmcM15]

where θ is a constant of proportionality parameter that we determine below. Equation (15) gives that:

[Formula ID: bmcM16]

If the virtual metabolic power is below Qmax and Q¯a=Q¯vir, we achieved from equations (5) and (16) that:

[Formula ID: bmcM17]

The alactic energy can be found by calculating the area between Q¯a and Qa in a power time diagram. This area is multiplied with θ/τa to find the alactic energy used. θ/τa can be considered as the effectiveness of alactic ATP production relative to the aerobic ATP production. Thus we set θ/τa=ηa/ηCK.

The mechanical efficiency was assumed to be η ≈ 0.5. For the chemical efficiency related to aerobic or lactic power ηG=ηaarebetween0:3and0:7. We used ηG=ηa≈0.6. This gave η×ηG=η×ηa=0.3. Since c was around 6.6 we found that ηr=1/c×η×ηG=1/6.6×0.5×0.6=0.51. Gonzales-Alonso et al. [26] concluded from experimental data that the heat per use of ATP was around two times larger for oxidative phosphorylation and anaerobic glycolysis compared to ATP from CK. For the alactic power we therefore forecasted that ηCK = 0.95. This gave that θ=ηa/ηCKτa=20s.

It has been shown that during intensities above the lactate threshold, ATP utilization increases and mechanical efficiency decreases at constant work rates. This may be explained by a change in fiber type recruitment, an elevated temperature, lowered pH or increased Pi levels [30]. We forecasted that τa, η, ηa, ηG and ηCK mainly depend on the lactate concentration and change in muscle pH, and assumed that a lower mechanical efficiency of muscle contractions and a lower P/O2 ratio could explain the drift in the steady state ATP consumption for a given work rate. To account for time varying parameters see the model in Appendix B.

Increased fitness after exercise is developed during the recovery period after the exercise. We defined an exercise which starts at time t0 and ends at time t1, and a subsequent recovery period which starts at time t1 and ends at time t2. Time t2 is defined as Qat2=Qat0 and C(t0) = C(t2). This means that EGt0,t2=0 and EGt0,t1=−EGt1,t2. We defined −EDF(t1, t2) as the energy depth EDB(t1, t2), to read:

[Formula ID: bmcM18]

We also defined the so-called excess post–exercise oxygen consumption (EPOC). In general EDB(t1, t2) is different from EPOC(t1, t2) unless Q¯vir=Q¯a. However, the latter is usually the case during recovery.

For Model 1 in equations (13) we achieved:

[Formula ID: bmcM19]

In this case the alactic energy used from time t0 to time t1 for Model 1 was therefore given by:

[Formula ID: bmcM20]

For Model 2 we have:

[Formula ID: bmcM21]

We defined EPOCAlt as an alternative variant of the excess post-exercise oxygen consumption:

[Formula ID: bmcM22]

Thus, finally:

[Formula ID: bmcM23]

Here, the alactic energy used from time t0 to time t1, could be found by measuring the aerobic power at time t1 and t2. The respective aerobic powers could be subtracted and then multiplied with θ to find the alactic energy.

It was also of interest to compare Model 1 and Model 2 more directly when Qvir = Qa and Q¯vir=Q¯a . When using equation (5) this gave:

[Formula ID: bmcM24]
IfQvir=QaandQ¯vir=Q¯a,then EPOCAltt1,t2=−τa∫t1t2Q·at′dt′=−∫t1t2Qat′−Q¯virdt′=EDBt1,t2=−∫t1t2Qat′−Q¯adt′=EPOCt1,t2

From equations (20)-(23) this gave:

[Formula ID: bmcM25a]

[Formula ID: bmcM25b]

These two methods gave different results due to fact that the second term for ECK1t0,t1 lacks in Model 2.

Hill et al. [34-37] have explained that the oxidative removal of lactate increases aerobic power. Using equation (25a) we achieved that:

[Formula ID: bmcM26]

If ηG = ηa, EPOC(t1, t2) corresponds to the total anaerobic energy, i.e. lactic and alactic anaerobic energy. Indeed, Margaria et al. [38] later on modified the concept of Hill et al. [34-37], and suggested that the increased aerobic power consisted of the fast alactic component and the slower lactic component. Finally, Gaesser and Brooks [39] introduced the term “excess post-exercise oxygen consumption”, which also included the much more prolonged increase in aerobic power that is observed for hours after exercise. Model 1 did not account for this very slow component that lasts for hours. Model 2 in equation (25b) said that EPOC corresponded to the alactic component. However, the development of the energy in equation (25b) depended of the time dynamics for the aerobic power modeled according to equation (4).

Experimental tests

The derived mathematical simulations were compared with experimental data from an elite cross-country skier while running with poles on a treadmill. The mass of the skier was m = 77.5 kg in all tests. All treadmill tests were performed on a 6 × 3 m motor-driven treadmill (Bonte Technology, Zwolle, The Netherlands). Inclination and speed were calibrated using the Qualisys Pro Reflex system and the Qualisys Track Manager software (Qualisys AB, Gothenburg, Sweden). The treadmill belt consisted of a non-slip rubber surface that allowed the skier to use his own poles (pole length: 80% of body height) with special carbide tips. Gas exchange values were measured by open-circuit indirect calorimetry using an Oxycon Pro apparatus with a mixing chamber (Jaeger GmbH, Hoechberg, Germany). Before each measurement, the VO2 and VCO2 gas analyzers were calibrated using high-precision gases (16.00 ± 0.04% O2 and 5.00 ± 0.1% CO2, Riessner-Gase GmbH & co, Lichtenfels, Germany). The inspiratory flow meter was calibrated with a 3 L volume syringe (Hans Rudolph Inc., Kansas City, MO). Heart rate (HR) was measured with a heart rate monitor (Polar S610, Polar Electro OY, Kempele, Finland), using a 5 s interval for data storage. Blood lactate concentration (BLa) was measured on 5 μL samples taken from the fingertip by a Lactate Pro LT-1710 t (ArkRay Inc, Kyoto, Japan).

Eight experimental protocols were performed, each on separate days with a minimum of 48 h between. The order of tests was performed as presented below. In order to investigate whole body exercise, running with poles was used in all tests. Before each testing session a standardized, test-specific 20-min warm-up was performed. Training on the days before testing was standardized, and the subject drank a standard fluid with sugar and electrolytes during all breaks while testing.

On the first test day, the skier performed six 5 min bouts with constant work rates at 0.105 inclines in radians. Five speeds at 0.25 m/s intervals below the lactate threshold were chosen, starting at 7.5 km/h = 2.08 m/s, followed by 2.33 m/s, 2.58 m/s, 2.83 m/s and 3.08 m/s. The sixth speed was increased by 0.125 m/s, giving 3.19 m/s, which is slightly above the lactate threshold. The starting speed was chosen based on experience from earlier tests of this athlete. 5-min bouts were used to obtain steady state conditions. A 2-min break with low-intensity walking was mandatory between each of the exercise bouts. Gas exchange values and heart rates were determined by the average of the last minute during each stage and blood lactate was measured directly after finishing each stage. The reason for using the last minute to assess respiratory variables was that the athletes are able to keep a more steady technique and physiological responses after 3-4 min.

In follow-up tests, we aimed to investigate whether our simulations fitted experimental data when varying the exercise time and during recovery after high-intensity exercise. Also, velocities slightly above or below the lactate threshold were tested in order to test the sensitivity of our model for lactate. Thus, the skier performed the following six tests on separate days with at least 2 days in between each test, at 0.105 radians incline with blood lactate measured directly after each stage:

– 2.87 m/s was performed for 800 s.

– 3.08 m/s was performed over 800 s, followed by a 5-min recovery phase at 1.67 m/s with blood lactate measured after 2 and 5 min. Additionally, a 3.05 m/s stage was performed to exhaustion.

– 3.08 m/s was performed over 2000 s.

– 3.19 m/s stages were performed for 400 s and to exhaustion.

– A 3.33 m/s stage was performed for 400 s.

– 3.88 m/s stages were performed for 150 s, 200 s and to exhaustion.

Finally, maximal metabolic power was tested on a separate day at an incline of 0.105 radians with a starting speed of 3 m/s. The speed was increased by 0.3 m/s every minute until exhaustion. VO2 was measured continuously, and the average of the three highest 10 s consecutive measurements determined VO2max and was used to calculate the maximal metabolic power. The test was considered to be a maximal effort if the following three criteria were met: 1) a plateau in VO2 was obtained with increasing exercise intensity, 2) respiratory exchange ratio above 1.10, and 3) blood lactate concentration exceeding 8 mmol/L.


Figure 2 illustrates that the experimental steady state blood lactate values showed good similarity with the simulated values based on Moxnes and Hausken’s [10] one compartment model for the average concentration in the total lactate pool. The blood lactate concentration equals the average concentration during steady state. At high lactate concentrations, the drain of lactate is saturated, which can be modeled by the function TanhχCt. No saturation was achieved when χ→0, since Limχ→0TanhχCt/χ=Ct. The lactate threshold QLT was achieved when 1−χp0d0Qmax−Qa=0⇒Qa=QLT=Qmax−χp0d0,QLTQmax=1−χp0d0Qmax. Qmax was known and we used the numerical value of p0 from Moxnes and Hausken [10]. Thus, we only needed to fit the two parameters χ/d0 and χ to the steady state measured values.

In the first step, we applied visual curve fitting, which means that we chose plausible values for the two parameters. We then plotted and compared the solution visually with the experimental data. The values of the two parameters were changed repeatedly until a good visual fit was obtained, while ensuring that the parameters had physiologically trustworthy numerical values. In the second step, a least square fit to the data was performed to produce ex post best fit estimates of these two parameters using the visual estimates as starting guess points and choosing a range around each starting point of the parameters. The method was performed separately for each parameter keeping the other parameter fixed. Steps 1 and 2 were repeated sufficiently until we were certain that we had obtained the optimal values for each of the two parameters. For robustness in the calculations, we applied the steady state solution. If we applied that Cs≈1+p0d0Qmax1−Qa/Qmax, the parameter χ got removed. The lactate threshold became Qa = Qmax. χ determined the lactate threshold by the exact relation. The rest of the lactate curve was dependent on p0d0Qmax1−Qa/Qmax. This gave a linear dependency of the lactate curve on p0/d0 and an inverse dependency of the lactate curve on QmaxQa.

Figure 3 demonstrates the measured and simulated lactate concentration during steady state work rate. In this simulation, we only applied visual fitting to determine the d0 parameter, and found that the simulated lactate concentration was somewhat higher than the measured (blood lactate) concentration for intensities above the lactate threshold. During recovery, the measured blood lactate concentration was higher than the simulated lactate concentration. Thus, we also calculated values according to Moxnes and Sandbakk [22] where the lactate concentration was separated into two compartments: muscles + internal organs and blood. In that case, Figure 3 shows that the simulation results are more similar to the experimental results. However, some discrepancies still appear for the higher work rates.

Figure 4 shows the different anaerobic alactic powers for intensities below the lactate threshold when using Models 1 and 2. The calculated alactic power for Model 1 showed a somewhat complicated behavior since it also depended on the model for the lactic power. A special feature was the overshoot that applied during recovery because the lactic power becomes negative. For Model 2 the lactic power was only different from zero when the rate of change of the aerobic power was different from zero.

Figure 5 illustrates the lactic and alactic storage of energy. We found different behaviors for Models 1 and 2. For Model 1, the alactic energy storage started to increase from a local minimum at the time the aerobic power reached a steady state. The aerobic power reached steady state before the lactic power did for v4 and v5. When the exercise was terminated for v6, v7 and v8, the alactic energy was fully restored. The alactic energy resources increased due to use of lactic power. Although the details of the time history of alactic energy depend on the chemical coupling coefficients, the overall scenario was much the same for different numerical values of the coefficients. It is important to note that the alactic power depended on the model for the lactic power. When O2 was excluded, the PCr store was rebuilt by anaerobic glycogenolysis or glycolysis [15-17]. In Model 2, the alactic energy showed a simpler behavior than in Model 1 since the alactic power only depended on the rate of change of the aerobic power. When comparing the solution in Figure 5 in the current paper with the invasive measurements shown in Figure 2 in the paper Jeneson et al. [23] it appeared that Model 2 gave a better fit than Model 1.

To study the kinetics of anaerobic powers more in detail we compared interval exercise with continuous exercise by using Models 1 and 2. These were simulated exercises matched for energy expended. Interval exercise changed in intensity with durations above and below the lactate threshold as described below. The continuous exercise maintained steady state intensity. The two exercises were given by:


T is the time variable, 0<=t<=T, and T = 360 s is a parameter for the time duration, i.e. how long the exercise lasted. Figure 6 shows the metabolic powers for interval and continuous exercises as functions of time, whereas Figure 7 shows the simulated lactate concentration for the two different exercises. The lactate values were reduced to the initial levels between intervals, whereas lactate continued to rise during continuous exercise. Figure 8 shows the lactic power and the alactic powers calculated by Models 1 and 2. Overall, the alactic energy stores were rebuilt during recovery with interval exercise, which allowed for the alactic energy to be used at the beginning of each interval. Figure 9 shows the lactic and alactic energy utilization and shows that during exercise the alactic energy storage decreased with time for Model 2. For Model 1 a local minimum appeared around the time when the aerobic power reached steady state. During interval exercise, the alactic energy storage as a function of time did not differ much between the two models.


The current study investigated two different mathematical models for the kinetics of anaerobic power during whole body exercise at different exercise intensities. The results indicate that blood lactate levels can be accurately modeled during steady state, and suggest a linear relationship between the alactic anaerobic power and the rate of change of the aerobic power. Overall, we propose that the current simulation models provide useful insight into how the anaerobic powers during practical training and testing should be interpreted.

When comparing the experimental blood lactate values with Moxnes and Hausken’s [10] model for lactic power, we found similar results during steady state exercise intensities. Uncertainties of less than ± 0.5 mmol/L using these mathematical models are indicated. For intensities above the lactate threshold, the measured blood lactate values were significantly lower than those calculated by the Moxnes and Hausken [10] model, whereas during recovery after high intensity exercise, the blood lactate concentration was higher than calculated. This could be explained by a possible delay in blood lactate concentration compared to the average concentration in the total lactate pool that applies theoretically. Therefore, we compared the experimental results with a two-compartment model by Moxnes and Sandbakk [22], where lactate concentration was separated into muscle and blood compartments. The employment of this model gave a better fit between simulated and experimental results, especially for the highest work rates and during recovery. However, some discrepancies still appeared for the higher work rates, which is a topic for further research.

The calculation of alactic storage of energy in Model 1, which is a generalization of the O2-deficit model, showed that the alactic energy storage increased due to use of lactic power from the time the aerobic power reached a steady state. However, the various inferences about the alactic storage of energy during steady state work rate have some shortcomings. For example, the results may have been influenced by the fact that the O2-deficit is composed of both lactic and alactic metabolisms. Thus, no definite conclusion could be drawn for Model 1 for alactic power. For Model 2, the alactic energy showed a simpler behavior since it only depended on the rate of change of the aerobic power. When comparing Models 1 and 2 with the invasive measurements by Jeneson et al. [23], we concluded that Model 2 gave the best behavior. Thus, a linear relationship between the alactic anaerobic power and the rate of change of the aerobic power was suggested. In future research, it remains a challenge to confirm the contribution of lactic and alactic anaerobic metabolism and power.

The simulations of interval exercises showed a drop in lactate levels between the intervals, and that the alactic energy stores were rebuilt during recovery. Practically, this allows athletes to work at high exercise intensities for longer time durations without lowering the pH values in the muscles compared to continuous exercise at similar intensities. This might be one reason for the effectiveness of interval exercise, as verified both by the scientific literature [40] and by elite athletes’ practical training programs where interval exercises are an important component. When the alactic anaerobic powers during continuous exercises from Models 1 and 2 were compared, Model 1 showed a local minimum appearing around the time when the aerobic power reached steady state for continuous exercise, whereas for Model 2 the alactic energy storage decreased with time. During interval exercise, the alactic energy storage as a function of time was similar for the two models.

Modeling in human biology is a challenge since one is confronted with conceiving a simple, but realistic representation of complex phenomena. The parameter values used in the current study will be dependent on the fitness level of the individual being tested and on the concentration of glycogen in the muscles. To construct valid simulation models during dynamic exercise such models need to be developed for each individual.

Appendix A: The lactate model

In equation (25) of Moxnes and Hausken [10] we find that:

[Formula ID: bmcMA1]

where p1() and d1() are dimensionless functions and α is a parameter to describe the saturation of the concentration of lactate C(t). We use TanhχCtχ instead of 1−Exp−αCtα to describe saturation. To read:

[Formula ID: bmcMA2]

We set that:

[Formula ID: bmcMA3]

and β0 is a parameter that must be fitted to the experimental data, as earlier done by Moxnes and Sandbakk [22] when studying a two compartment model of lactate. The steady state solution is given by:

[Formula ID: bmcMA4]

We use 1/β0p0=0.6Qmax, where 0.6Qmax is around peak fat metabolism. All together we set:

[Formula ID: bmcMA5]

Appendix B: Model for parameters varying with the lactate level

To account for situations with varying τa, ηa, ηG and ηCK we can use a first order differential equation of the virtual aerobic metabolic power, to read:

[Formula ID: bmcMB1]

Equation (B1) allows the calculation of the aerobic power when the efficiencies and the time parameter τa are time dependent. This might appear due to changes in lactate levels. We assume that η × ηa is constant through time. This gives from (B1):

[Formula ID: bmcMB2]

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JM performed the mathematical simulations, ØS performed all laboratory testing and KH supported in finishing the paper in its final form. All authors contributed with important intellectual content in all parts of the manuscript. All authors read and approved the final manuscript.

Cabrera ME,Saidel GM,Kalhan SC,Lactate metabolism during exercise: analysis by an integrative systems modelAm J PhysiolYear:  277R1522R1536
Stirling JR,Zakynthinaki MS,Saltin B,A model of oxygen uptake kinetics in response to exercise: Including a means of calculating oxygen demand/deficit/debtBull Math BiolYear: 200567989101510.1016/j.bulm.2004.12.00515998492
Murias JM,Spencer MD,Kowalchuk JM,Paterson DH,Influence of phase I duration on phase II VO2 kinetics parameters estimates in older and young adultsAm J Physiol Regul Integr Comp PhysiolYear: 2011301R218R22410.1152/ajpregu.00060.201121490368
Billat VL,Bocquet V,Slawinski J,Laffitte L,Demarle A,Chassaing P,Koralsztein JP,Effects of prior intermitted runs av v VOmax on oxygen kinetics during an all-out severe run in humansJ Sports Med Phys FitnessYear: 20004018519411125760
Gaesser GA,Poole DC,The slow component of oxygen uptake in humansExerc Sport Sci RevYear: 19962435708744246
Brooks GA,Lactate production under fully aerobic conditions: the lactate shuttle during rest and exerciseFed ProcYear: 198545292429293536591
Brooks GA,Current concepts in lactate exchangeMed Sci Sports ExercYear: 1991238959061956262
Brooks GA,Lactate shuttles in natureBiochem Soc TransYear: 20023025826412023861
Brooks GA,Link between glycolytic and oxidative metabolismSports MedYear: 20073734134310.2165/00007256-200737040-0001717465603
Moxnes JF,Hausken K,A mathematical model for the training impulse and lactate influx and outflux during exerciseInt J Mod Phys CYear: 200920114717710.1142/S0129183109013522
Carlson FD,Siger A,The mechano-chemistry of muscular contraction I. The isimetric twitchJ Gen PhysiolYear: 196044336010.1085/jgp.44.1.3313690828
Rossiter HB,Ward SA,Doyle VL,Howe FA,Griffiths JR,Whipp BJ,Interference from pulmonary O2 uptake with respect to intramuscular (phosphcreatine) kinetcs during moderate exercise in humansJ PhysiolYear: 199951892193210.1111/j.1469-7793.1999.0921p.x10421675
Margaria R,Biomechanics and energetic of muscular exerciseYear: 1976Oxford: Oxford University Press
Morton RH,A three component model of human bioenergeticsJ Math BiolYear: 19862445146610.1007/BF012368923805905
Crowther GJ,Kemper WF,Carey MF,Conley KE,Control of glycolysis in contracting skeletal muscle. II. Turning it offAm J Physiol Endocrinol MetabYear: 20022827479
Jubrias SA,Esselman PC,Price LB,Cress ME,Conley KE,Large energetic adaptations of elderly muscle to resistance and endurance trainingJ Appl PhysiolYear: 2001901663167010.1063/1.137445211299253
Lanza IR,Wigmore DM,Befroy DE,Kent-Braun JA,In vivo ATP production during free-flow and ischaemic muscle contractions in humansJ PhysiolYear: 200657735336710.1113/jphysiol.2006.11424916945975
Medbø JI,Mohn AC,Tabata I,Bahr R,Vaage O,Sejersted OM,Anaerobic capacity determined by maximum accumulated O2 deficitJ Appl PhysiolYear: 19886450603356666
Bangsbo J,Quantification of anaerobic energy production during intense exerciseMed Sci Sports ExercYear: 19983047529475643
Gastin PB,Energy system interaction and relative contribution during maximal exerciseSports MedYear: 20013172574110.2165/00007256-200131100-0000311547894
Noordhof DA,de Kroning JJ,Foster C,The maximal accumulated oxygen deficit method: a valid and reliable measure of anaerobic capacity?Sports MedYear: 20104028530210.2165/11530390-000000000-0000020364874
Moxnes JF,Sandbakk K,The kinetics of lactate production and removal during whole-body exerciseTheor Biol Med ModelYear: 20129710.1186/1742-4682-9-722413898
Jeneson JAL,Schmitz JPJ,van den Broek NMA,van Riel AAW,Hilbers PAJ,Nicolay K,Prompers JJ,Magnitude and control of mitocondrial sensitivity to ADPAm J Physiol Endocrinol MetabolYear: 2009297E774E78410.1152/ajpendo.00370.2009
Mahler M,First order kinetics of muscle oxygen consumption, and an equivalent proportionality between Qo2 and Phosphorylcreatine levelJ Gen PhysiolYear: 19858613516510.1085/jgp.86.1.1354031824
Meyer RA,A linear model of muscle respiration explains monoexponential phosphocreatine changesJ Appl PhysiolYear: 198875648656
Gonzalez-Alonso J,Quistorff B,Krustrup P,Bangsbo J,Saltin B,Heat production in human skeletal muscle at the onset of intense dynamic exerciseJ PhysiolYear: 2000524260361510.1111/j.1469-7793.2000.00603.x10766936
di Prampero PE,Ferretti G,The energetics of anaerobic muscle metabolism: a reappraisal of older and recent conceptsRespir PhysiolYear: 19991182–310311510647856
di Prampero PE,Factors limiting maximal performance in humansEur J PhysiolYear: 20039042042910.1007/s00421-003-0926-z
Losnegard T,Myklebust H,Hallen J,Energy system contribution as a determinant of performance in elite skiersYear: 201110.1249/MSS.0b013e3182388684
Bangsbo J,Krustrup P,Gonzalez-Alonso J,Saltin B,ATP production and efficiency of human skeletal muscle during intense exercise: effects of previous exerciseAm J Physiol Endocrinol MetabolYear: 2001280E956E964
Bertuzzi RC,Franchini E,Ugrinowitsch C,Kokubun E,Lima-Silva AE,Pires FO,Nakamura FY,Kiss MA,Predicting MAOD using only a supramaximal exhaustive testInt J Sports MedYear: 201031747748110.1055/s-0030-125337520432195
Cerretelli P,Pendergast DR,Paganelli WC,Rennie DW,Effects of specific muscle training on VO2-on response and early blood lactateJ Appl PhysiolYear: 197947761769511683
Binzoni T,Ferretti G,Sechenker K,Cerretelli P,Phosphocreatine hydrolysis in 31P-NMR at the onset of constant-load exercise in humansJ Appl PhysiolYear: 199273164416491447115
Hill VH,Lupton H,Muscular exercise, lactic acid, and the supply and utilization of oxygenQ J MedYear: 19231613517110.1093/qjmed/os-16.62.135
Hill AV,Long CNH,Lupton H,Muscular exercise, lactic acid, and the supply and utilization of oxygen: parts I-IIIProc Roy SocYear: 19249643847510.1098/rspb.1924.0037
Hill AV,Long CNH,Lupton H,Muscular exercise, lactic acid, and the supply and utilization of oxygen: parts IV-VIProc Roy SocYear: 1924978413810.1098/rspb.1924.0045
Hill AV,Long CNH,Lupton H,Muscular exercise, lactic acid, and the supply and utilization of oxygen: parts VII-VIIIProc Roy SocYear: 19249715517610.1098/rspb.1924.0048
Margaria R,Edwards HT,Dill OB,The possible mechanisms of contracting and paying the oxygen dept and the role of lactic acid in muscular contractionAm J PhysiolYear: 1933106689715
Gaesser GA,Brooks GA,Metabolic bases of excess post-exercise oxygen consumption: a reviewMed Sci Sports ExercYear: 19841629436369064
Laursen PB,Jenkins DG,The scientific basis for high-intensity interval training: Optimising training programmes and maximizing performance in highly trained endurance athletesSports MedYear: 200232537310.2165/00007256-200232010-0000311772161


[Figure ID: F1]
Figure 1 

The steady state aerobic powerQ¯aand the virtual steady stateQ˜¯in J/s as functions of work rate P in J/s when running with poles at the angle of inclination of 0.105 radians for 5-min steady state work rates. The lower horizontal line is the lactate threshold = 1654 J/s, and the upper horizontal line the maximal aerobic power Qmax=1886J/s. ■: Experimental data, v1=2.08m/s,Q¯vir=0.59Qmax, v2=2.33m/s,Q¯vir=0.66Qmax, v3=2.64m/s,Q¯vir=0.74Qmax, v4=2.83m/s,Q¯vir=0.80Qmax, v5=3.08m/s,Q¯vir=0.87Qmax, v6=3.19m/s,Q¯vir=0.90Qmax.

[Figure ID: F2]
Figure 2 

Steady state lactate concentration C in mmol/L as a function of the fraction of the maximum aerobic power Qmaxwhen running with poles at the angle of inclination of 0.105 for 5-min steady state work rates.v1=2.08m/s,Q¯vir=0.59Qmax, v2=2.33m/s,Q¯vir=0.66Qmax, v3=2.64m/s,Q¯vir=0.74Qmax, v4=2.83m/s,Q¯vir=0.80Qmax, v5=3.08m/s,Q¯vir=0.87Qmax, C=Cs=12Ln1+p0χd0Qmax1−Qa/Qmax1−p0χd0Qmax1−Qa/Qmax/χ, χ=2.54m3/kg,p0=10−5kg/m3s/J/s,d0=1.110−7/J/s2/s,Qmax=1886J/s,QLT/Qmax=0.88.

[Figure ID: F3]
Figure 3 

The theoretical and experimental blood lactate concentration C in mmol/L as a function of time t in seconds when running with poles at the angle of inclination of 0.105 radians. The eight symbols from below are experimental values from the different tests v1-8. The straight curves represent results from a one-compartment model by Moxnes and Hausken [2], whereas the three stapled curves represent results from the two compartment model by Moxnes and Sandbakk [25] for v5, v6 and v8. Lines from below: black =v1=2.08m/s,Q¯vir=0.59Qmax, orange =v2=2.33m/s,Q¯vir=0.66Qmax, cyan =v3=2.64m/s,Q¯vir=0.74Qmax, red =v4=2.83m/s,Q¯vir=0.80Qmax, blue =v5=3.08m/s,Q¯vir=0.87Qmax until a recovery phase starts at t = 800s with v = 1.67m/s, Q¯vir=0.47Qmax, green =v6=3.19m/s,Q¯vir=0.90Qmax, yellow =v7=3.33m/s,Q¯vir=0.94Qmax, pink =v8=3.88m/s,Q¯vir=1.09Qmax.LT=0.88Qmax.

[Figure ID: F4]
Figure 4 

Powers Q in J/s as function of time t in seconds when running with poles at the angle of inclination of 0.105 radians.v = 3.08 m/s, Q¯vir=0.87Qmax when t < 800 s. v = 1.67 m/s, Q¯vir=0.47Qmax when t > 800 s. ____ ___: Aerobic power Qa=Qb+Qaw, ________: Lactic anaerobic power (QG), -----------: Alactic anaerobic power (QCK) for Model 1, .........: Alactic anaerobic power (QCK) for Model 2.

[Figure ID: F5]
Figure 5 

The lactic and alactic anerobic energy utilization E(t)/E(0) as a function of time t in seconds for different work rates when running with poles at the angle of inclination of 0.105. ______: Lactic energy utilization EGt/EG0,EG0=870J/kg×m, - - - - -: Alactic energy utilization ECKt/ECK0,ECK0=830J/kg×m using Model 1, ......: Alactic energy utilization ECKt/ECK0,ECK0=830J/kg×m using Model 2. From below in different colors: pink =v8=3.88m/s,Q¯vir=1.09Qmaxyellow =v7=3.33m/s,Q¯vir=0.94Qmax, green =v6=3.19m/s,Q¯vir=0.90Qmax, blue =v5=3.08m/s,Q¯vir=0.87Qmax, black =v5=3.08m/s,Q¯vir=0.87Qmax until a recovery phase starts at t = 800s with v = 1.67m/s, Q¯vir=0.47Qmax, red =v4=2.83m/s,Q¯vir=0.80Qmax, The exercises are terminated for v8, v7 and v6.

[Figure ID: F6]
Figure 6 

The metabolic power as a fraction of the maximum aerobic power Q/Qmax as a function of time t in seconds for an interval and a continuous exercise. The horizontal lines represent the lactate threshold. -------: Continuous. _____: Interval.

[Figure ID: F7]
Figure 7 

The lactate concentration C in mmol/L as a function of time t in seconds for interval and continuous exercise during the protocol showed in Figure6. ....: Continuous. -------: Interval.

[Figure ID: F8]
Figure 8 

The anaerobic lactic and alactic powers (Q) in J/s as a function of time t in seconds for interval (black lines) and continuous (red lines) exercise during the protocol showed in Figure6. ______: Lactic anaerobic power (QG). ---------: Alactic anaerobic power (QCK), Model 1. ......: Alactic anaerobic power (QCK), Model 2.

[Figure ID: F9]
Figure 9 

Calculations of the energy utilization E(t)/E(0) as a function of time t in seconds for (black lines) and continuous (red lines) exercises during the protocol showed in Figure6. _______: Lactic anaerobic energy utilization EGt/EG0,EG(0)=870J/kg×m. - - - - - -: Alactic anaerobic energy utilization ECKt/ECK0,ECK0=830J/kg×m Model 1. ......: Alactic anaerobic energy utilization ECKt/ECK0,ECK0=830J/kg×m Model 2.

Article Categories:
  • Research

Keywords: Blood lactate, Cross-country skier, Mathematical model, Oxygen uptake, Phosphocreatine.

Previous Document:  Solvent-Dependent Gating Motions of an Extremophilic Lipase from Pseudomonas aeruginosa.
Next Document:  Breaking Bad News in a Southeast European Population: A Survey among Cancer Patients in Albania.